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ABSTRACT 


Approximate solutions are presented for the linear 
isothermal response to symmetric dynamic loading of viscously 
damped thin homogeneous isotropic elastic plates and shallow 
spherical shells having regular polygonal boundaries. Normal 
modes of vibration are determined with the aid of the boundary 
collocation technique, and modal participation functions, 
necessary for forced solutions, are calculated by employing 
numerical integration. The solutions are applicable to 
elastically built-in plates and shells with time-dependent 
edge conditions. Numerical solutions are provided for a 
large number of polygonal plates and shells subjected to 
uniform and central point loads for a variety of edge 
Peng tions... Lbeoretical and experimental -resuits gare 
compared for a square plate and a shallow spherical shell 


with an hexagonal boundary subject to blast loading. 
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A list of all important symbols in the text is 
given here. Symbols which are introduced in some sections 
Ol vie Lexy but which are not referred to later are not 
included. Some of the symbols may have different meanings 
in -differentwsectionasrnowever; these have, been.clearly 
defined in each section and should not be confusing to 


the reader, 


a radiussofecirclesncircumscribing regular 
polygonal boundary 


a radius of circle inscribing regular polygonal 
© boundary 

R radius of curvature of middle surface of shell 

h thickness 

p number of boundary sides 

E modulus of elasticity 

Vv Poisson's ratio 

Y mass per unit volume 

m mass per unit area, m= yh 

D flexural rigidity, D = Eh*/[12(1 - v*)] 

r radial polar coordinate 

p dimensionless radial polar coordinate, p = r/a 

6 angular polar coordinate 


5,0 polar coordinates of boundary points 
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number of collocation points on the 
characteristic segment boundary 


central angle of characteristic segment, 
ae = t/p radians 


Cartesian boundary coordinates normal and 
Peaneenvial tosune edge, respectively 


transverse coordinate normal to the middle 
surface; also used as an integer ranging from 
0 to © in the functional series solutions 
which describe each mode of response 


time 


transverse load intensity per unit area, 
dn= q,a4(p)cos pm@ Q(t) where m is any integer 


reference transverse load intensity per unit 
area; for blast loads q. is the peak reflected 
overpressure; for static loads q, is the 
static overpressure 


normalized transverse load distribution in 
the radial direction 


normalized transverse load transient 


concentrated central transverse load, 
P.=.P,Q(t) 
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eirculam Cnequency. of. .-—th mode..of vibration, 


w, = (k, /a) ?¥D/m for plates, 
w, = C1/y) [(D/h) (k, /a) * + E/R?] for shells 


Shell constant associated with i-th mode, 
A, = Eh/(maRw; ) 


damped circular frequency of i-th mode of 


d 
vibration, w Se ae 


coefiicient. of -viscous damping associated. with 
the lth mode, of. vibration, 64/Ws is. thenratio 
of damping to critical damping 


transverse displacement 
dimensionless transverse displacement, n = w/a 


transverse displacement eigenvector of the i-th 
mode of vibration associated with the homogeneous 
edge conditions 


stress function 
dimensionless stress function, & = F/Eha? 


stress function eigenvector of the i-th mode 
of vibration associated with the homogeneous 
edge conditions 


dimensionless transverse displacement and stress 
function, respectively, associated with the 
homogeneous edge conditions 


time functions associated with the &-th 
nonhomogeneous edge conditions in functions 

of transverse displacement and stress function, 
respectively 


dimensionless transverse displacement and stress 


function, respectively, associated with e® and 
fet eS pecul vey. 
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integration constants associated with the i-th 
mode of vibration and determined from the 
homogeneous edge conditions 


integration constants determined from the &-th 
nonhomogeneous auxiliary edge condition in 
transverse displacement 


integration constants determined from the -th 
nonhomogeneous auxiliary edge condition in 
stress function 


modal participation function of the i-th mode 
Of vibration_for uniform loads associated with 
the homogeneous edge conditions 


modal participation function of the i-th mode 
ot vibration fori uniform, loads associated with 
the homogeneous edge conditions and calculated 
using only the zero-order terms of At 


modal participation function of the i-th mode 
of.vibration for central point loads associated 
with the homogeneous edge conditions 


modal participation function of the i-th mode 
OfevVibravion Lorsuniform loadssassociated with 
the homogeneous edge conditions and calculated 
using nodal areas 


modal participation functions of the i-th mode 
of vibration associated with the homogeneous 
edge conditions and the £-th nonhomogeneous 
edge conditions 


in-plane displacement component in the ay 
direction, where a, = r,0,v,s 


stress resultant component per unit length 


Whichsis poSivive on the posivi ys say etement 
face in the a, direction, where a, ,a,=r,0,n,Vv,5s 
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stress couple component per unit length with a 
vectorial component (right hand rule) which is 
positive on the positive a, element face in the 


a, direction, where a,,a, = r,0,v,s 


strain.component.which is.positive on the 
positive.a .element.face,in thea, direction, 
WREPewd. Oo, c= rn, V,s and a, = m5b, (blank), 
indicating membrane, bending and membrane plus 
bending, strain, .respectively 


dimensionless static central transverse 
displacement for uniform loads, W= 2 W, = 
= wD/(q,8¢) 46 ba 
dimensionless static central transverse 
displacement for uniform loads associated 
with the i-th mode of vibration 


dimensionless static central transverse 
displacement for central point loads, W* = 
= 5 WP = wD/(P_a?) 

Th we Cae 


dimensionless static central transverse 
displacement for central point loads 
associated with the i-th mode of vibration 


dimensionless static central transverse 
displacement for uniform loads, W® = wD/Lq (pra) "J 


dimensionless static central transverse 
displacement. for central point, loads, 
WPse = wD/[P.(p&a)*] 


dimensionless bounding radius of a surface 
having the same area as that of the surface 
bounded by a@ particular regular polygon 


coefficient of edge clamping relating 
dw(p,8)/dv to Mie (5,0) for the i-th mode 
of vibration 


coefficient of edge clamping relating Fys(B,8) 
to E55 (P59) for the i-th mode of vibration 
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standardized n-th order Bessel function 
of the first kind of argument x 


standardized n-th order Bessel function 
of the second kind of argument x 


standardized n-th order modified Bessel 
function of the first kind of argument x 


standardized n-th order modified Bessel 
function of the second kind of argument x 


linear homogeneous differential operators 
of order 2 or lower defined in Eqs. (3,20a) 
throughecs:206) 
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CHAPTER 1 


INTRODUCTION 


fe UD eculves 


The purpose of this thesis is to present approx- 
imate theoretical solutions for the free and forted 
linear transverse vibrations of thin homogeneous isotropic 
elastic plates and spherical shells of constant thickness, 
both having regular polygonal boundaries. The theory 
gives the linear viscously damped response of these plates 
and shells to arbitrary transient symmetrically distrib- 
uted overpressures for a wide variety of homogeneous and 
time-dependent edge conditions. 

Theoretical frequencies, mode shapes and 

modal participation functions for uniform loads and cen- 
tral point ioads have been tabulated: for many different 
simply supported and clamped polygonal plates. Since 
the numerical application of this theory requires ex- 
tensive use of a digital computer, it was impractical to 
tabulate results for shallow shells having all of the geo- 
metrical properties and edge conditions which might be of 
LHVErest. However, a sufficient number of different 


shells have been analysed to establish some general 


£ ASTIARS 













UOITOUGOATUI 


-xo1ggs Snmece tq oF el eteedty etd? Io ssoqiugq sAiT 
besi012 bas satt sdt tot smoltuloe I[softesosd? ovamt : 
otgottoet avyosnsegzomod nidt to arottsidiv setevenstd soon 
,eaendolds tnustanoo to ellede Lasiresdiqs bas esisiq ottveste 
vtosds sdiT seicnierenniaee [snogylog tsluget ynived tng 
setsiq seers to samnielivatie beqmsh ylevooetv asenki edd aevig 
-dittetb yilsotatemmye tustenatsd yisttidis ot efieda bas 
Das suyoensgzomod to ytettsvy sbiw s 10l esivesstgIsvo sage : 
.enioltibnos sybe inebnegebuamid | . 
bas esqade sbom ,astonsupett L[sottetosctT t 
ees bas ebsol miotiny tot enoltonuvt noidagqtotiasq £ 
Ynasisttib yasm 101 Ssvsludst need svaed absol talog Ist — 
sont2 .aetslq Ismogylog bequisio bas betroqqua yl 
-%9 sextupe1 yroon? etdt Io mottsoliqas [sottemwa ¢ it 
ov isottoaiqmt sew Jt ,tesuqmoo [stigth eg Io sav 
joa ro {is achved affeds wolisde 10% naluewt:# 
| neers sabe brs set 
ys y 











nite sna sence sip 





_ 
Pas ty a 


se dell 


effects of variations in the shell geometrical proper- 
ties and edge conditions on the natural frequencies and 
modal’ participation functions. Computer programs which 
calculate free vibration and forced response of these 
peaves and shelis are included. These programs are 
written in Fortran IV language for use on IBM 360/67 
and 1130 (with modifications) digital computers. 

tt is a major purpose of this*thesis to pre-= 
sent detailed comparisons between transient experi- 
mental and theoretical stress resultants (forces) and 
stress couples (moments), obtained from the response 
to air blast overpressures of a square plate and a 
shallow shell with an hexagonal boundary. These com-= 
parisons reveal some of the limitations of the theory 
and some of the difficulties which may be encountered 
when an attempt is made to correlate simultaneously 
both magnitudes and frequencies of the sectional re- 
SUuLcaive fOr precuical applications. The erfecu of 
variations in shell geometry from the ideal, espec= 
deliy in ragius- o. curvature, is considered = oo play van 
important part in the forced response of thin shallow 


shells. 
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wee SCUSBS LOnPOor sb teratvure: for Plates 

Literature representing some of the different 
analytical approaches to the solution of the linear 
Peenoverce st exural vibrations of elastic plates is 
Summarized in this section, From,this survey,.it is 
evident that while much work has been done on free 
wibrations of plates,,theoretical forced.response has 
been computed for only a few plate shapes and detailed 
experimental results for the forced response charac- 
teristics of a complete plate are lacking. 

ioe e per uienial Results. Lassiter and 
Hess [1]! studied the response of a clamped rectangu- 
Pam be ev OnraucOm oCcousticC. loading measuring strain 
P@eovoeOOsc. VlOb Gully walet ne middie, of. the shorteside., 
They concluded that for flat panels there 1s 4 non= 
Pivesr ay elaivOlvineed esti fening spring .constanv and 
Porscuryea panels the nonlinearity involves a jde— 
Creaciieesprinie. CONS lant, they, found, also, suhas 
combined structural and air damping was independent 
of panel thickness but increased rapidly with stress 


atc higher stress levels, 


1 Numbers in square brackets [ ] refer to references 


in Bioliography. 
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Crocker [2,3] has studied the response of rect- 
angular panels to normal and travelling sonic booms and 
to step and oscillating shock waves. He showed that 
contributions of higher modes to bending strains (in the 
direction of the longest sides of length a) are sig- 
nificant, particularly as the panel aspect ratio a/b 
(where b is the length of the shortest sides) is increased. 
He compared the experimental and theoretical center 
transient, bending strain of a clamped 4.5°x 2.5 in. panel 
subjected to a N wave and found fair agreement. 
Arising from the present work, Walkinshaw and 
Kennedy [4] have given the approximate theoretical sol- 
Ution tor the forced response of polygonal plates subject 
to homogeneous edge conditions and compared experimental 
and theoretical transient central displacements and 
moments at different positions for an elastically clamped 
square plate subjected to air blast loads. 
woe eukecu .neorevical SoluDlons a haacy 
solutions for the free linear vibrations of circular 
plates have been obtained for various edge conditions. 
Radially symmetric free vibrations were analysed first 
Dye Pol esor 5 |. Later Kirchoff [6] calculated non- 


axisymmetric frequencies for a circular plate with free 
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edges. Airey [7] and Rayleigh [8] analysed free vi- 
brations for various edge conditions. Mindlin [9] 
enowear the erfect of neglecting rotatory inertia and 
transverse shear, termed the thickness-shear mode effect, 
Onan zheremodes of vibration. Coupling of the thickness- 
shear mode with the higher flexural modes is shown by 
Deresiewicz and Mindlin [10] for the axisymmetric vi- 
brations of a circular plate with free edges and by 
Deresiewicz [11] for the axisymmetric vibrations of a 
circular plate with clamped edges. Raju [12] calcu- 
Jated the natural frequencies of annular plates with a 
variety of edge conditions at the inner and outer bound- 
aries. Reid [13] analysed free vibration motion of a 
Circular platej subjected to an initial displacement 

and velocity distribution and Harris [14] solved the 

free vibrations of a circular plate having a lenticular 
PaLeKticss Varilav von. 

Forced movion olf circular plates wundergoing 
symmetrical deformation has been analysed by Sneddon 
[15] for an artificial boundary condition using in- 
tegral-transforms. Response of clamped circular 


plates has been analysed by Flynn [16] for impulsive 












=? é : ae 
a wv age ae 4 apy f 3iT & , 3] siglel ve ae Sis ut 
- . 


y gobi tate snotgibsioo e_7os suo lisv haa anol. Br 


7 Oly Gee 4: t br 74 BOR " Q Z 597 Is ig ‘bever iB. 


ie 


ett sbom saesrofothd sfv bolwset . sede osteVene 


ap te 
"Pr 7 t. 7 ‘Pero! totustdiv ta sebom  wSeniein ae 
d ‘ . od > a - 
- «ie 
wode el esbom Istwxs/ 1 rweAigin ot déiw sho 6she 
| i ae 
map e wtve tT 3 aig t6 1 Os iLiba lM Dos sotwelee ted - 
wd - 
. +’ a 
ee saet dtiw etal rafvorlo 8 lo enolvsia = 







7 efit gifiey mofsfOnoo yrsbruvod istolti¢«se as. t02 











a Be _ r %& 7 Sof pears 
“sluotio begmalo te satogest eon at 





tol [OL] tereltt vd beayisns | rasd aad & 





\ 


loads and by Reismann [17] for concentrated harmonic-= 
eliveosc! lating loads. McLeod and Bishop [18] have 
calculated the forced symmetrical and nonsymmetrical 
Vibmatlonssorlclpcularm plates subjected to a central 
force Or a line force at any radius .and of annular 
plates subjected to moment excitation at free and 
Pinneasecarvest@andstorce excitation at free and sliding 
edges. Weiner [19] has analysed the forced motion 
of an elastically built-in circular plate under con- 
Cent rvce rane Voading, central loading and uniformly 
distributed loading. Reismann and Greene [20] have 
analysed the forced motions of circular and annular 
plates subject to general stationary or time-depend= 
ent boundary conditions using the "improved theory" 
which includes the effect of transverse shear de- 
formation and rotatory inertia. Their comparison 

of central displacements and central and edge moments 
calculated by the classical and improved methods for 
a clamped circular plate reveals the increasing im- 
portance of the shear deformation with increasing 
ChLekness=-to-radius ratio h/a and decreasing ratio 

of central area uniformly loaded to plate area. They 


concluded that the improved theory should be used for 
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uniformly loaded plates having a thickness ratio h/a >0O.1 
and showed that the central displacements predicted by 
thes@lagsicalsGhneory will be too small both for static 
and dynamic loads. 

HeacClmeOluvlolsiorethne linear yiprations of 
rectangular plates with simply supported edges have been 
given by Rayleigh [8] and Timoshenko [21]. Voight [22] 
Obtained an exact solution for the free vibration of a 
rectangular plate simply supported along one pair of 
parallel edges and free at the other two edges. Mindiin, 
Shacknow and Deresiewicz [23] investigated the influence 
of rotatory inertia and shear deformation on the flexural 
vibrations of rectangular plates with all edges simply 
Supported and with two parallel edges simply supported 
andetne.OUner two. [ree. 

1.2.3 Approximate Theoretical Solutions. 
Approximate solutions have been obtained for the linear 
transverse flexural vibrations of elastic plates having 
complex geometrical and material properties subject to 
a variety of boundary conditions. A more complete 
summary of approximate procedures than is presented by 
this review is given by Desai [24] and expanded by 


Meirovitch [25]. 
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Ritz [26] introduced the energy approach, in 
Wiich the calculus of variations is used to find the 
minimum of potential energy for self-adjoint systems, 
for free vibrations of rectangular plates with edge 
eenditions other than having all sides simply supported. 
Natural frequencies of a square plate with free edges 
were found by Ritz [27] using the energy technique. 
Using the energy approach, Iguchi extended his earlier 
work [28] to find natural frequencies of rectangular 
plates with various length-to-width ratios and edge 
conditions [29)]. Young [30] used the Ritz method to 
calculate natural frequencies of rectangular plates 
for various edge conditions and Ota and Hamada [31] used 
this approach to obtain fundamental frequencies of sim- 
ply Supported rectangular plates clamped along portions 
of the edges. 

Nowacki [32] analysed the free vibrations of 
a clamped rectangular plate by obtaining an integral 
equation from compatibility of a series of clamped 
strips. Nowacki [33] has given a solution for the 
linear forced flexural response of a rectangular plate 
with simply supported edges subjected to uniform or 
concentrated loads using Greene's function and the 


orthogonality condition. 
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Tolke [34] introduced the boundary colloca- 
Cron vecnnigques, woien can pe*used To satisfy arbitrary 
edge conditions Tor*piates ofvany shape, for the static 
solution of square plates with clamped, simply suppor- 
ted and free edges. This technique has been employed 
by Conway [35], who introduced the nomenclature "point- 
matching", and by Leissa [36] to obtain the fundamental 
frequencies of some regular polygonal plates. Conway 
and Farnham [37] used this technique to obtain funda- 
mental frequencies of some polygonal plates and Chen 
and Pickett [38] have developed some computer programs 
vo ODvLaIn Navural™ irequenctes ‘for plates” of any shape 
subject to any boundary conditions. 

Free vibrations have been analysed using the 
Ritz energy method for rectangular and skew cantilever 
plates by Barton [39], for triangular cantilever plates 
by Anderson [40] and for some irregularly shaped plates 
with square and circular holes by Rao and Pickett [41]. 

Reipert [42] extended the Levy method to 
analyse vibrations of polygonal plates resting on el- 
astic foundations subject to transverse and in-plane 


loads and various boundary conditions. 
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Solecki [43] employed Fourier transforms to ob- 
tain the free vibrations of freely-supported four sided 
plates with two opposite parallel rectilinear sides and 
two other sides of arbitrary shape. 

Shahady, Passarelli and Laura [44] used com- 
plex variable theory to map conformally various plate 
shapes onto a unit circle to obtain the fundamental 
frequencies of vibration of these plates. 

Pandalai and Patel [45] analysed natural fre- 
quericies*of*vibration of ‘circular’ plates having diff- 
erent material properties in the radial and circun- 
ferential directions. Salzman and Patel [46] analysed 
the asymmetric natural vibrations of orthotropic cir- 
cular plates with variable thickness using the method 
of -Frobenius'. 

Dawe [47] applied the discrete element dis- 
placement method to predict changes in natural fre- 
quencies of lateral vibration due to axial and biaxial 
in-plane loading for some simply supported and clamped 
rectangular plates. 


Les -Discuss On onebacverature forushallow Sphericale cheiis. 


Literature representing some of the analytical 
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solutions for the vibrations of shallow spherical shells 
is summarized in this section. A more extensive review 
of the analytical methods of analysis for free vibration 
and transient linear elastic response of general thin 
shells is given by Kalnins [48]. It is evident that 
although much work has been done on computing natural 
frequencies of shallow spherical shells with circular 
boundaries, theoretical forced response of these shells 
has been computed for only a few different edge con- 
ditions and thinness and shallowness parameters. 

Since most available experimental verification of shell 
frequencies and mode shapes have been done during 
stability investigations, some literature on the buck- 
ling of shallow spherical shells is included. Dis- 
crepancies between the theoretical and experimental 
results given for the buckling of shallow shells in- 
dicarewthati itheweffects of “imperfections: in shell 
geometry and the inability to satisfy idealized edge 
conditions exactly Houta be studied experimentally 
horebiiie ar qv.brataon si: No detailed experimental and 
theoretical comparisons for the linear forced response 


Ofeshalstowns heliswcould! be: found’ inthe literature, 
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Deol LimearsVibrations: An exact solution 
for the asymmetric free vibrations of shallow shells 
with clamped edges was formulated by Reissner [49] who 
evaluated an approximate fundamental frequency using 
the Rayleigh=Ritz procedure. Later, Reissner [50] 
showed that longitudinal or in-plane inertia could be 
neglected in solving for the transverse free vibra- 
tions of shallow shells, thereby reducing the problem 
to the solution of two simultaneous differential 
equations involving a stress function and the trans- 
verse displacement rather than three simultaneous 
differential equations in the three displacement com- 
ponents. Using this approximation, Reissner calcu- 
lated the frequencies of simply supported shallow 
shells of rectangular plan with sliding edges. 
Neglecting longitudinal inertia, Reissner [51] de- 
termined axisymmetrical natural frequencies of 
shallow shells“with*both*clamped*and free edges, and 
forced vibrations of an tnifiinivewshaldvow *shelrewi on 
a=poanty Lload*av'thetapex’ Free transverse vibra-= 
GiGho Or shallow spherical, shells with -Tfreesedces were 


analysed by Johnson and Reissner [52]. 
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Naghdi [53] derived the differential equations 
for the linear response of thin shallow elastic shells 
including the effect of transverse shear deformation, 
Kalnins and Naghdi [54] gave an analysis for the 
coupled longitudinal (torsionless) and transverse 
asymmetric vibrations of shallow spherical shells. The 
ditfferentialeequations “governing vibrations of shallow 
spherical ishells elineluding jthe coupling effects of 
bongitudinalwatransversepand rotatory inertia, and 
transverse shear deformation were given by Kalnins 
[55], who calculated the frequency equation for a 
clamped shallow spherical cap from these differential 
equations. Free nonsymmetric and axisymmetric 
vibrations of clamped shallow spherical shells, 
including some mode shapes, were computed by Kalnins 
[56]. He showed that neglecting longitudinal inertia 
Pesusvedeinathe omissionpof ran infinite numberof 
in-plane inertia modes but that transverse vibration 
frequencies computed ae the classical method were quite 
accurate and essentially uncoupled from the in-plane 
frequencies, 


Van Fo Fy [57] introduced two auxiliary 
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functions’ for shallow spherical ‘shells with longitudinal 
inertia neglected, thereby reducing the two simultaneous 
differential equations in stress function and transverse 
displacement to two simpler uncoupled differential equa- 
tions, the solution of which reduces to the solution ob- 
tained by Reissner [50]. 

Lee [58] analysed the axisymmetric vibrations 
of a shallow spherical shell with a heavy mass at the 
apex and Pandalai and Dym [59] showed that for practical 
purposes the transverse shear force can be neglected in 
the calculation of the axisymmetric natural frequencies 
of a shallow spherical cap. 

Koplik and Yu [60] studied the axisymmetric 
vibrations of homogeneous and sandwich shallow spheri-= 
cal caps with clamped edges including the thickness- 
Shear effect. They showed that the frequencies of the 
higher modes of vibration of homogeneous and sandwich 
Capercen be calculated from the uncoupled equations of 
clamped circular plates as the shell curvature affects 
only the lower frequencies, and that the thickness- 
shear effect can be ignored for homogeneous caps. 

Koplik and Yu gave an approximate solution for the 


axisymmetric vibrations of clamped spherical caps [61] 
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and analysed the torsional vibrations of homogeneous 
and sandwich spherical caps and plates with clamped 
and free edges [62]. 

Lock, Whittier and Malcom [63] showed that 
the nonsymmetric transverse vibrations of a clamped 
Shallow spherical shell depend upon a single shell 
geometric parameter, and that the eigenvalues of the 
higher modes converge to the clamped plate eigen- 
values. 

The response of shallow viscoelastic 
spherical shells to arbitrary time-dependent axi- 
Symmetric loads has been studied by Naghdi and 
Orthwein [64]. Neglecting longitudinal inertia, 
they solved the differential equations for unlimited 
shallow shells and shallow shell segments using 
Laplace and Hankel transforms, Van Fo Fy and Buibol 
[65] obtained a solution for shallow spherical shells, 
neglecting longitudinal inertia, subject to a harmonic 
axisymmetric ete ee load. Kraus and Kalnins [66] 
derived a general solution for arbitrary elastic 
shells subjected to time-dependent surface loads, for 
which the free vibration characteristics are known, 


using spectral representation or normal modes, They 
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16 
also showed that for a shallow spherical shell with 
simply supported edges having a base circle-to-thickness 
ratio of cO and a half angle of opening of 15 degrees, 
three symmetric modes usually are required for the 
transient response to a suddenly applied uniform 
pressure while for a hemispherical shell with the same 
thinness parameter, fifteen modes are required. The 
method of analysis for obtaining the response of 
shallow spherical shells” subjected to transient surface, 
edge and thermal loads using spectral representation 
has been outlined by Kraus [67]. Reismann and 
Culkowski [68] have obtained a solution for a clamped 
shallow spherical shell subjected to an impulse uni- 
formly distributed over a central area including the 
effect of shear deformation and rotatory inertia. 

1.3.2 "buckling: Hossack [69], in the 
course of an experimental investigation on the non- 
linear behaviour of shallow spherical shells subjected 
to static overpressures, Suggested that shells manufac- 
tured by spinning are subject to residual stresses 
and variations in thickness. He stated that the 


clamped edge condition is almost impossible to achieve 
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in practice and that attempts to attain it will induce 
random forces and moments and, in any case, in most 
engineering applications, the. behaviour. of a, shallow 
cap would approximate that of a freely supported shell 
even, when. some measure of restraint is provided. 
Hossack found for some freely supported caps that ex- 
perimental bending strains agreed more closely with 
theory than did the membrane strains. 

ins experiments. on, dynamic. buckling of 
clamped shallow spherical shells under uniform press- 
ure pulse loading, Humphreys, Roth and Zatlers [70] 
found that maximum prebuckling deflections exceeded 
the theoretical predictions of Budiansky and Roth 
[71] by as much as 100 per cent. 

Huang [72] suggested that when the shell 
height-to-thickness ratio exceeds a certain value, 
initial geometric imperfections or disturbances 
during deformation would cause the axisymmetric de- 
formation of clamped shallow spherical shells to bi- 
furcate to asymmetric deformation at some prebuckling 
pressure. Using this assumption, Huang, predicted 
lower buckling pressures for these shells which, 


however, were still not conservative in comparison to 
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experimental values. 

Archer and Famili [73] adopted the vibration 
method of stability analysis to theoretically predict 
static buckling pressures for clamped shallow spherical 
shells. This method is based on the fact that the 
frequency of the buckling mode, which can be. predicted 
by Huang's theory [72], diminishes from the free vibra- 
tion values given by Kalnins [56] to zero frequency 
when the shell is loaded at the buckling overpressure. 
Okubo and Whittier [74] applied this nondestructive 
method of determining the static buckling pressure to 
some clamped shallow spherical shells. They found 
that the experimental frequency of the axisymmetric 
buckiing mode, predicted as the failure mode by Huang, 
was slightly lower than the value predicted by the 
metnod of Archer and Famili at all static overpressures 
for all shells tested. They suggested that radius of 
curvature variations in the experimental shells as 
well as thickness Ben iactons might explain the dis- 


crepancies. 
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CHAPTER 2 


PLATE THEORY 


2.1 Governing Equations. 


The fundamental differential equation of motion 
given by Poisson [5] and Kirchoff [6] for thin elastic 
plates subjected to linear isothermal deformation by a 
normal dynamic massless loading function q(r,®,t) with 
viscous damping included, neglecting the effects of 


rotatory inertia and transverse shear, is 


ow(r,6,t) aewlr,6,t) 
D V*w(r,0,t) + 2¢m —————_ + m————- = q(r,@,t). (2.1) 
at ot? 


In this equation w is the transverse displacement of 
piresmiug2eeplaneé or the platey c is the coefficient of 
viscous damping, m is the mass per unit area, 
Peeeone eles) 1s the flexural rigidity, ’h is the 
plate thickness, E is the modulus of elasticity, v is 
Foveson s ratio; r and 9 are the radial and angular 
9? linye Eye 
polarescoordinaves, tuis time,) and V7 = —> +7) 
antng neeorme oO” 
when operating on w(r,6,t). 
Introducing the dimensionless parameters 
W r 


Nn =— and po =-—, where a is the radius of the circle 
a a 
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circumscribing a regular polygonal plate, Eq. (2.1) 


becomes 
' dn(p,0,t) d*n(p,8,t) 
D V'n(p,6,t) + 2tma* —————— + ma* ——————- = 
at at2 
=52. Gi.pacnis)s G2 #2) 


: 9? 19 dip gd 
where V“ = —— + — — + — —— when operating on 
Cpe ope 0-90 - 
nip,0,t). 

In the interest of obtaining the most 
general solution, the possibility of time-dependent 
nonhomogeneous edge conditions is included. A 
solution of Eq. (2.2) for nonhomogeneous boundary 


conditions which depend on functions of time a*(t) is 


assumed in the form [25] 


(e 
n(p,0,t) =n (0,8,t) + alanis vn (2.3) 


L=1 
where c = 4 for a plate with interior and exterior 
boundaries and c = 2 for a plate possessing only an 
exterior boundary. The functions 5*(0,8) are chosen 
to make the boundary conditions expressed in terms of 
n (p,8,t) homogeneous. Substituting Eq. (2.3) into 


Eq. (2.2)yields 
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_ Odmt (E.5) .o8 yotivditedy® .asoensyomod (tee don 


abletyi $.8)% 
7 - 


dn (p,6,t) , 2-7, (p,8,t) 


Pevens(o,0,t) +, 2¢ma" + ma 
: ot at? 
Q 
Cee 2 Q Curent) pace 
= a*q(p,0,t)-D } & (t) V*E (p,6)-2¢ma* } ——— & (9,6) - 
Q=1 fel dt 
Q 
ey cul toe 8 
- ma* ——— ¢£€ (9,6). (2.4) 
gal) dt 


2.2 Free Undamped Motion 
2.2.1 General Solution. The homogeneous 


equation of the undamped system 


mage any 0,05 ¢) 


Ven (p,6,t) 1 tees 
i D ote 


= 0 (275) 


can be solved by the separation of variables technique 


n,(es8,t) = n(p,0),exp(iut), 


where w is the undamped circular frequency. Sub- 
Beatie von Into. bq. (2.5) gives two simpler 
differential equations of the form 

CV 7 4ak*) n, (58) = 0 (OAS, 
and 

(v2 - k?) 7 (p,6) = 0, (oa) 


where 
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For linear problems 
n(p,®) = n (058) + n, (9,8). 
Love [75] showed that solutions for equations 
of this type can be given as 


Co 
n(p,9) = V (p) ,cos nO . 
& 2 n {sin n@ 


Consequently the general solution of Eq. (2.3) is 


given by 


@ 
nose )eemye PA Ia(ko)+ Bol (kp) + -Caywvdiko) + 


n=0 ni on noe) ai tl 
cos né 
+ ue Sg 2 aah sors (2.8) 


where A,B, C and D_ are integration constants, 
ae n°) on n 


Lis and a are Bessel functions of the first and second 
Kit respectively p-and I, and K, are modified Bessel 
fUlevrTOons=otsine first and second kind, respecuively. 
Restricting the solution to plates possessing no 
interior boundaries yields 

C =D =0 

n n 
and applying the solution to plates of regular polygonal 
shape possessing p-ply rotational geometric periodicity 


1 
as shown in Fig. 2.1 yields for symmetric response 


LL 


‘This assumption of symmetric response restricts the 
loading function to be of the form q(p,0,t) = q,a(p)cos pmo Q(t), 


where oe is a reference pressure and ™ is any integer. 


(8.9) p fl + ( ‘<< ; 
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24 
n(p,8) = ie emeko)et BT U(ko) ] "cos pnd, (259) 
n=0 pn pn Pn pn 
6 is 0 and t/p radians, respectively, along radial 
lines which pass through the boundary center and 
eorners. Consequently, the response is symmetric 
about these radial lines if n(6) is chosen to vary as 


cos pné, since 0n(p,0)/06 = dn(p,t/p)/a6 = 0. 


2.2.2 Edge Conditions, The homogeneous edge 
conditions which must be satisfied for an elastically 


built-in plate are 


n(6,6) = 0 (2,10a) 
and 4 
=~. dn( 6,98) 
M0 0s) ei 8 (2,10b) 
Vs 1 dv 


where 0 and 6 are the coordinates of boundary points. 
Tnebow.(2.10b) M (P29) is the flexural stress couple 
vector tangential to the boundary and 8, is the edge 
spring constant of proportionality ranging from zero 
for a simply supported plate to infinity for a clamped 
plate. The boundary and interior coordinate systems 
are shown in Fig. 2.2. The associated nonhomogeneous 


time-dependent edge conditions will be dealt with in 
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FIG. 2.2 Vector Diagram Showing Relation of Boundary Coordinates v, s to 
Interior Coordinates r, 0 
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Chapter 3. 

Since at the present time an exact solution 
for this problem is not available for the edge 
restraints in Eq. (2.10a) and (2.10b), Tolke's 
boundary collocation method [34] is employed, whereby a 
rigorous satisfaction of the prescribed boundary 
conditions is collocated for a number of discrete 
points located on the boundary of: one of the rotationally 
periodic segments of the plate. 


Expressing Eq. (2.10b) in polar coordinates 


yields 


M (6,6) cos?6 —-M (6,8) sin?e +™M (p,6) sin 26 = 
rg 6r re 


—_ 
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ah 1 anke,e9 
: eye) 0 96 


clekcy Gh ey | (20c ) 


The flexural stress couple components Megs 


and M_ , the torsional stress couple components Woe 
Or 


and M6 and the transverse shear stress resultant 


components F and as , which are: shown dnwhieiec.>. 
rn n 


can be expressed by 
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FIG. 2.3 Vector Diagram of Plate Element Showing Sectional Resultants 
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Substitituon of Eqs. (2,9) and (2.11) into Fqs, (2,.10a) 
and (2.10c) produces linear homogeneous equations, the 
Solu onsOle Which Williavield an infinite set of 


eigenvalues k and their associated eigenvectors. Thus 


CO 
yee feet be yal = 0, (2,10a*) 


CO 5 a 
at 7. a. 
yA We 8 —v) +B (wy - 8B —¥,)] = 0. (2,10c*) 
5 1 3 6 
D D 
The coefficients in the boundary equations (2.10a*) 
and (2.10c*) are 
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v, aoe = COR 8 cos pné Fee (fi) sin @ sin pné, 
1 pn 0 pn 
v =k I (fi) cos 6 cos pnd +22 1 (fi) sin 6 sin pné, 
4 nh Veye) 5 pn 
ee er ee inte < Ks ~ 25 A 
sys com w Cl) = ioe eG) ot Vi). 0? (ft) cos-8 cos pnd. + 
ah Mele! 0 pn p pn 
n Ks x a 
FLV TPC) = (EM 7° (fi) + (—)o' (7) sin?6 cos pnd - 
1 pn 6 pn peep 
a (Ch es OS Ty er ere (Pe) a! (fi)] sin 26 sin pné 
~2 q\% 3 
(she SAE Pp pn 
n Ks - ~ 
es [kj edd - v(Be)? I (fi) + v(—) I' (ji)] cos*6 cos pné + 


6 pn fe) pn 


+ 


k A ~ 
[vki T2.CH) - (2B)? 1 (i) + (4) 1t (a) ] sin?6 cos pné - 
£ p 


(1 - vals Ipn(fi) - ky (BR) Tr (fi)] sin 26 sin pne, 
(@) 


where ti = k 6. 


2.3 Forced Motion 

2.3.1 General Solution, The classical method of 
spectral representation in which the dependent variables are 
expanded in an infinite series of normal modes of free vi- 


bration is assumed for n,(p,8,t), defined by Eq. (2.4), in the 
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form 
OO ° 
Omer eaee en cp, ) eT A(t), (2,12) 
: i=l 
where 
i ee Al ct 
Recs Ope) LAL eRd (kee) + BOI ({kap) Jeog pné feeise 
n=0 pn pn pn pn 


and tit) is usually referred to as the modal participation 
Ge cor me oubst i tUting Ed. (2,le) into Hq. (2.4) and 


Simplifying yields 


co @@*pt(¢) att(t) 
>» [—— + 25, —— + w(t) ]At(,6) = 
ia) edt. + at = 
ie D 2 
=—- G(p) cos pmo Q(t) -— y 8*(t) v'z*(p,0) - 
ma ma* £=1 
2) de-(f) a, gece (0) eee, 
a et; ). Saas (p,8) = Ds = eto g (0, u26 (2,14) 
g=1 © dt g=1 = at 


Here it is assumed, as is done by Meirovitch [25], that 
the eigenvalues ks and the eigenvectors a associated with 
the homogeneous equation of the undamped system, Eq. (2.5), 
satisfy the nonhomogeneous equation with viscous damping 
Dmohuded ss Eqe. 2.4) 

Multiplying both sides of Eq. (2.14) by AI (p,8), 
integrating over the plate characteristic segment area and 
applying the condition of orthogonality of the eigenvectors 


pointed out by Rayleigh [8] and shown in Appendix B.1, 
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Eq. (2.15) has been simplified by arbitrarily choosing 
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ug 
& (p,®) as solutions of the biharmonic equation in the form 
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This is permissible since these auxiliary edge condition 


functions are not unique, as Meirovitch [25] points out, 


2.3.2 Initial Conditions. The solution of 
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t 
Lyi igen a2a*(7) 
@,7 (t) = (w2/a®) i ——-— expl-t,(t=0)] x 


0 
a 
XeS 1 Deut TdT (2.21c) 
i 
d 2 2 
and a = ju? = oF is the damped circular frequency of 


the i-th mode, 
If the plate is initially assumed to be at 


rest and to have undergone no transverse displacement, 


then 
n(p,6,0) = 0, (2a 22a) 
dn(p,8,0) 
ot 


and the free part of the solution, Eq. (2.20), vanishes 


leaving 
T(t) = C1/(w? yxt)] C(aj/ma) xt at (t) - 
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participation functions 


a al : ul j ah Gul 
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In evaluating xP and Xe terms involving 
products of functions of different orders n have been 
neglecved, inis Szmplification 1s exact “only for 
circular boundaries and introduces errors which will 
increase in magnitude with the deviation of the nodal 
lines from circles concentric about the plate center. 
The accuracy.of this simplification is* checked for’ some 


numerical examples in Chapter 4, 
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Evaluation of the coefficients $...6 | 
involves integration over op and 6. The integration 
over p can be carried out exactly as shown by McLachlan 
[76], for example, however the integration over 6 must 
De carried out numerically. If the terms involving 
products of functions of different orders n had not 
been neglected in Eq. (2.22), then formulas for their 
associated %-coefficients suitable for polygonal plates 
with any number of sides could have been solved only by 
numerically integrating over both coordinates p and 6. 


For a load distributed uniformly over the 


plate surface G(p) cospm@ = 1 and the coefficients 
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plate center, PS replaces q. in Ea. (2.19) and the 


ceoerficientiss fon Xs become 


The values for 9 and . contain the divisor 2pa? 


since x" is evaluated for the plate characteristic 
segment area with a circumscribing radius of unity 
rather than the complete plate area with a 


circumscribing radius a. 
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CHAPTER 3 
SHELL THEORY 


3.1 Governing Equations 

The differential equations governing the linear 
transverse vibrations of shallow spherical shells sub- 
jected to linear isothermal viscously damped deformation 
by a normal dynamic massless load q(r,6,t), neglecting 
the effects of transverse shear and rotatory and longi- 


tudinal inertia, are 


Eh 
V'F(r,6,t) - =I Vow n,0.c)) = 0, (3% 1a) 
' ae, dw(r,6,t) 
DAVAWGL Ont mte—FV LCN .O9o) + 2on ———_— 
R at 
d?w(r,6,t) 
+m ——— = q(r,6,t). (3.1b) 
at? 


In Eqs. (3.la) and (3.1b) w is the transverse displacement, 
Pereeune stress function, © is the coefficient of viscous 
damping, m is the mass per unit aréa,.D =,Eh® / [12(1-v7)] 
DomulewMLexUraLerieiolcy |v viS the modulus of Stasti cacy. 
hois thegsheblathicknéss § viis: Poisson, siiravto kh is: the 
radius of curvature of the shell middle surface, r and 6 
are radial and angular polar coordinates, t is time and 
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Introducing the dimensionless parameters 


_ W AGA 2ir 
Dates: ae ore and er where a is the radius of the 


circle circumscribing the shell's regular polygonal 


base, Eqs. (3.la) and (3.1b) become 


V°E(He6, tyes = V?n(p,0,t) = 0, (3.2a) 
u Epa \ on ben at? 
D V'n(p,0,t) + —a°V2E(p,6,t) + 2a*tm ————— + 
R at 
a°n(p,9,t) 
+ ma* S telgis in” coe Orety C3820) 
ot 


9? 19 Wenge 
where V7 = — + —- — + — — when operating on 
dp? 0 909 On 007 
bunceuLons. of op and 6. 
A solution of Fqs.,(3.2a) and «3e2b). for*™non- 


homogeneous boundary conditions which depend on 


functions of time aX(t) and P*(t) is assumed in the form 
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where c=4 for a shell with interior and exterior 
boundaries and c=2 for a shell possessing only an 
exterior boundary. The functions aCpnen and h* (0,6) 
are chosen to make the boundary conditions expressed 
in terms of n,(p,6,t) and E 6928,t) homogeneous. 
BUDSCICUCcIN@etds, 3.54) and (3,36) Into mas. 


(3.2a) and (3.2b) yields 
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3.2 Free Undamped Motion 


3.2.1 General Solution, The homogeneous 


equations of the undamped system 


TeeeC pa ,t) = = V2n,(p,9,t) = 0, (3.5a) 
Eh a nein se gc) 
D V"n,(p,0,t) + — a® V7E,(p,0,t) + ma* ——————— = 0, 
R at? 
3. Sbp 


can be solved by the separation of variables technique 
n (9,85) = n(p,0) exp(iuwt), (3.6a) 


E(p,0) exp(iwt), (3.6b) 


ba Vo s0, ©) 


and by introducing the auxiliary functions used by 


Van Fo Fy [57] 


nip,o) = V°bUCo,98) + AVC on6) Jor Cau 7a) 


E(p,0) V(p,0) += UCp,8) (3.7b) 


where \=Eh/maRw? and w is the undamped circular 
frequency. Substitution into Eqs. (3.5a) and (3.5b) 


gives 
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V*C(V"= k*)U(p,6) + AV"V(p,6)] = 0, (3.8b) 
a Eh 
where k* = —[mw? - at 


If k=0, then w= /Eh/mR?, \=R/a and Eq. (3.8a) 
is identically satisfied for all. V(0,6) while Eq. (3.8b) 
reduces to 
R 
V°[U(p,8) + a ae = 0, k=0. (359) 
The solution of Eq. (3.9) for a shell having no interior 
boundary and a regular polygonal exterior boundary 


possessing p-ply rotational geometric periodicity is 
R 
U(p,.8) + — V(p0,8) = 
cos pné 
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= ss cos pné 
= n 
ECO. 0) a=a.)) E aie + D2 oP { » k=0, (3,11b) 
n=0 L pn BY sin pné 
Sa 2 a 
where C =m—C* and D*? #-=p' , 
pn R pn pn R pn 


If k#0, then Eqs. (3.8a) and (3.8b) reduce to 


two uncoupled differential equations in U and V 


V"v(p,0) = 0, (Seize 


Vilv" —.K"4U( 9,8) =,0% (20a 2b) 


THEecOlULLON Of Eqs.) (3.lea) and (3,l2b) for a shell 
having no interior boundary and a regular polygonal exterior 
boundary possessing p=eply rotational geometric 


periodicity is 
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Eas. (3.7a) and (3.7b) and restricting the solution to 


symmetric response! yields 
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These solutions for n(p,6) and &(p,@) are 
valid for all k since they reduce to the solutions 


given in Eqs. (3.lla) and (3.1lb) for k = 0 where 
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1 The assumption of symmetric response restricts the 
loading function to be of the form q = fhe a(p) cospmé Q(t), 


where q is a reference pressure and m is any integer. 
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3.2.2 Boundary Conditions. The theoretical solution, 
which is based on supplementing the force and moment equili- 
brium equations by the force and moment compatibility rela- 
tions, yields the transverse displacement and stress function, 
In-plane displacements uy, and u. must be determined by integra- 
tion of the strain-displacement relations, This integration 
yields expressions for the in-plane displacements in terms of 
the stress function and integration functions which could not 
be determined,! Consequently, edge conditions which could be 
expressed in terms of the in-plane displacements if the alter- 
nate theoretical approach [54] had been used, are expressed 
in terms of other known functions, 

Time-dependent edge conditions which can be satisfied 
approximately by the collocation technique are: 

(a) Elastically Clamped Edge? 
L (n(B,@,t)] = n(6,6,t) = s(t) , (3,15a) 
1 It was discovered at a late date that expressions for the 
in-plane displacements had been obtained for static loads by 
Reissner [97] and Fettahlioglu [98]. 
2 The edge conditions used for the stress function are 
equivalent to satisfying ug (5,0) =#0: butinot u_ (6,8) = a 


as will be explained in Section 6.1. 
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Li Cn(6,8,t)] = L(n(6,6,t)] - BL [n(5,0,t)] = 


une an(p,6,t) e * 
= Mg (Ps %st) - 8, eee = 67(t) - Rie toe eat ty 
(3.15b) 
eh _~ me 
L [&(6,6,t)] Sere. O sv) a= f(t), (37150) 
VS 
L C6(5,6,t)] = € (Ps 0st) meredt):, (3.15d) 
(b) Elastically Clamped Edge on Rollers 
ws 1 
Iinte,o,tideee (t), (3.15a) 
¥ me a 
L'tn(6,8,t) ] et): (3.15b) 
L,(E(6,6,t)] = FL (8,6,t) = F(t), (3.16a) 
L*CE(6,6,t) J = L Ce(5,8,t)] a BL Ce(5,8,t) J =e 
2 
=P -(pgegtvepeers Gpyespt) = f(t) + 
vs ss 
-3 2 
bem, fei ( i) (356 b) 


In Eqs. (3.15a) through (3.15d), (3.16a) and (3.16b) 

Lise is the flexural stress couple vector tangential to 
the boundary, ae is the in-plane shear stress resultant, 
Foy is the in-plane normal stress resultant, Fey LS tie 
in-plane normal strain tangential to the boundary and 

6 and 8 are the coordinates of boundary points. The 


boundary and interior coordinate systems are shown in 
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Wigyee.e, son Hay (3 25b) 6, "ranges from*zero*for"a 
shell unrestricted against edge rotation to infinity 
for a shell completely restricted against edge rotation. 
meds 0b) B, ranges from zero for a shell 
unrestricted against "edge extension to infinity for a 
shell completely restricted against edge extension, The 
operators L, to Lg are linear homogeneous operators of 
order two or lower, 

Expressing the functions in Eqs. (3.15b) through 


(3.15d), (3.16a) and (3.16b) in polar coordinates yields 
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Some of the stress resultant and stress couple components 
are shown in Fig. 2.3. The other vectorial components of 
the stress resultants and strains use the same kinematic 


subscript notation indicated in this figure. 
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For plane stress, strain can be expressed as 
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Stress resultants and stress couples can be 
expressed in terms of the dimensionless transverse dis- 


placement and stress functions as 








97F 
66 30? 
1 9E i et 
F = Ehl—-—+— —l], (3519b) 
Es Pa Og wpe 10.07 
9 1 0& 
F =F =-Eh| —/{|—-—II, (Bao cy 
ré er 9p \p 96 
Df d7n eno 27 Laon 
M = - -| — ye Se Ss GS2lod) 
re a ldap? eee oe 0 90 
D.Gat 377 ae on 97n 
M = -| = — + —- — + V — ~ (3.19e) 
Br sa tape 36° 0 30 9p? 
D Pp) 1 an 
Na lama =— (lev)| — |—- —| : (3.19f) 
rr 80 a 9p \p 386 


Substituting Eqs. (3.17a) through (3.17e), (3.18a) 


through (3.18c) and (3.19a) through (3.19f) into Eqs. 
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es inee On| C3r20f) 


Substituting Eqs. (3.3a) and (3.3b) into Eqs. 


(3.15a) through (3.15d) and setting 


L C#(5,8)] = 1, L* Cz (6,8) ] =eo) (3.21a) 


L,(# (B,6)] = 0, Late (B,6)] = 1, (3,21b) 
L,Ch (6,6)] = 1, uh'(,8)] = 0, (3,210) 
L th (6,6)] = 0, Ch (8,8)] = 21, (3.214) 


yields homogeneous "elastically clamped edge" conditions 


leg a and aA 


L,{n,(6,6)] = 0, (3,22a) 
L[n (6,6) = 0, Groep 
L Ce, (6,8) ] = 0, (3,22c) 
L Ce (6,9)1 = 0, (3.224) 


Similarly, substituting Eqs. (3.3a) and (3.3b) 
ipo Bqse\s.15a), (3.25b), (3.16a) and (3.165) and setting 


LCE (B,6)] = 2, LELE'(5,8)1 = 0, (3.21a) 
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L(87(6,6)] = 0, 4% (#2(3,0)] = 1, (3,210) 
L ({ht($,6)] = 1, Lith*(6,6)] = 0, (3.23a) 
L Ch?(8,6)] = 0, L*Ch*(6,0)] = 1, (3.23) 
yields homogeneous "elastically clamped edge on 
rollers" conditions in n and E 
Ltn (6,6)] = 0, (3.224) 
L*Cn (,8)] = 0, (3.22b) 
L Ce (6,6)] = 0, (3. 24a) 
L'CE (6,0)] = 0. (3.24b) 
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3.2.3 Auxiliary Edge Conditions, The choice 
-& - 
of functions g (9,8) and h*(p,6) which satisfy the 


auxiliary edge conditions formulated in Eqs, (3.2la) 
through (3.21d), (3.234) and (3:23b) is not unique, as 


Meirovitch [25] points out. It is convenient to choose 
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pn 
x vs a . 
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Substitution of Eqs. (3.26a)4through (3.26f) into 
Eqs. G6e2la) through €3F210)s"(3. 23a) "and 3.i23b). "yields 
linear nonhomogeneous equations, the solution of which 
yields the functions £(p,9), &7(p,0), h'(p,0) and h?(p,6) 
which are required for the time-dependent edge conditions 
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[¢ 8) 
n,(,0st) = 2 A*(p,0) T(t), (3.27a) 


E*(o,0) T(t), (3.27b) 


E(9,8,t) 
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The term v'z*(p,0) vanishes in Eq. (3.29) since B* (0,6), 
Eq. (3.25a), has been chosen as a solution of the 
biharmonic equation, Multiplying both sides of Eq. 
reo Ja Dy AY (p50), integrating over the characteristic 
shell segment area and applying the condition of 


orthogonality of the eigenvectors shown in Appendix B,2, 


yields 
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Using the restriction applicable to shallow 
shells, that fi-p?(2)*21 forQ0< p < a, permits the 
integration to be carried out over the projected area 
of the shell. Thus modal participation functions can be 
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products of functions of different orders n are neglected 
in the integration. This property is true only for 
shells with circular boundaries; The errors introduced 
for shells with polygonal boundaries should increase, 
Just as for) plates, with the deviation of the nodal lines 
from circles concentric about the apex. This deviation 
of the nodal lines from the circular shape generally 
increases with decreasing number of boundary sides and 
increasing mode number, The resulting modal partici- 


pation functions are 
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q(e) cos. pm6 = 1, and the coefficients deeds are: 


ad 











f i £ ha 
, + a ; 
+19 —. oe " er” o 
f B ££. ee 
+ hie ae" aa Di WOE 
I M4 ae 
+ 6 SS @ a] a (L * aq)’ } 
cs 6nd ie nq {sr 
L 2s & 
> > al “ais + aga + | i) 9 + ’ 
as 4 es 
nq rq Nq 
- tI bs 
BEE 2 + * 9 3 
— . ti es tt os” rr 
i - 
z,¢ i 4 f at 5,8 
: Fai Lec? ce ee ae gk A] = oF 
ft 
’ re? Ph , ne” . rs ea” igs : 
L M £ ry S 
+f o OS © @+ @ Al } a 
“*e ng as ag $2 ong ag 4f<n. 
te i 7: t @24e- 
(D¥ELE) {| o a + ® fF + 6 a} a + 
88 Be "Song as” nq 


baot edd 10 daebaegeh 915 | 4... > ednotottsos edt. 
7 — 

betudiadelrh baof s =o% .(8,9)5 aottenut aolsudia 

39 nouruin Iasi Lieds ont tevo 


2 





67 





























Q6 
cos? fe fi A fs dé 
PreGioeres 4 cu, | 
a 7 8 cos?0.- 
0 
84 
cos nth | ; . dé 
ee) at i ) | ; 
Sapa oes cos*é 
cars 1 tee O aoe 
2 O 
at 
) 
, 08 8, i dé 
ois J I + ] 
Ki | Gee Gy + 1, Gg), Gey, 
0 
6 
a dé 
ee COS 0 i ay Mes 
Pp cos 6 
@) 
86 
ous, 
oeA COs. 4” | iL Cu.) . 
sh 0 cos 6 
6 
O 
cos76 pn \? ; 
K 2 oy aig) #[2 = [=] Jor “a0 
a 2 pn u pn 
0 Oo 
cos’*pné 
soe Gy, 
cos78é 
s] 
cos’ 6) : ; pny 
ep eee ff ott Mig) peat || Pe AG 
a = pn u pn .o 
O 
0 
cos*pné 
db 
cos760 
8, 
2pnt+2 cos*pné 
‘ape MGengh ao 1d) | Gately she) i ares es 
+ © (cos g)tPnt 












7 











¥9 | - th) 
oe 7 
| ae a ofa yo : 
; 6teae Lua ? (ut ht o 5 i} 
e069 F 8 } 7 al _ 7 
° 
Ob @* x00 ; 
ft. u)eT + fa, nyt - ] TE 9 
1* 200 b a Br 
i 
©, s 
« ,o5 site — . Ps 
6” 
8b n° 209 
‘ [(a? %. (0) + 7 AH? 1\ ae 
§ #00 : 
’ oh 
i‘ Cu) % 8 k 46 ¢ 
Reverie ee ui? 7 BOD --=s 
‘¢ ap ° 4 \ . tt . 





f) 
ob | 
‘ a). I \ _& soo k 14 «= ,? 








a 802 9 tt 
9 
i 
£ . ELT £ 0 *a09 
x Cu) | oe rl Cw 'L \2 ‘y = % 
mq ot rg S i i 
0 _ 
ong*s0o 
2D ——inme 
OF eGR ae 









[=|+«]- (yin)? ts anal alae 
taht ie) es ys 


10 


16 


+ 


68 











s) 
; cos on 3 
- x’ HT eleet Gl) Et (ig) qb iaug dea ety) | x 
2h 2 . pn pn pn pn 
cos*pné 
dé, 
COSmd 
86 
2 
pnt cos “pné 
setkiaAeCon—+ locos 10) iI J, (tae) ay 
ti. 9 ; pnt, (cos 9)?” 
86 
; pntl cos*pnéd 
kA Gon. 1)) (cos 76, ) i i . 
+1 
4 ai 2 f pn+i1 ° (cos 6)?” 
OS 
dé 
- k.cos 6 if ee Gils) 
i % - gos 90 ; 
0 
8 
“ de 
k= cos, 0 | I (uo) ; 
x 9 : cos 6 
0 
Ay UST.) &, 
i O 
8 
© Z cos pné 
- 2k cos 6 ef as (u.)-__——_ de + 
ai © a=0 A pnt2zati ° cos 6 


8 
os O 
4 y. (a + 1) fl Ail (Lt diGO.Se Dio Or. 
a=0 ; pnt2zat2 ° 


de, 


dé, 













8s 


t [CoH Se ee 2 Goi) Ei) “ie 








Oaq*soo, 
Pats) ere ete 
6 #00 
o® 
6ag* BOO t+nq ba 7 
6b t=" (uy) u ( 86 goo) (f + aq) AaB ew o 7 
; "Ae goo) ° teat tt A ? if *s ) 
9? 
2 
Gitq" B00 ioRG ‘ 
9b. 25 ( g) I \ {8 goo) (£L + ag) A dh @ ¢ 
~ "Ae ao9) © ang a rt sf: 
0 i 
R ; 
metas ; = 6 
yw) | \ 6 soo Ae 6 
“8 800°! = t os 
0 
8 
80 aT F 5 ‘y 
——aee ? T n 
“2 ee eo aary ae 
- a 


«08 ante ke 7.7 
° t ‘@ 


> 
vos 


a 
6 goo WS = = 
° ); "40 gana 

oa 


(re of go 
; ; ve 


17 


23 


24 


af 


85 
CO 
4 YY (-1)% (a + 1) if 
a=0 0 


DuAtaL 


| 
Lag Dine t ec 


3 
kK cos ye 


aes 
ei 


=yke (cos. 6.) 
A O 


aL: 


ab 





O 
2pnt2 
2d (cos SA } 


0 


86 


© fo} 
CyeecOs"O me) 1(=1) I I 
O a=0 


pntzar2 


cos pné 
de 





Oe) 


pnt2aty cos 6 


(uo) COS pne edo, 


6 


O 
pnt2 cos pné 
(cos ha \ aero = are a oe 








2 
a sin 20, (1 + 2 cos*@,) , 


shake wel 
apn cos .c..) / i Gime) 
‘ £ pntlees 


(cos 9)Pnt? 


69 








86 
pnti J cos*pné 
1 UP Nt reemenmnmrceresteny ielkc).5 
; pnt, (cos 6)” 
8 
0 cos*pné@ 
ntl do, 
; (cos 0)? 
s) 
cos*pné 
dé 
liecs g)2pnt2 f 


0 


ck 86 
+ 
1 Yohe y, 2 cos*@, ice ae 
0 0 
@ 
° 
de ; © do 
=ONcos- 6 Le Ges) 
cos*6 s ip pa acos.e 


eo 











o® 
Onq Bac . 
me eneeent 5 tf. 


Ae ; 
Cay) t \ (OG 
@ soo 9 ¢kwsterg \ Oo 


oe 
,05 Omq Boo Cy). 


7 _ 
I \ (io + oh en7 a¢ z 
0 s#0 strig. ‘ 











8 
Srtq. BO. sng {+ nq . 
o——— ( & gos) (= kt = Vi 
S*19(6 200) o $ + aqi ft 
0 
hy 2 
€b . 6b 
e (5 L \ 9" 809 S$. + (oH). % \ dfeos # = = ast 
6* 00 . WT 6* Bo9 : " tL i 
58 of 
8h *. : 0b 58 a F 
ewe ( i). IT 6*s09 S «= (.u) I aco 4A =. 8 
; 8% 800. © \ 5* 805 ” 9 f oe 
; 
s rh 
« (58*e00 S + £) 20S mie ae «? : 
9 ' 
Gnq*eon °,  .#nag ; 
: 6b reo) \ ( 6 amo) a + 8 
\ ‘ | 
‘ aT 208) Oia : ; #4 





8 
Cha* 209 





70 


8 
pnt3 ° cos*pné 
W) = - k (cos 0 ) f J GUE 3) o Seteeremmre gem els) Fy 
< 1 e 0 pnt: © (cos 6)Pn*® 
re) 
pnt2 cos *pné 
tee CCOS Oe) f J re) apni eoneoremmnenirry Ae heh 
. 0 pnt2 (cos g)Ph 
s) 
pnts = cos*pné 
d =k (cos 6 ) | ah G0) eee 
aie 4: P 0 pnti © (cos 9)?” 
8 
pnte cos7pné 
- 2(cos 6) f I (le) amet naenyins a0 
es 0 Diteee oe COs 0). 
£5 
pn+i epnty cos*pnée 
¢ = 21 |———J| (cos 9 ) i a ae 
aed ilpn + 2 Oo (cos 9)?Pnt* 








For a concentrated load P = Piatt) at the shell 
apex, Po replaces oo Ere O.s (3.56), ana themcoc: cients 


a 
for x become 
1 


¢ =-o = - k*/(2pa?), 
14 i 


o- 
W 


Wreytepa oe, 
nt 


Oy 











> 
on — 


bat* eon ; : \ — o-nak) Mee : 
+ OD opener . fap: 0 8G2)' > ea? 
“TAGs eos) © rng. g ° r = sen 


i] 





: a 7 
ang" 209 . gta = ae 7 
ry a em ee L 
2° i «Oe ; i? Seavey - 
(8 p03) Siig 0 >: i. 
°° 
diig*eox. * ia on a ’ 
Ot Cane. ( u) L \ 8 e0g : = ‘ 
"T(G edo) 9 t4aq ; Ean tf 
) 
Oi2* aoo stg 
«9D a Fan (ww) ‘(fT \ (9 200) = . 
“(8 @05) % iShng y al : 
9 ' 
Onq* ade i. ets i + ag 
~ \ 9 Bow) ——| AS = sf 
“TAT (O Boo) =f e S + agit . 
' 0 - 
iferie off! +8 (71004 = 7 bsol bodsisaesaes & 10% 7 


. 
alielollisoos sid Gna’ ($82.6) sod at 3? aesoalqet o «*9q8 


emooed x sot 
. | e aS 
4S ane Say = * re : 
me a oi? nal aia | ; 4 


ila al = a 





(es 


CHAPTER 4 


PLATE NUMERICAL RESULTS 


4.1 Higenvalues, Eigenvectors and Modal Participation 


Functions 

in the practical execution orstne.poundary 
Gollocav@on technique, itfis necessary to assume 
convergence of the series solution for the transverse 
displacement eigenvectors, Eq. (2.13). Using this 
assumption, symmetric eigenvalues for various clamped 
and simply supported polygonal plates were obtained by 
employing machine computation using 16 figure accuracy 
Gomsearcn Loreche roots of a series of linear homo-— 
geneous equations of the form of Eqs. (2,10a*) and 
(2,.10c*) where the summation with respect ton is 
truncated, 

Collocation points were regularly spaced on 
the boundary of each plate. A Cypical seu of nine 
boundary collocation points for a square plate is shown 
lieth ree cond. td On) gS. Go, dUas) send ACe Uce® 
Werewsavisticd. at pounce. 1.5). VY. 5.016) vendo oe Gh 
(2,10c*) was satisfied at point 2 and Eq, (2,10a*) was 


satisfied at point 9. Edge condition Eq, (2.10c*) was 
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@ POINT ON CHARACTERISTIC PLATE SEGMENT AT WHICH BOUNDARY 
CONDITIONS ARE SATISFIED IN NUMERICAL SOLUTION. 


O POINT AT WHICH BOUNDARY CONDITIONS ARE AUTOMATICALLY 
SATISFIED DUE TO PLATES ROTATIONAL PERIODICITY. 


FIG. 4.1 Boundary Collocation for a Square Plate 
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not satisfied at point 9 since in practice it is often 
difficult to predict its value at the plate corners, 
Edge conditions were satisfied in a similar manner for 
the other plates. Consequently, an equal number of 
Simultaneous boundary equations and unknown integration 
constants was obtained with L collocation points on the 
boundary of the characteristic plate segment for n 
ranging from 0 to (L-2). 

The convergence of the symmetric eigenvalues 
k, to the values given in Tables 4.1 and 4,2 was 
checked for each polygonal plate by increasing the 
number of collocation points used in the truncated series 
solution until the minimum number required to obtain 
these values was determined, It was found that funda- 
mental eigenvalues for simply supported plates with from 
7 to 12 sides as well as some of the higher modal eigen- 
VemUeowmespeclally Tor 98and gl2psides, qwere less 
convergent than the eigenvalues obtained for other plates, 
sometimes having significant variations in the third 
figure. Errors in these eigenvalues might have been 
caused by the Poisson's ratio effect which is explained in 
Section 4.2. Fundamental eigenvalues previously obtained 
by Conway [35], Leissa [36] and Conway and Farnham [37] 
for some of the plates using the boundary collocation 


technique agree to two or three figure accuracy with the 
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Drecehusvaluecs. Ao vne number of plate sides increases, 
these eigenvalues converge to the values given by Flynn [16] 
and Weiner [19] for circular plates. The eigenvalues given 
in Table 4,.3b for a clamped square plate all lie within the 
previous best limits given by Bazley, Fox and Stadter [102], 
and the eigenvalues given in Tables 4,4a and 4.4b for simply 
Supported triangular and square plates are accurate according 
to the exact solutions [103,21]. 

The values of the modal participation functions 
xi 7x7 for uniformly distributed loads and oe? Be for central 
point loads, given in Tables 4,1 and 4,2, were obtained using 
the eigenvalues in the same tables and their associated 
eigenvectors where ae = 1, Integral values for the coefficients 
ore! Fe were determined using the generalized Simpson's 
formula with 56 equal increments of 9 for the equilateral 
triangle, 46 increments for the square and 36 for the other 
polygons. The summation with respect to a was truncated when 


the values—-of and ees Converced vor ti ves del ceaccuracy.: 


The eigenvectors At associated with the eigenvalues 
for clamped plates in Table 4.1 are given in Tables 4.3a to 
4,3h and the eigenvectors associated with the eigenvalues for 
simply supported plates in Table 4,2 are given in Tables 44a 
to 4.4h, The elements ie and be of the eigenvectors 
converged to the accuracy shown in these tables for the number 
Of collocation points indicated. Some of the higher order 


elements which are given to only one significant figure may 


be accurate only to an order of magnitude. 
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Eigenvalues, Modal Participation Functions and Eigenvectors for 
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Table 4.3a 
Eigenvalues, Modal Participation Functions and Eigenvectors for a Clamped Plate 
With 3 Sides 
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Table 4.3b 


With 4 Sides 
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-0.2224 E-03 


p20" 305) AO, 250 4055, 60 


-0.6E-06 
-0.627E-03 


0.100000E 


mamamamMmmMmmmm 


0.100000E 01 
0.13977 E 01 
E-02 
E 01 
E 01 
E 02 
E 03 


E-03 
E-03 
E-~02 
E-01 
E01 
E 02 
E 03 
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-0.2E-07 
0.0178 


0'.100000E 01 
E 00 
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Wr 
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1 
CO000000 
eG See is 
OW ADWNHNNDOD = 





a Clamped Plate 


0.100000 
-0.1118 
-0.7869 

0.193 
-0.9 
=O. 


-0.9 


0.156 

0.939 

0.134 
-0.106 
-0.341 
-0.8 
-0.1 
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Table 4.3c 


Eigenvalues, Modal Participation Functions and Eigenvectors for a Clamped Plate 
With 5 Sides 


W Collocetion.Polnterst G7) 5°4010°,820°,..26°,. 347, 36° 


3.78397 


1.61783 
1.6288 


0.100000E 01 
-0.39435 E 01 
E03 
E 06 
Eat 
E 16 


68067E-01 
6549 E 01 
97, a 02 
E 06 
ali 
Es6 


7.39506 


-1.2945 
-1.369 


0.100000E 01 
0.42284 E 00 
-0.44603 E 01 
0.17 E03 
-0.8 E 05 
-0.8 E 09 


-0.276570E-02 
-0.12898 E-01 


10.5049 


0.100000E 01 
-0.1422 E 01 
O.2052° 6 07 
-0.13 E 02 
0.1 E O04 
0.5 E 06 


0.12218 E-03 
-0.646  E-03 
-0.406 €E-01 
-0.13 E 01 
-0.9 E 02 
-0.2 E 06 


Table 4.3d 


With 6 Sides 


11.5889 


0.3780 
1.202 


0.100000E 01 
0.15750 E 01 
-0.72990 E 00 
-0.13 E07 
-0.6 E 02 
-0.6 E 04 


0.10561 E-03 
E-03 
E-02 
E 00 
E 02 
E O4 


F Col Location: Points at.0°; 5°, 10°; 15°,.20°,; 25°) 30° 


3.57640 


1.61710 
1.6198 


0.100000E 01 
-0.11747 E 02 
E 05 
E 10 
E 16 
E23 


ooo 
. 8 @ 


1 1 
bboboo 

: 

nD 


7.0267 


-1.3048 
-1.323 


0.100000E 01 
0.4811 E 00 
-0.338 E 02 
0.5 E 05 
-0.8 E 09 
-0.1 


-0. 26439 
-0.25649 
0.486 


-0 
0 
0 


10.3921 


0.100000E 01 
=O). 472) 8 EL00 
0.2208 E 01 
-0.19 EsO3 
0.3 E 06 
0.2 E 10 


0.12418 E-03 
0.205 E-03 
-0.542  E-01 
0.3 E 01 
-0.2 E05 
-0.6 E 09 


11.857 


0.02880 
0.440 


0.100000E 01 
0.3909 E 01 
-0.397 £01 
E 02 
E 05 
E 09 


E-04 
E-03 
E-01 
E 01 
E O4 
E 08 


-0.6005 
-1.083 


0.100000E 01 
-0.13175 E 01 
-0.54046 E 00 

0.290 E01 
-0.24 E 02 
-0.7 E 03 


-0.7175 
0.5610 
0.1907. E-03 

-0.141 E-01 

—Oeilitvae Oil 
0.4 E 02 


E-05 
E-05 





Eigenvalues, Modal Participation Functions and Eigenvectors for a Clamped Plate 


-0.88766 
-1.097 


0.100000E 01 
-0.634 E 00 
-0.206 €£ 00 

0.61 E 01 
-0.8 E 03 
-0.4 E 06 


E-05 
E-O4 
E-O4 
E-01 
E 01 
E 05 


-0.5952 
-0.1488 
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Table 4.3e 


Eigenvalues, Modal Participation Functions and Eigenvectors for a Clamped Plate 
With 7 Sides 


7 Collocation Points at 0°, 5°, 10°, 15°, 19°, 22.5°, 25.7° 


1 2 3 4 5 
k; 3.4632 6.8181 10.16 12.621 13.646 
xipgi n=O, 4.8 1.6160 1.3007 1.099 -0.107 0.853 
OY hao 1.6169 -1.3067 4.125 -0.684 -1.028 


0.100000E 01 0.100000E 01 0.100000E 01 0.100000E 01 
E02 Oo 33s OOM =O cr mm OO E 01 
EO7 -0.46 E 03 0.69 2 (eh E 01 
BWA 0.13 Eos -0 E O04 E 03 
Br2250 -0.1 E 14 0. E 09 Clee 6 504 
E30 -0.4 Ee20 0 = WA ° Eee 


E-O1 -0.2593 E-02 E-03 . E-07 
E02 =Os579 = e=Oi E-03 E-03 
E07 0.95 E 02 E 00 . E~01 
EMoe 041 E 07 2 8) E 00 
B22 0.6 Eas E 08 E 06 
E 30 0.2 E20 SA End 


0 
1 
2 
3 
4 
5 
0 
1 
2 
3 
4 
5 





Table 4.3f 


Eigenvalues, Modal Participation Functions and Eigenvectors for a Clamped Plate 
With 8 Sides 


6 Collocation Points at.0°, 5°,. 10°, 18°, 19s 2255 — 


1.6154 -1.297 
1.6158 -1.300 


0.100000E 01 0.100000E 01 0.100000E 01 0.100000E 01 
-0.1816 E 03 O.13770NEnO1 69-0.220 £400 E0O -0.297 €£ 01 
0.29 Eo9 -0.91 E 04 0.41 E02 E 01 0.51 E01 
-0.71 Sali 0.8 E10 -0.1 ESO”, EO4 -0.2 E 04 
-0.3 Ey 0.1 E18 -0.4 Beii2 E09 -0.4 E 08 


0.5592 E-01 -0.2568 E-02 0.1150 £E-03 E-05 -0.410 E-05 
0.1025 £03 -0.1465 E 00 0.132 £-02 E-0o4 -0.109 E-03 
-0.20 E 09 0.24 EOL -0.23 E 01 E-0O1 -0.320 E-01 
0.6 E17 -=053 E 10 0.1 E 06 E 03 0.2 E 02 
0.2 E27 -0.6 EZ 0.1 iS We E 08 0.3 E 07 


0 
1 
2 
3 
4 
(6) 
1 
2 
3 
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Table 4.3g 


Eigenvalues, Modal Participation Functions and Eigenvectors for a Clamped Plate 
With 9 Sides 


ae (Oh, Gas 
X4/%o 
n=0 


rwWwn-oO P~WNH-O 


& 


6 Collocation Points at 0°, 5°, 9°, 13°, 17°, 20° 


1 


3.349 


1 2615 
1.615 


0.100000E 01 
-0.891 E 03 
0.36 Eds 
-0.4 Bea 
-0.7 B32 


E-01 
E 03 
E 11 
E 21 
E932 


0.5585 
0.53 
-0.27 
0.3 
0.6 


2 
6.603 


-1.295 
-1.297 


0.100000E 01 
EO 
E 06 
= le) 
E 22 


E-02 
E 00 
E105 
Eel3 
E22 


0.100000E 01 
-0.285 E 00 
E 03 
E 09 
E 16 


E-03 
E-~02 
E 02 
E 08 
Ss 16 


Table 4.3h 


13.12 


-0.956 
-0.979 


0.100000E 01 
0.2574  Ea0e 
-0.95 E 01 
0.2 E 06 
0.8 ral 


-0.504 
-0.266 
0.118 
-0.8 
-0.1 


E-05 
E-04 
E 00 
E 04 
E11 


14.71 


0.0006 
-0.04 


0.100000E 01 
-0.726 £01 
0.249 E02 
-0.6 E 05 
-0.1 ala 


-0.108 

-0.126 

-0.124 
0.1 
0.7 


E-05 
E-03 
E 00 
E 04 
E 09 





Eigenvalues, Modal Participation Functions and Eigenvectors for a Clamped Plate 
With 12 Sides 


4, Collocation Points at harp 10°; US 


0.100000E 01 
-0.21 E 06 
0.3 E 18 


E-01 
E 06 
E 18 


0.5574 
0.14 
-0.24 


-1.292 
-1.293 


0.100000E 01 
0.77 E 02 
-0.2 E 11 


-0.2536 £-02 


E 02 
ait 


-0.164 
0.1 


0.100000E 01 
-0.140 E 01 
0.2 Es 07 


0.1115 
0.432 
-0.3 


E-03 
E-01 
E 06 


te) 


0.100000E 01 
0.194 — 00 
-0.5 E 04 


-0.487 
-0.346 
0.2 


E-05 
E-03 
E803 


0.100000E 01 
-0.173 E 00 
0.7 02 


E-06 
E-05 
E 00 


0.212 
0.33 
-0.4 
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Table 4.40 


With 3 Sides 


10 Collocetion Points at 0°, 5°, 10°, 20°, 30°, 40°, 45°, 50°, 55°, 60° 


1 2 3 4 5 


k; 4.18879 8.37758 114 .08250 12.56637 15.10290 
yg tire ee 1.3900 0.79467 0.000031 -0.0000011 0.00002 
Teens 1.2134 0.50687 -0.2522 -0.39 0.2993 


Or 

0 0.100000E 01 0.100000E 01 0.100000E 01 0.100000E 01 0.100000E 01 

1 -0.115470E 01 0.115470E 01 -0.623483E 00 -0.1059 E 07 0.862226E 00 

2 -0.19997 E 01. -0.199999E 01 0.833820E 00 0.7274 £00 -0.230316E 00 

3 0.110 E01 -0.115475E 01 -0.114334E 01 0.1059 E 07 0.663646E 00 

ai 4 0.9 E 01 0.20007 E01 -0.13048 E 01 -0.97 E01 -0.194695E 01 
pn 5 -0.2 E 05 0.1 E01 0.5908 Ee s00R 0.11068 EnO7 meu—Oeg O51 ceOn 

6 =0.24 E 08 0.2 El02m 0.20 E O10. 4, E 04 0.679 E00 

7 -0.5 EB it 0.5 EO, ae =O65 E 01 0.9 E06 -0.43 E 00 

i} -0.8 Er 0.6 EvO6mm=Os2) E03 =0.1 E 08 0.19 E 01 

0 0.1 E-05 0.5 E-08 -0.34 E-09 -0.55 c-04 0.8 E~11 

1 -0.13 E-04 -0.2 E-07 0.13 E-08 0.19 E-03  —022 E-10 

2 OnL7 E-03 0.2 E=Oomer= On, E-08 -0.8 E-03 0.8 E-10 

gi 3 =Oe47 E=01) 0-4. E-05 Onn E-07 Our E02 On E-09 
pn ii 0.2 E 02 0.2 E=03e O62 E=O5mueO nl E 00 0.6 E-08 

5 0.3 E03 -0.5 E-01 O02 E-03 0.8 EOlmeOns E-06 

6 0.3 Eos -0.9 = (oy 0.18 E201 0.5 E03 -0.1 E-04 

7 0.8 Elite =08d E 04 0.1 E 01 0.3 E05 -0.4 E-03 

8 0.2 E 15 0.4 EXOGme-O 62 E03 -0.3 E 07 0.2 E-01 

Table 4.4b 


Eigenvalues, Modal Participation Functions and Eigenvectors for a Simply Supported Plate 


Eigenvalues, Modal Participation Functions and Eigenvectors for a Simply Supported Plate 
With 4 Sides 


$ Collocation Points tat:0°, 5° 10°, 20°,.30°, 35°, 40°, 45° 





(0) 
1 
2 
3 
4 
5 
6 
0 
1 
2 
3 
4 
5 
6 


3.14159 


1.58955 


0.100000E 01 
-0.200003E 01 

-2009 E 01 
E 02 
E 06 
E 10 
B14 


fo) 
RUSH AROnD OND 


E-06 
E-05 
E-02 
E 02 
E 06 
E 10 
E 14 


b0b0b085 obo 


7.02482 


-1.0564 
-0.87236 


0.100000E 01 
0.559996E 00 
-0.168636E 01 
-0.1507 £ 01 
0.16 E 01 


Oooo0°o oo 
° 


SRORwWUmD Wo 


Le 


9.42478 


0.17677 
0.2728 


0.100000E 01 
-0.200001E 01 
0.199999E 01 
-0.199982E 01 
0.199 E01 
Bp E 01 
ii E 03 
E-08 
E-08 
é E-06 
4 E-05 
E-03 
E 00 


0 
-0 
0 
-0 
0 
-0 
0 
-0 
0 Broz 


11.32717 


0.60157 
0.4905 


0.100000E 01 
0.140828E 01 
~0.167371E-01 
-0.143184E 01 
-0.200006E 01 

E 01 

E 01 


E-10 
E-09 
E-08 
E06 
E-05 
E-02 
E-02 


12.95312 


-0.24331 
0.0316 


0.100000E 01 
-0.111418E 01 
-0.758587E 00 

0.195940E 01 
-0.142468E 01 

E 00 
E 01 


E-10 
E409 
E-09 
E-07 
E06 
E-03 
E-02 
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Table 4.4c 


Eigenvalues, Modal Participation Functions and Eigenvectors for a Simply Supported Plate 
With 5 Sides 


7 Collocation Points at 0°, 6°, 12°, 18°, 24°, 30°, 36° 


1 2 3 4 D 


2.8212 6.4282 9.5642 10.65 13.0870 
Ts ee 1.6049 =10777 0.514 0.3535 -0.5043 
Xie n=O 1.5776 -0.9900 0.7033 0.5678 -0.509 
Order_n 

0 0.100000E 01  0.100000E 01 0.100000E 01 0.100000E 01 0.100000E 01 

1 SOLS17~2 E01 ~Oc4259° E00 <0.13472°E 01 0:164008F Of ~-0.12864°E 01 

ai 2 0.17 EO, -0.557 £01 0.2051 E01 ~-0.80032 £00 -0.517632E 00 

pn 3 -0.4 Bos "0.66 9 F 03 20476 E02 =0:131. E01 - 0.3903 £01 

2 0.1 E43 “co.2 E06 0.1 EOL. “OW E02 0.2 E 02 

5 -0.8 E19 0.2 E14 «£06 EO7 “0; E06 “OH E 05 

0 0.6 etn) ehors7 “BESS 30567 ESR 0r39 «E07 | 0 E-07 

1 On E-02 0.5 BAO, -fohi4 “ES05: FOlS. “EEO: 0545 Eee 

B 2 0.2 E03 " Gols Eso1 fo 33 E-03 -0.1 EO, <0.4 E-05 

Pr 3 0.6 E07 ~-0.2 Eton =fo73 E-01 -0.6 E-03 -0.5 E-o4 

4 0.6 Ei 3 065 E06 O04 E403 =0.7 Eo1 -0.6 E 00 

5 0.1 E16 uls-0),4 E09 0,3 E06 1-034 EOL) 0,3 E 03 





Table 4.4d 


Eigenvalues, Modal Participation Functions and Eigenvectors for a Simply Supported Plate 
With 6 Sides 


6 Collocation Points at 0°, 6°, 12°, 18°, pI sie 


-1.077 
-1.032 


0.100000E 01 
0.49 E 00 
0.5 E02 
-0.1 E 05 
011, GEN 


-100000E 01 
E 00 
E 01 
= (os) 
E307 


0.100000E 01 
E 01 
E 01 
E 02 
E06 -! 


0.100000E 01 
-0.627 £E 00 
-0.197 E00 
0.54 E 01 

E 04 


° 
Wr 


OoOO000 
e e 


E-06 
E-O4 
E-03 
E 02 
E 06 


E-06 
E-05 
E-03 
Sake 
E 05 


E-04 
E-03 
E 01 
E 05 
E 09 


E-07 
E-06 
E-O4 
E-01 
E 02 


0) 
1 
2 
3} 
4 
0 
1 
2 
3 
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Table 4.4e 


Eigenvalues, Modal Participation Functions and Eigenvectors for a Simply Supported Plate 
With 7 Sides 


6 Collocation Points at 0°, 5°, 10°, 15°, 20°, 25.7° 


0.100000E 01 0.100000E 01 0.100000E 01 0.100000E 01 0.100000E 01 
-0.60 E 02 On E00 -0.260 €£ 00 Eo1 -0.684 € 00 
0.5 E08 -0.9 E 03 0.76 E 01 E01 0.68 E-01 
0.2 Bm6 Eos -0.1 E 04 E703 6) Emo? 
0.3 Bez E14 0.2 E 10 Eos -0O E 06 
8 E-03 E-O4 0) E-05 E~06 0 E-~07 
95 E 01 ° E-02 -0 E-04 E-O4 0 E-~06 
02 E 08 . E 02 0 E-01 E-02 -0O. E-O4 
3 E 16 Eos -0O E 04 E 02 0 EO 
4 E 24 iS, Ve 0 E 09 ° E 06 6) E105 
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Table 4.4f 


Eigenvalues, Modal Participation Functions and Eigenvectors for a Simply Supported Plate 
With 8 Sides 


6 Collocation Points st 0°, 5°, 10°, 15°, 19°, 22.5 
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Table 4.49 


Eigenvalues, Modal Participation Functions and Eigenvectors for a Simply Supported Plate 
With 9 Sides 


6 Collocation Points at 0°, 4°, 8°, 12°, 16°, 20° 
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Table 4.4h 
Eigenvalues, Modal Participation Functions and Eigenvectors for a Simply Supported Plate 
With 12 Sides 


4 Collocation Points at 0°, 8°, 12°, 15° 
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Values for modal participation functions O/x 
calculated for uniformly distributed loads using only the 
6- independent zero order terms in Eqs. (2.22) and (2.23) 
and henceforth referred to as (xi 7x") are included together 
with the values calculated using the higher order terms 
which will be referred to simply as eye The difference 
between the values of Syne and (een. iS 2nd Cavive vO 
some extent of the Slenificance of the errors introduced 
by the omission of the products of the terms of different 
orders n in er EQ. \e.ec)) Results in Tables 4,3” and 4.4 
show that the possibility of error in Me increases with 
increasing mode number and decreasing number of plate sides, 

The Fortran digital computer program used to 


obtain these results is given in Appendix C.l. 


4,2 Effect of Boundary Shape on Symmetric Eigenvalues 


The effect of boundary shape on the symmetric 
eigenvalues k©® for plates of equal surface areas is shown 
i 
in Tables 4.5 and 4.6, These eigenvalues are calculated 


from the eigenvalues k given in Tables 4.1 and 4.2, using 
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and o~ 18 the dimensionless radius of the circular plate 
having the same area as the regular polygonal plate, 

These tables reveal a pattern of fluctuation 
of Sy with the number of the sides of the polygon and the 
mode number which, with the exception of a few discre- 
panctes, appears to be identicalefor both simply 
Supported and clamped plates. The direction of decreasing 
elgenvalue magnitude for either increasing or decreasing 
number of sides is marked by arrows for each mode, These 
arrows indicate the polygon having the maximum modal 
frequency, 

Exceptions to this pattern occur in Table 4.6 
where the fundamental and some of the higher modal eigen- 
values appear to be too low for the simply supported plates 
with 8 and 9 sides. The same discrepancies also were 
noticed in-these eigenvalues when calculating modal 
Gisplacements witich aretdiscussed later. The magnitude of 
these discrepancies couldspe faccounted for by the numerical 
convergence error previously noted for the related eigen- 
values in Section 4.1, Table 4.2, 

Since the frequency varies as the square of the 
eigenvalue, these tables show that the fundamental fre- 
quency of vibration of polygons of equal area having simply 
Supported or clamped edges decreases with increasing number 
of sides. The relative difference between the fundamental 


eigenvalues for a circular plate and a plate with 12 
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equilateral sides is larger for the simply supported 
edge condition than the clamped edge condition. This 
seems reasonable since eigenvalues are independent of 
Poisson's ratio in the case of simply supported poly- 
gonal plates as well as clamped polygonal and circular 
plates, but are dependent on this ratio in the case of 
simply supported circular plates, having their smallest 
values for v = 0, Rao and Rajaiah [104] have discussed 
this-effect of Poisson's ratio for statically loaded 
Simply supported plates. They showed that inaccuracies 
can arise for p>5 on account of the slow rate of 


convergence of the solution, 


4,3 Symmetric Mode Shapes and Modal Participation Functions 


The symmetric mode shapes for the characteristic 
segments of some clamped and simply supported polygonal 
Ploves. computed using the values for the eigenvectors 
given in Tables 4,3 and 4.4, are shown in Figs. 4,2a 
through 4.2d. Mode shapes for the complete plates are 
inferred by the rotational periodicity of the solution. 
The boundary collocation points used to obtain the eigen- 
vectors are shown as small circles along the plate 
characteristic segment edge, 

The dotted portions of the riodal lines in Figs. 
4,2a through 4,2d are considered to result from errors in 
the satisfaction of the edge conditions, More errors 
occurred for the clamped plates than the simply supported 
plates since the satisfaction of the condition of zero 


edge slope naturally tends to introduce more incorrect 
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nodal lines near the boundary than the condition of zero 
edge moment. Since this solution,is exact only for 
circular plates and since a polygonal boundary deviates 
most from a circle’at the corners and this deviation 
inereeses with decreasing number of sides, dtlis 
consistent that most numerical errors occurred for the 
triangular plate near the boundary in the corners along 
generally circular lines emanating from boundary 
collocation points. 

The deviation of the nodal lines from concentric 
circles increases with decreasing number of sides and 
increasing mode numbers, The degree of this deviation 
is reflected in the degree of variation between the values 
of OG7x7) and be given in Tables 4.3 and 4,4, 

An attempt was made to relate the nodal areas, 
defined as the areas between the nodal lines and denoted 
as ratios of the plate segment area in Figs. 4.2a through 
4,2d, to the modal participation functions for uniform 


loads using 


(xt7yt) = om (271) , (4,3) 
1p oe 1 0 pe12 
where 
Mw at/iai) , (4,4) 
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FIG.4.2a- SYMMETRIC MODE SHAPES FOR THE CHARACTERISTIC SEGMENTS OF REGULAR 
POLYGONAL PLATES WITH 3 AND 4 SIDES 
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In Eq. 423) OG7x2) are the values of ‘Cae for the 
p=12 
12-sided polygonal plate given in Tables 4,.3h and 4.4h. 
af 


In Eq. (4.4) a is the difference between the sums of 
the i-th nodal areas which respond in the same direction 
under a uniform load. 

The modal participation functions xX for 
the 12-sided plate were chosen as the basis for obtaining 
OC/x,), for the other polygons since the evaluation of 
a from Eq. (2.22) is most accurate for the polygonal 
boundary most closely approximating a circle, 

The modal participation functions C/K) as 
calculated from Eq. (4.1) and given in Tables 4.7 and 4.8, 
agree reasonably well with XX with differences 
generally increasing with decreasing number of sides and 
increasing mode number. The largest discrepancies occur 
for the fourth mode of many of the polygons and for all 


modes of the triangular plates. 


4.4 Central Displacement 


The static central deflections for the symmetric 
modes of response of polygonal plates with the same 
a71I8cribing radius a, were calculated for uniformly 
distributed loads using 
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and ai i i 
a A +B Xs 
W els SS Or 0 e 5% 5) 9 (4.5b) 
sil WEEP geveyeven) Na 1a 
4 0 
and for central point loads using 
ak Ai i 
p ao 8 7 xa), 
W = ° ‘eas er e ao Ca o6) 
eeecos 0 X 
1 0 0 


The convergence of the sums of the modal 
displacements for polygons with clamped and simply 


Supported edges subjected to uniform pressures Sha 
W= )W = wD/q_a® (4,7a) 
roy: P xe 


and 


a a : 
we Dw (4.7) 
ak. (50 


to the values given by Leissa, Lo and Niedenfuhr Le eze 
is shown in Tables 4.9 and 4,10 for the first em 
modes. 

The agreement between the values of W, W@ 
and the values given in [77] is quite good for all 
polygons in Tables 4.9 and 4,10 except the triangular 
plates and the simply supported plates with eight and 
nine sides. The disagreement between central dis- 
placements calculated for these simply supported plates 
with eight and nine sides by the static approach in [77] 


and the present dynamic approach could be eliminated 
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by using the slightly larger fundamental eigenvalues 
indicated by the discrepancies previously noted in 

Table 4.6. Good agreement could have been obtained 
between the static central displacements computed by the 
two methods for the simply supported plates with seven, 
eights and jninewsides: 1f whe fundamental eigenvalues were, 
respectively, 2.59, 2.53 and 2.49 rather than the values 
2.60, 2.50 and 2.45 given in Table 4.2, Consequently, 

it seems that the modal participation functions for these 
Dlaleoesnould tervreltable.y However; for; the triangular 
plates the values of Xa“ appear to be incorrect while 
the values of OxE/%9 5 provide good convergence of the 
sum of static central modal displacements W* to the 
values given in [77]. 

From these results it appears that the values 
of GUARIE which were calculated from the nodal areas 
Provide |veryereliable estimates fof, at least. ther first two 
modal participation functions for the symmetric response 
of uniformly sloaded Gegular podygonal ae ore 

The convergence of the sums of the modal 
displacements for polygons with simply supported and 


clamped edges subjected to central point loads es 
te iy 
W = >) W = wD/P_a? (4,8) 
anes Saat 
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modes, As far as is known, static deflections for 
centrally loaded regular polygons have been calculated 
only for triangular, square and circular shapes, These 
values have been given by Timoshenko and Woinowsky- 
Krieger [78]. Agreement with present values is poor 

for the simply supported triangular plate and fair for 
the square plate reaffirming that errors exist in 

values of x calculated from Eq, (2.22) for these plates. 
Better agreement for these plates could be obtained by 
using values of x for the first two symmetric modes 
which are derived from values of (xe/x) given in 
Tables 4.7 and 4.8 along with the values of ee 
calculated from Eq. (2.23). However, values of vs 
calculated in this way provide unacceptable values for 
the modal participation functions ere for centrally 
loaded triangular plates for i > 3. Values of WE should 
be reliable for the polygons with p > 4 except for the 
errors previously mentioned in the eigenvalues for the 
first mode of response of the simply spponeed plates 
with eight and nine sides. 

Comparison of the static central displacements 
for the first five symmetric modes for each plate in 
Tables 4.9 to 4,12 reveals that usually only the first 
two modes must be considered for computing the response 


of these plates to most uniformly distributed dynamic 
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loads while the first five modes may be needed for 
central dynamic point loads. Naturally, the response 
of the higher modes becomes more significant as the 
load duration becomes short compared to the lower modal 
response periods, 

Static central displacements of polygonal 
plates having the same surface area and subjected to 
uniform and central point loads, termed W® and Woe 
respectively have™been ‘calculated from the values 
given in [77], [78] and from present values where 


existing ones are not available, using 
Wwe = W cos"6 /(p°)" = wD/q.,(p©a) * (4,9a) 


and 


Pye [ 2 2 €.\2 
W = W cos 6/6 o*) = wD/P (pe ayn, (4.9b) 


Comparison of the values of W® and iets given in 

Tables 4.9 through 4.12 shows that central displace- 
ment increases with the number of polygon sides with 
the exception of the value of w'’° for the simply 
supported plate with nine sides which can be ignored 
for the reason given previously. This implies that 

for simply supported and clamped polygons of equal area 


the effective plate flexural bending stiffness increases 


as the number of boundary sides decreases, 
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The convergence of the central displacements 


We ana wrsé 


to the values for circular plates as the 
number of sides of the polygon increases is not as 
rapid for simply supported plates as it is for clamped 
plates... This can be explained partially_by~the-fact 
that the displacements W and wP are dependent on 
Poisson's ratio only for simply supported circular 
plates, decreasing with increasing values of v. 

The undamped forced motions of some clamped 
plates with the same circumscribing radius a and 
subjected to the initial conditions given in Eqs. (2.20a) 
and (2.20b) are depicted in Fig. 4.3 for a typical blast 
wave of incident load intensity Vo and duration of the 
positive loading phase toe These displacements were 
computed using the values of x7 /x" given in Tables 4,3 
and 4.4, The smoothness of the ERCOTStHM ET response 
Curves is indicative of the small finfluence Of the 
higher modes on displacement for this transient load. 
However, the influence of the higher modes on stress 
couple response is more significant, as is shown in 


Chapter © for a square plate, 
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UNIFORM EXPONENTIAL 





te BLAST LOADING 

m= Gi 76-426 Aaa 

a = 0.2310 oh aitl=q(i-t)e '/%o 
t, t 









l2 sides 
6 sides 
“sides 
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FIG. 4.3 Theoretical Undamped Center Deflections of Some Clamped 
Polygonal Plates vs. Time 
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CHAPTER 5 
BRXPERIMENTAL COMPARISON WITH THEORY FOR A SQUARE PLATE 


5.1, Procedure 

bela expermmentviat Plate, A 34,0 in (4070 an.) 
square aluminium alloy 65S-T6 plate having a uniform 
Cabekness oO: Oyo in, (20.004 in.) wasemounbedtin the 
rolling section of the DRES six-foot diameter shock tube 
agnSNOWnelineties-5.1. Young's molulus Eiand Poisson's 
ratio v were measured statically for two plate material 
Samplesmasy 10.05) x 10° pst and 0.31, and 10,10 x 10° psi 
and 0,306, respectively. Since the small differences 
between these values and the Alcan specified values for E 
and v of 10.0 x 10° psi and 0,33, respectively, theoreti- 
cally would have little effect on the plate response, it 
was decided to use the manufacturer's specified values 
throughout the theoretical analysis. Strain rate effects 
for aluminium, as pointed out by Kolsky [79], are not 
Significant for elastic vibrations and consequently were 
ignored in the theoretical analysis. 

The plate was held in position by one inch 
thick steel clamping plates and rested on a wooden block, 
These clamping plates partially restricted plate edge 


rotation while permitting in-plane motion subject to 
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6' dia. shock tube - 


8" x 4" T beams 1/2'' thick cast iron 





34 Ot x 540 0. 75.71 
aluminium (65ST6) plate 






1" thick front plate 


channel 
weld 
3/4"' aluminium spacer 
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1'' steel clamp 


aluminium plate 
1/2" steel 


3/4" al. facing 


1'' steel front plate 





8" x 4" T beam 


1'"' steel clamps - 3/4"' aluminium spacer 


1" steel front plate 


channel brace 


aluminium plate 
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Back View 


FIG. 5.1 - Plate Mounting Assembly 
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resistance by frictional forces, The edge clamping 
pressure was made as uniform as possible by shimming 
and by adjusting the tightness of the bolts passing 
through the clamping plates and the 3/4 in. aluminium 
spacer plates, By means of this arrangement the plate 
was elastically clamped along its boundary with fair 
uniformity and had an exposed surface area 30.0 in, 
Gad. ein. Jasouare, 

As demonstrated by Lassiter and Hess [1], 
initial curvature in plates increases the complexity 
of analysing deformation by introducing middle surface 
membrane strains even for small deflections, The 
experimental plate initially had a maximum deviation 
promecnestlat,ot 0,013 in. scvergan 3 inch span or 4 
minimum radius of curvature R of approximately 2500 in. 
An estimate of the effect of this initial curvature on 
the transverse vibration frequencies can be obtained 
from the equations for the first two frequencies OF; 
transverse vibration of. an infinitely long barrel 
vault simply supported along its lateral edges which 
are a distance ea | apart. These equations are given by 


Hubbard and Houbolt [80] as 
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and 
Toa 1 2a 
w? =).16|— — [1 - — ee a 
2a,!/ m De ae 
al, a 
Ss | : Cone) 
ia ala 
For the experimental plate 
pene a OO 
and 
2a,/n = 39,63. 
Hence 
eal. oD 
ps ea | See 15,.01.4 | 
: 2a,! m 
and 
an fp 
w? = 16 |——} —[1.000]. 
2 2apike m 








It appears from these values of a and W, that the 
maximum deviation from the flat plate fundamental 
frequency due to initial curvature effects is less than 
one per cent and that tiie deviation in higher frequencies 
is even smaller. These calculations indicate that the 
deviations of the middle surface of the experimental 
plate from a flat surface were not sufficiently large to 
introduce significant membrane strains. 

The effect of transverse shear and rotatory 
inertia on the plate's linear ees could be neglected 


according to Reismann and Greene [20] since the plate 


Lit 












dag 


m 1,35 


oe 
BS | Z 
Be a 4 =} ——— +, 
‘i |} & "rol 
stelq Latneniveqxe sit 16% ri nV OmgaEE 
$f£0,0.= H\ 2S jis aa t 
brs 
sEB,SE = A\, eS “1 bat 
[0,09 ~ = ef | ‘¢. ae | 
“es! i 
nyt 
Gg 2 ale W J a 
Aes of = *y : 
fr ig 8S 8 mo Ase 
oe, ay 


git ted3 “” brs lo ceulav eeoddy mott a a 


isinemsinw? siete vali odd motl nottatveb 

“eis ezol af eatostte syydsvivo {sisint ot eyb yee 
e2lonsupst? isrigid al selssiven ad? ted bas ‘oon y 10 
ee 
bisa dant adeotbak eroteatvelgo seat Bibehrsrs = * at 


oe 6 Am. 


ar sitiioece wet einen vit ai a 
ee ae — salt 


x ; =. 


o. ae = 









~~ = ‘ cs 


in 
or 


del 


thinness ratio h/a of 0.05 was less than 0,1, 

One hundred strain gauges were fixed to the 
plate using Eastman 910 cement. These gauges were 
placed at twenty unique r,6@ locations in each of two 
characteristic plate sectors along the radial lines 
Gee e0 el >a, 30, 7 and 45°, measuringastrain in the 
radial anducircumierenvial directions as shown in Fig. 
B.c. Gauges, were placed on bothesidesr of thejiplate at 
five of these locations in each sector. Budd strain 
Pauses ityoe Cle-levaA, having gauge factors of 2,07 and 
2.08 (+ 0.5%), a gauge resistance of 120 ohms and a 
Baucee enegunl Om, ic 5ein.,, “were used ati most of the 
positions, Micro-measurement strain gauges, type 
EA-06-250-BF-30 having a gauge factor of 2.105 (+ 0.5%), 
a gauge resistance of 350 ohms and a gauge length of 
UPe vUesn., weresused at the positions nearess the 
boundary on the four radial lines where strains were 
relavl Ve lLyasmaLler, 

5.1.2 Static Response. Static pressures 
ranging from 0 to 20 psi’ were applied uniformly over 
the front surface of the experimental plate by bolting 
a steel pressure plate to the one inch steel front 


phate Showneinerie. 5.1, sealing all joints. and 


ST Sp INE LE LTTE ES ET EET RT TI TE 2 NO, IR FE BRE TO I ET IE CIR TE TE RE TCE 


1 All pressures refer to pressures relative to 


atmospheric pressure (i.e., psig) unless otherwise 


stated, 
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8"X 4" STEEL I BEAM 


I" STEEL CLAMPING PLATE 
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DISPLACE MENT 


(NOT TO SCALE) 
a Ne 
5.6" TOP 


F }~>——_ 32" —_4 


ACCELEROMETER 

LINEAR POTENTIOMETER 
PRESSURE GAUGE “CENTRAL DISPLACEMENT 
DISPLACEMENT GAUGE WOODEN MOUNT 


DISPLACEMENT GAUGE FRONT AND BACK 
RADIAL AND CIRCUMFERENTIAL STRAIN GAUGES 
RADIAL AND CIRCUMFERENTIAL STRAIN GAUGES FRONT AND BACK 
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FIG. 5.2 Plate Instrumentation 
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connecting a pressure line from high pressure air 
evyvlinders to the resulting pressure vessel, In this 
WaVemun form pressures were applied only over the 39 x 30 
in, plate surface and the mount was supported against 
displacement at a distance approximately three inches 
from the plate edges. Pressures were measured with a 
Wallace and Tiernan bourdon tube gauge, model FA223, 
having a range of 0 to 50 psi and an accuracy of 

TOE OSes, 

Deflections were measured over the plate 
profile and around the boundary using Starret model 
656-617 dial gauges, each having a 0.0001 in. graduation 
and a 0.400 in. range. These gauges were mounted as 
shown in Fig. 5.3. Strains were measured with a Budd 
P-350 strain indicator and two Budd SB-1l switch and 
balance units. Strains and deflections were recorded 
ateoe Seow ncrements, from 0 «tor720 psi for both the toad 


and the unload portions of the loading cycle, 


S Leo) [ ransLent Response, The cross 


section of the shock tube shown in Fig. 5.1 was 
subjected to air blast waves having normally reflected 
peak overpressures qo ranging from approximately 5 to 
30 psi. The blast waves were formed by the detonation 
of charges of 60/40 RDX/TNT positioned at the central 
axis of a 16 inch bore Naval gun barrel which served 


as the compression chamber for the shock tube. The 
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waves travelled through a 15 degree conical expansion 
chamber into a 6 ft, diameter 1/2 in. thick cylindrical 
tube, The plate test section was located approximately 
160 ft. downstream from the compression chamber, A 
more detailed explanation of the characteristics of 
this shock tube is given by Campbell, Jones and 

Watson [81]. 

This test section was located in the rolling 
portion of the shock tube, Impingement of the blast 
waves caused the section to move very slowly along 
rail tracks. Noticeable motion of the rolling test 
section started approximately 12 to 20 msec, (sec. x 1073) 
after the blast wave arrived at the plate, It was 
considered that these rigid body accelerations pro- 
duced negligible forces on the experimental plate in 
comparison to the forces of the blast loads, 

The blast waves at the test section were 
characterized by a shock front rising to a peak 
magnitude in a few nanoseconds (sec. x 1079) and 
subsequently decaying exponentially to zero magnitude 
in 20 to 90 msec, depending upon the charge weight. 

The positive overpressure phase, having an initial 
normally reflected pressure at the plate Qo and a 
duration ee was followed by a smaller amplitude 
negative overpressure phase, This characteristic shape 


of an air blast shock wave is the result of a flat 
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compression wave pulse being reduced in magnitude 
approximately Srnonene rel iT oy an ensuing rarefaction 
wave, The ratio of the initial reflected blast over- 
pressure an to which the plate actually was subjected, 
LO the initial incident or free fleld overpressure, 
increases with incident overpressure from a value of 
2.0 for normal reflection, A succession of smaller 
amplitude shock waves followed the initial shock wave, 
These secondary waves were caused by longitudinal 
reflections of the original shock wave within the 
shock tube. The period of interest for the plate 
response, including the oscillation of maximum 
response plus several smaller amplitude oscillations, 
always occurred well within the positive phase Gs of 
the initial shock wave. 

Four Ampex 14-track instrumentation magnetic 
tape recorders, having a total of 56 data channels, 
were used in this experiment. The recorders use 
frequency modulation electronics with d.c. response 
range of 0 to 20,000 cps and amplitude accuracy and 
linearity of + 1 per cent. Signals were recorded at 
a tape speed of 60 inches per second and were reproduced 
at 1 7/8 inches per second thereby expanding the time 
scale and reducing the apparent recorded frequencies 
by a factor of 32. Visual reproduction was obtained 


using a Honeywell model 1508 moving light paper chart 
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oscilloscope with a variable chart speed giving a time 
scale ranging from 0,125 msec,/in. to 0,26 msec./in, 
and galvanometers with a response range of 0 to 

2,000 cps. 

Some paper chart oscilloscope records were 
reproduced directly on to graphs while other records 
were digitized with random time-spacing on to computer 
tape using a Telecordex chart reader, The data from 
tne=chart reader were read into the DRES IBM’ 1130 
computer and then plotted with a Calcomp model 565 
X-Y digital plotter. Radial and circumferential 
digitized strain data from the front and back surfaces 
of the plate were interpolated linearly in a computer 
program to yield data with equal time-spacing. Two 
hundred equitime-spaced points were obtained from each 
BULaa Mere cOrumrOreciewtinse 1.5 msec, of plate response. 
The program then correlated strain data from the 
different gauges to yield membrane and bending strains 
and sectional resultants at the 20 locations on the 
plate surface, It was estimated that plate response 
data associated with frequencies of up to 6,667 cps 
would be reproduced using this technique, 

Pressure variations at the test section over 
the central 30 in. square were measured for a series 
of shots using five shear tube piezoelectric pressure 


transducers of DRES design [82] which were shock mounted 
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in a dummy plate. These gauges were calibrated after 
each series of shots to determine gauge signal-to- 
pressure factors. The gauges had a response capability 
of up to 10,000 cps with pressure magnitudes accurate 
TOe2) >) per cent. Each pressure signal was amplified 
and recorded on magnetic tape. Three seconds before 

a shot, a known electrical charge was applied to the 
input of the pressure gauge amplifier and after 
amplification was recorded on each tape channel. In 
conjunction with the gauge factors,this procedure 
provided known pressure calibration signals on each 
pressure channel. It was found that time and spatial 
variations in the pressure over the 30 in. square were 
within the possible gauge errors. Consequently, the 
transient load was assumed to be uniformly distributed 
over the experimental plate surface, 

During the experimental plate tests, 
reflected overpressures were recorded by four shear 
tube gauges shock mounted in the one inch steel frone 
Dleaveras snowteinetie, 5,2. As a further check wprecsures 
were computed from the shock front velocities (as 
measured by side-on pressure gauges mounted in the 
tube wall upstream from the test section) used in con- 
junction with atmospheric temperature, pressure and 
relative humidity measurements taken immediately 


upstream from the plate, 
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Fifty strain gauge signals were amplified 
and recorded on individual tape channels for each shot 
involving strain measurements, The plate's response 
Symmetry was checked for some shots by recording strain 
gauge responses at the same 'characteristic segment' 
radial and circumferential positions in the plate 
Seevorsedesvenatredel@and e2ainwiganss2) weStrain gauge 
responses were recorded at all twenty positions in 
either sector 1 or sector 2 for other shots. Known 
calibration voltage signals, which took into account 
gauge and line resistance of the four-armed strain 
bridge circuit used, were produced with resistance- 
capacitance shunts of DRES design [83]. The calib- 
ration signals were amplified and recorded on each 
tape.channel three seconds before a shot, 

The central displacement was measured with 
a Bourns model 156 Aligno-pot linear potentiometer 
Rhavingvean accnracyeoctea: On/Stperrcent« forse comieinch 
travel and a maximum range of 11/2 inches. This 
potentiometer was mounted on the wooden box beam shown 
injFig. 5.2, which was bolted to the back clamping 
plates so that displacement was measured relative to 
the edge. The travelling shaft of the potentiometer 
was attached to a small circular disk which was glued 
to the plate center, Edge displacements were measured 


at the positions shown in Fig. 5.2 using six Bourns 
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model 141 Linipot linear potentiometers with a range 
of 7/16 inch. These potentiometers were mounted on 

6 in, long 1 1/2 x 1 1/2 in. steel angles which were 
welded to the 8 x 4 in, I beams and their travelling 
shafts were held against the clamping plates by 
compression springs. Linear potentiometer signals 
were recorded on separate tape channels without 
amplification. Signals were calibrated over the range 
of interest by displacing the travelling shafts known 
amounts using feeler gauges, 

Accelerations were measured at the positions 
shown in Fig. 5.2 with two Endevco model 221D piezo- 
electric accelerometers having a frequency response 
rangesOtslromee touw/,J008cps. Accelerometer signals 
were amplified and recorded on magnetic tape. Cali- 
bration signals were obtained in the same way as for 
the pressure transducers, 

A portion of the instrumentation bunker is 
shown in Fig. 5.4. <A complete instrumentation block 


diagrametsesnownInebipy 5.5. 


5.2 Theoretical and Experimental Results 
5.2.1 Theoretical Elastically Clamped Square 


Plate. Eigenvalues, modal parti- 
cipation functions and static central deflections are 


given in Table 5.1 for the first five symmetric modes 
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of vibration of a uniformly loaded square plate having 
uniform elastic edge clamping ranging from the simply 
Supported to the clamped condition, These values were 
obtained using seven collocation points on the boundary 
of the characteristic segment. Values for the 
coefficients do eeeP oe required for, modal participation 
functions,were obtained using 46 equal increments of 6 
for the numerical integration according to Simpson's 
rule, 

Theoretical symmetric frequencies are given 
in Table 5.2 for the experimental plate subject to 
various degrees of elastic clamping. These results show 
that modal frequencies increase relatively more rapidly 


with values of B for decreasing mode number, 


5.2.2 Transient Pressures, Peak reflected 
overpressures 18 and associated durations of the 
positive overpressure phase oe aregpgiven inelabie 5.35 
for 24 separate blast loads on the experimental plate 
resulting from the detonation of charges of RDX/TNT 
Weighing from Oo to; ,0e@lios, Average values of Io 
and t,, obtained from the four shear tube pressure 
records for each shot, were quite consistent for each 
charge weight as were values of q,g obtained from the. 


shock front velocity measurements using [84] 
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Table 5.3 ay 


Plate Transient Pressures 





q, (psi) q. (msec) t | (msec) 
Average of Calculated 
4, Measured from Shock 
Values Front 

Velocity 

4.8 21 

8.6 33 

li2e9) Bi 

20.9 73 


29.8 9 
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ReM* - M2? + 2 ~— 2k 





ee ys ok 
eee “5 Tar ER alee eer cone Coi3)) 
o atm KM? - M? + 2 K +1 
where 
MeV 
RES (5.4) 
and 


ae =fkRrT (5.5) 


mg Eqoy #5. 3) through (545) “k= 8 No and es and ste are 
the specific heats of air at constant pressure and volume, 
respectively, R is the gas constant, T is the absolute 
temperature, erie is the velocity of sound and en is the 
atmospheric pressure, all measured immediately upstream 
BEOmeune Pilave sill the ‘undisturbed ‘air ,” V isithe velocity or 
the incident shock front and M is the shock front Mach 
number 

The experimental values of th measured by the 
pressure transducers were consistently larger by approximately 
1 psi than the values calculated using measured shock front 
velocities, Initial values of q predicted using Eq. (5.3) 
tended to be less reliable at lower Mach numbers due to the 
lower accuracy of shock front velocity measurements. They 
usually did not represent the peak overpressures, in any 
case, since boundary layer effects introduced a small 


secondary rise in the reflected compression wave overpressure 


before it began decaying approximately exponentially, 
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Values of positive phase duration wae 
measured from the pressure transducer records were 
quite consistent for different shots having the same 
charge weight. Variations between the values of qQ, 
obtained from either pressure transducer records or 
shock front velocity measurements for different shots 
having the same charge weight were usually smaller than 
variations between the values of aS obtained for each 
shot using the two different techniques. Consequently, 
it was considered reasonable to use one value for a5 
and Le for each charge weight in the theoretical plate 
transient response calculations, with Io being chosen 
nearer to the values measured by the pressure trans- 
ducers. It is seen from Table 5.3 that the possible 
errors in theoretical values chosen for qa, range from 
approximately + 19 per cent to — 8 per cent for a-0.5 
POeSnOlsCOss OL per centbsto — J per cent for ali0.001b. 
shot. 

In view of the errors possible ingq,, it 
seemed reasonable to assume a transient load history of 
Che a0Orm 


q(t) = qo (1 = t/t) exp(— t/t), tis toe Ge) 


for the theoretical response calculations rather than 
to numerically integrate the Duhamel integral, kq. 
(3.33a), for individual pressure records. Typical 


experimental and theoretical transient blast pressures 
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TOrec@l oe bweanOurare sshown in Fig, 5.60, “At ileast 
some of the higher frequencies present in this 
experimental pressure-time record resulted from gauge 
"'ringing’. For this shot, an initial reflected over- 
pressure as of approximately 8.0 psi would have been 
obtained from accurate shock front velocity measure- 
ments. However, a value of 8.6 psi was used for the 
theoretical initial pressure do Ine orders lo sopsain a 
more reasonable pressure impulse for the period of 


interest which was approximately 15 msec. 


5.2.3 Damping. Damping introduces non- 
conservative forces which dissipate energy in a system 
undergoing cyclic stressing, usually by converting 
mechanical energy to heat. It was considered that both 
structural and air damping would affect the transient 
response of the experimental plate, 

Structural damping is composed of material 
damping, also termed internal, hysteritic or visco- 
elastic damping, and interface or Coulomb sliding 
damping, both of which are described in detail by Lazan 
and Goodman [85]. Generally, structural damping is 
nonlinear and is not representable by a viscous damping 
model. 

Air damping, as discussed by Hubbard and 
Houbolt [80], is composed of a component in-phase with 


acceleration producing an apparent "added mass" and a 
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component in-phase with velocity termed "radiation 
resistance", According to Hubbard and Houbolt, the 
acceleration component of air damping is negligible 
for the experimental plate since 2a)/h = 39.6 < 500. 
However, particularly for blast waves, the velocity 
component of air damping may not be negligible. 

The method of calculating the radiation 
resistance 1s dependent upon the value of Wan / oh es 
where eye. Lse=the-velocity sof sound) inthe air 
surrounding the plate. The values of ee ranged 
from approximately 1100 to 1300 fps for the air blast 
overpressures associated with the plate experiments, 
Radial strain vs. time records at p = 0.377, 9 = 15° 
forei a5 and 6.0 lb, shots, shown in Figs, 5.7 and 
5.8, respectively, indicate that the fundamental 
frequency of vibration of the experimental plate was 
approximately 165 cps for the 1.5 1b. shot and 160 cps 
for the 6,0 lb. shot. Relating these frequencies to 
the theoretical fundamental frequencies for different 
edge clamping values given in Table 5.2, shows that the 
experimental edge clamping approached the simply 
supported condition, 


Thus for the experimental plate 
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FIG. 5.7 - Plate Experimental Radial Strain at P = 0.377, @= 15°, sector 1 
and Pressure vs. Time for a 1.5 1b. Shot 
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FIG. 5.8 Plate Experimental Radial Strain at p= 0.377, @ = 15°, sector 1 
and Pressure vs. Time for a 6.0 lb. Shot 
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According to Hubbard and Houbolt, for a 


simply supported square plate having a Poisson's ratio 
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air 
inekasea (54 7addand (547b) ica/ued ea dsathearatiovoft 
viscous damping of the fundamental mode to critical 
damping computed for air on one side of the plate 


Only, Y andy are the densities of air and plate 
Bair 
material, respectively, and c = E/y is the velocity 


of sound in the plate material. 
At atmospheric pressure and 70°F., (c ) 


ar atm 


2 =6 
a t pmmere) d Ovl2 Xxqad 
S approximately 11 fps an OreT une is 


1b.sec.* ai For the aluminium plate c is 16,600 fps 
and y is 0.25 x 1073 Woe ona For a normally 


reflected air shock wave, (y ) and (c ) are 
a1 rece l Sir) pet 


given by [84] as 


Cy ) 
air ref G7M> 242))3M2 
= (5.8) 
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(c ) (cam? - 1)¢M? + 2) 
air ref a 


ee ) 3M 
air atm 


° (5.9) 


Theoretical values of plate viscous damping 
due to air coupling are given in Table 5.4 for a range 
of blast loads. Assuming one side of the plate was 
exposed to atmospheric pressure (q = 0 psi) while the 
other was subjected to air blast loads, theoretical 
initial values of alr damping of the fundamental mode, 
obtained by adding the damping on each surface given by 


Pomme av ,eraigeaerrom 0,50 per cent of critical for a 


blast load having qo 4.8 psi to 0.67 per cent for a 
blast load having sar 29.8 psi. Even for a blast load 
having q, = 463 psi, initial air damping would have been 
Oye eperrcent Oo. critical, If Equet5s.7b) had 
governed, then air damping would have been slightly 
larger for all blast overpressures and it would EME 
increased more rapidly with blast overpressure, However, 
even assuming that Eq. (5.7b) governed, air damping never 
exceeded 1,6 per cent of critical for a blast load with 
ad, as large as 29, 8Pps ik 

Experimental values of viscous damping computed 
from the strain transients in Figs. 5.7 and 5.8 using 


c/w it m 
lh Fens ing (5.10) 
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Table 5,4 


Plate Theoretical Air Damping On One Surface 


) 
) 


( 


Cc ) 
air ref 


Ke ) 
air atm 
Bq. €58/a)) shq.05. 1b? 


0.0023 0.0044 


(y 


air’ref 


(y 


air’atm 


0:50 0.27, 0.0054 
0.0031 0.0067 
040033 0.0078 
0.0039 0.0098 
0,0044 0,012 


0,011 0,080 





areseiven in Table 5,5, In Eq. (5.10) Xn is the 
amplitude of oscillation m for free vibration. Values 
of ne used in the viscous damping calculations were 
obtained from the forced response strain transients by 
measuring the strain amplitudes between successive 
minimum and maximum strain values of the fundamental. 
mode, where the maximum value is labelled oscillation m, 
In this manner the effect of the time decaying pressure 
transient on vibration amplitude was partially removed 


and values of the ratios xf %, for the different strain 
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amplitudes were thought to approximate closely the 


related free vibration values, 


Table 5,5 


Plate Experimental Viscous Damping 


Peak m to Peak n | 1,2 sly 1,4 15 2,4 335 
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Damping measurements for the 1.5 lb. shot 
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commencing at the maximum strain levels occurring for 
oscillations 1, 8, and 11 which were associated with 
shock waves arriving at’ 0, 42 and 111 msec., indicate 
that viscous damping decayed with the stress amplitude 
but was considerably larger than the values estimated 
for air damping at all blast overpressures, Comparison 
of the damping values for the 1.5 and 6.0 lb. shots 
when dg was 8.6 and 2059 psi, respectively, reveals the 
nearly linear increase of initial values of viscous 


damping with stress amplitude, The rapid decay with 
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time of viscous damping for both shots indicates that 
initial values of damping could have been considerably 
larger than the values measured from the first two 
oscillation peaks. This linear increase in viscous 
damping with stress amplitude did not. occur for the 
10.0 lb. shot when plate damping levels were nearly 
equivalent to those for the 6.0 1b, shot. . It was 
thought that this was a result of a decrease in the 
plate edge clamping friction, This is discussed 
furthemrin tection S,i2.5. 

The strain response magnitudes of the higher 
modes, as shown in Fig. 5.7, appear to have remained 
at least as significant relative to the fundamental mode 
magnitude after several oscillations of the fundamental 
mode as they were for the first oscillation. This 
indicates that the higher modes were subjected to 
considerably less viscous damping than the fundamental 
mode, because of the smaller amplitudes associated with 
the higher modes. 

It was considered from these results that 
most of the experimental damping resulted from the 
interface damping which occurred at the plate edges in 
the mounting assembly. Since the experimental values 
of viscous damping were so large, it is possible that 
these values represent not only damping but also the 


effect of the presence of in-plane tensile strains which 
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are discussed in Section 5.2.5, However, Lassiter and 
Hess [1] found that combined structural and radiation 
damping increased rapidly with stress at higher 
maximum stress levels for a rectangular panel subject 
to random acoustic input and obtained similarly large 
values of damping at similar stress levels, 

Coupling between the symmetric transverse 
modes due to damping was not considered to be important, 
even at these high damping levels, since the modal 


frequencies were 'well separated', 


5.2.4 Displacement. The theoretical 
variation of static dimensionless central displacement 
with edge clamping factor 8 is. shown =inelid g..5i9eLonr 
the values given in Tables 5.1. Static central 
deflections measured at 10 and 20 psi show that values 
of Bg decreased slightly with increased load, While 
the value of B a/D was approximately 8 for the static 
loads as indicated by central deflections, it 
decreased to between 0 and 0,523 for the blast loads 
as indicated by values of fundamental frequencies. The 
reduction of edge stiffness for the blast loads resulted 
for two reasons, Firstly, for the dynamic loads the 
mount was supported against displacement by the 
8 x 4 in. I beams at a distance of approximately 11 in. 


from the plate edges, while for the static loads the 
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FIG. 5.9 Theoretical Variation of Static Central Displacement with Clamping 
Factor a for a Square Plate 
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mount was supported only 3 in. from the plate edges. 
Further, the complete 6 ft. diameter section was 
subjected to blast loads while only the 30 in. square 
plate was subjected to static loads. 

Theoretical and experimental central 
Gieplacements vs. vame areyshown; in Figs, 5.40 and 


Date Org... 5s andeal.5 lb. eshots, mrrespectively, Since 


- maximum response occurred during the first cycle of 


the fundamental mode, theoretical and experimental 
results were compared for the first few cycles only. 
Theoretical transient displacements computed using 
Bi a/D = 0.523 and c/w, =0.1 agree well with the 
Gesu OF sbne /0,5 1b. shot, however! for the 1,5 1b. 
shot the experimental clamping was slightly more 
Poen Die sande tne wiscous damping exceeded 10 per cent 
Ofecritacawyor the first cycle, 

Dimensionless central displacements are 
Snow meieties., 5.12 for.tatic and blast Loads.) Since 
ratios of maximum dynamic-to-static displacement 
exceeded the maximum dynamic load factor of two for a 
flat pulse, it. ids evident that the elastic edge 
clamping was more flexible for the blast loads, Non- 
linear effects were almost unnoticeable for static 
central displacements ranging up to one-tenth the plate 
thickness although there appeared to be slight non- 


linear softening, likely resulting from an edge 
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1D to 215) EXPERIMENT 
m “ae i 
Wiiees) audience ic THEORY B'a/D = 0.523, £./w,=0 
q.a 
5 — 0.1130 a | ss ome Ba/D= 0.523, Ci7 w@; =O! 
FIG. 5.10 


Plate Theoretical and Experimental Central Displacement vs. 
Time for a 0.5 1b. Shot 
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FIG. 5.11 Plate Theoretical and Experimental Central Displacement vs. 
Time for a 1.5 lb. Shot 
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flexural softening. However, there was a general hardening 
Oey ave wespense With Increasing blast. load fori all, values 
of w/h with one exception, This hardening was observed 
previously by Lassiter and Hess [1] for a rectangular plate 
subjected to acoustical loading, Theoretical results given 
by Bauer [86] for a simply supported square plate with a 


stress-free boundary and n/a = 0,06 (h/a, = 0,05 for the 


'. experimental plate) indicate that nonlinear hardening effects 


Morea silat pressure pulse are nov significant until the 
central displacement-—to-thickness ratio exceeds 0,3, 
Consequently, this hardening nonlinearity is probably a 
result of a combination of the increase in damping with 
Weiasct toad and the=sottening in the ratio of edge-to-—central 
transverse displacement discussed below, 

Ratios of maximum transverse edge displacements- 
to-central displacements for the blast loads are shown in 
Pip. 5,13 for the edge positions labelled 3"and 6 injFig, 5.2. 
Displacements measured at the other boundary positions l, e 
end 5,'had values ranging-up to 30™per cent smaller than che 
welues measured at position 3. However, the increase in the 
proportion of edge-to-central displacement with the load was 
observed at all boundaries, This increase was noticed for 
the static loads as well, but the static edge displacements 
were approximately one-tenth the value of the dynamic edge 
displacements as a result of the two. different loading 


arrangements, These results indicate that an elastic edge 
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*NUMBERS IN BRACKETS REPRESENT RATIOS OF MAXIMUM CENTRAL DISPLACEMENT 
TO PLATE THICKNESS whh. 


FIG. 5.12 Plate Experimental Static and Maximum Dynamic Central 
Displacements vs. Load 
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LOCATIONS OF EDGE DEFLECTION MEASUREMENTS ARE SHOWN IN FIG. 5-2 


FIG. 5.13 Plate Experimental Ratio of Maximum Dynamic Edge Displacement — 
to-Central Displacement vs. Load 
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condition relating transverse displacement to Kirchoff's 
effective transverse shear force should also be considered 
theoretically for the dynamic case, 

Integration of accelerometer records corroborated 
the values of edge motion obtained from the linear 
potentiometer records and indicated that the wooden mount 
for the potentiometer which measured central displacement 
‘followed the edge displacement motion fairly rigidly. 

Fundamental frequencies measured at the plate center 
end, thes boundary positions labelled: 1 and’ 5 in Fig, 5,2 are 
given in Table 5.6. The fundamental frequency at the plate 
center for the first peak-to-peak cycle was consistently 
lower than the other measured frequencies, A freqency 
reduction to 161 cps could be explained by a reduction in 
the elastic clamping constant B, so that the boundary 
approached the simply supported condition, Large values 
of viscous damping could reduce the frequency further, 
Bowever wo reduction in the, fundamental, frequency, from 
161 to 147 cps, for example, requires a linear viscous 
damping of 41 per cent of critical, which is considerably 
larger than measured values. Consequently, this low 
initial plate frequency probably resulted from the. fact 
that initially the edge displacement did not vanish completely, 
as noted previously, and also that the edge conditions were 


time-dependent, 
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Table 5.6 


Plate Experimental Fundamental Frequencies 


Frequency w fen (cps + 4 cps) 


aka 1) 143 143 1136 134 
167 













Plate Center 











Edge 1 





167 














Edge 5 169 169 








Plate center 


Bagel 


Cycle peak-to-peak 


Edge 5 







ogo) L259 20792900 


q (ps) 
O° 






5.2.5 Strain, Typical experimental bending 
and membrane strains vs, time are shown in Fig. 5.14 for 
ail.5 1b. shot. The frequency of the first peak-to-peak 
cycle of the fundamental mode was approximately 143 cps 
for the bending strain while for the membrane strain it 
was 166 cps. Comparison of these frequencies with the 
frequencies given in Table 5.6 for q, = 8.6 psi indicates 
that the membrane strain frequency had the same value as 
the plate edge transverse motion frequency. Similar 
correlations between membrane and transverse edge motion 


frequencies occurred for the other blast loads. 
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q,= 86 psi seececcceccerceee EXPERIMENT, SECTOR | 
to = 33 msec ————_ EXPERIMENT, SECTOR 2 
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FIG. 5.14 Plate Experimental Bending and Membrane Strains vs. Time at 
p = 0.047, 6 = 0° for a 1.5 lb. Shot 
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Typical values of maximum bending and membrane 
strains recorded at the ten positions on the plate surface 
having strain gauges front and back which are shown in 
Picmore.manremclvenein Taple 5,{/ for a static load of 10 
Poimancu@ite able 5,O 1O0r a blast load Of O,O ps. Initial 
reflected overpressure, These results show that although 

membrane strains were quite small in proportion’to asso- 
piaved bending strains for static loads, they had signi- 
ficant tensile values for blast loads even when the max- 
imum center displacement-to-plate thickness ratio was 
Lesecmcran 0,1, 


Table 5.7 


Plate Experimental Bending and Membrane Strains 


forcunrormMeotavicrhboadwor lO psi 


POSLUTGH EOC Sector 
(in.-deg.) 
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Table 5.8 


Plate Experimental Maximum Bending and Membrane 
Strains for Uniform Blast Load, quns 8.6 psi 


and oe = 0,033 sec, 


Positioner. Oo Sector 
(in.-deg. 3 





The large Teniarion in the values of membrane 
strains recorded at the same r,9 positions in sectors l 
and 2, in comparison with the reasonably symmetric values 
for bending strains, and the correspondence of the values 
of the plate membrane and transverse edge frequencies, 
indicate that the membrane strains principally resulted 
from a nonsymmetric boundary Serepehing motion associated 


with the boundary transverse motion, In contrast to this 
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conclusion, Eikrem and Doige [87], who also observed membrane 
strains for the small deflection response of a simply sup- 
ported square plate to blast loads in a recent test conducted 
at DRES, seemingly attributed these membrane strains to some 
coupling effect with the transverse motion of the plate itself. 

Results given in Section 5.2.4 which showed that 
transverse edge motion was not symmetric and was approximately 
een times larger for blast loads than comparable static loads, 
further corroborate the postulate that the membrane strains 
result from in-plane boundary motion, Using this postulate 
these results explain both the increased magnitude of the 
membrane strains for the blast loads in comparison to the 
static loads and the nonsymmetry of the membrane strains, 
Since a simultaneous in-plane outward displacement of only 
0.00075 in. at each of two opposite boundaries was required 
to produce a tensile membrane strain of 50 ue, it is quite 
conceivable that this plate stretching could result from 
boundary rotation in conjunction with an absence of free 
slippage of the experimental plate between the edge clamping 
plates, 

PlovsmMinerig. >, 1S tof ithe radial tandi-circumperential 
bending strains near the plate corner vs, load reveal the 
general nonlinear softening of the edge support resistance 
to rotation normal to the boundaries with increasing blast 


loads. 
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Variation of membrane strains recorded near 
the plate corners with blast loading is shown in Fig. 5.16. 
Generally, the membrane strains reduced nonlinearily 
with increasing blast load. However, the in-plane 
motion of the edge clamping plates probably increased 
nonlinearly with the load, in conjunction with edge 
transverse motion. Therefore, the boundary clamping 
friction must have decreasedwith increasing load as was 
eoncluded#in isectiton_ 5.24 3 

The present theoretical solution can pre- 
dict correct flexural mode shapes for the experimental 
plate only if the membrane and bending stress actions were 
completely uncoupled, Theoretical results given by 
Dawe [47] for the response of rectangular plates show 
that the degree of coupling between the membrane and 
bending stress actions depends not only upon the magnitude 
of the membrane stresses and their distribution but also 
upon the plate edge conditions. Even if there was 
complete uncoupling for the experimental plate response, 
the presence of the membrane stresses would change 
the modal frequencies of transverse vibration. Results 
given by Dawe for v = 0.3 show that transverse vibration 
frequencies decrease with the ratio of the membrane 
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FIG. 5.15 Plate Experimental Maximum Dynamic Bending Strains vs. Load 
at p = 0.778, 0 = 45° 
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FIG. 5.16 Plate Experimental Maximum Dynamic Membrane Strains vs. Load 
at p = 0.778, @ = 45° 
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buckling stress factor K.,, where K is positive for 
compressive membrane stress o, A modal frequency 
approaches zero when K approaches ee forthe mode, 
as was pointed out for shallow shells in Section 1,3.2. 
Dawe indicates that mode shapes remain essentially 
unchanged if the modal frequencies vary linearly 
Dxith (Kes 
| According to Dawe, the mode shapes of simply 
Supported square plates subjected to uniform uniaxial 
stress remain unchanged for all values of K < Kae. 
and the modal frequencies vary as fl - K/Ko, -deForethe 
fundamental mode of a simply supported square plate 
eee = 4, Consequently, for a uniform uniaxial strain 
of 50 pe in the experimental plate, K = 0.085 and the 
theoretical fundamental frequency increases by 1,05 
per cent from the value given by the present theory for 
simply supported edges in which the effect of membrane 
Stresses 16 not kKtonsidereds 

Dawe showed that the mode shapes of a clamped 
square plate subjected to uniform biaxial stresses are 
affected only slightly by these stresses for |K| < 5. 
His results indicate that the fundamental flexural 
frequency of the experimental plate with clamped edges 
would increase by approximately 0.8 per cent from the 


value given by the present theory for a uniform biaxial 


strain of 50 we. 
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From these results it appears that the 
flexural response of the experimental plate to blast 
loads can be predicted reasonably accurately using 


the present theory. 


5.2.6 Stress Couples. Experimental 


values of static and dynamic flexural stress couples 
recorded near the plate center vs. load are shown in 
Fig. 5.17. These stress couples were computed from 
the 'true' bending strains with membrane strains 
removed. Although there appears to be some plate 
elastic hardening for incident blast loads between 
4,8 and 8.6 psi, most nonlinearity is within the 
pressure measurement errors, Consequently, the 
effect of the nonlinear decay of membrane strain and 
the elastic softening of the edge clamping with blast 
load must have been offset by increased viscous 
damping. 

Dimensionless static experimental and 
maximum dynamic Be eirortal and theoretical flexural 
stress couples formulated in Appendix A.1l are plotted 
along the radial lines 6 = 0°, 15°, 30° and 45° in 
Figs. 5.18 to 5.21, respectively, The experimental 
values were obtained for a 10 psi static load and an 
8.6 psi reflected overpressure blast load. Since all 
dynamic membrane strains meceuvensi en dynamic stress 


couples calculated from 'true' bending strains were 
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Static and Maximum Dynamic Plate Experimental Flexural Stress 


FIG 5.17 


Couples vs. Load at p = 0.047, 9 = 0° 
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always smaller than stress couples calculated from 
"apparent' bending strains which were measured on the 
back of the plate remote from the blast, The uneven 
variation in the experimental maximum dynamic stress 
couples with radius appears to result entirely from 
the é€limination of membrane strains from the stress 
cOupie calculations at some positions but not av others. 
The smooth variation in experimental maximum dynamic 
stress couples shown in Fig. 5.19 is likely represen- 
tative of the true variation assuming uniform membrane 
strains. However, the true magnitudes of the dynamic 
stress couples would be approximately 0.5 ordinate 
divisions smaller in magnitude than the stress couples 
given in this figure, 

Theoretical stress couple variations along 
the four radial lines are shown for edge clamping values 
gt a/D = 0.523, 1.57 and 10.46 and for viscous damping 
values 6 ON FAO 0, anu. U,c. LNeOrevlcalearesulus 
calculated using an edge rigidity of Bo a/D = 10,46 shown 
in Figs. 5.18 and 5,21 appear to give the best agreement 
with the static experimental results. Stress couples 
computed for gt a/D = 0.523 and ees = 0,2 shown in 
Figs. 5.19 and 5.20 give the best agreement with the peak 
magnitudes of stress couples recorded for the blast 


loading. 
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Theoretical and experimental flexural stress 
couples) vay Gime are compared at four different positions 
OnmenempLave 1m Hies, 5.22 through 5.25 for. the blast load 
having q, = 8,6 psi and t, = 33 msec. Experimental stress 
couples were calculated from 'true' bending strains in 
Pie ee cok conde eat LDeOreci cal flexuralistress 
couples could not be computed accurately in the region 
O0<p<0.06. However, experimental results in Fig. 5.18 for 
uniform ystGatic sehoads indicate that the stress couples 
Varied #iyt re with radius for p < 0,1. Consequently, 


comparisons in wie, 5.22 of itheoretical stress couples 


computed at 0 0.0707 with experimental stress couples 
computed at p = 0.0471 should be valid for uniform dynamic 
loads, a Membrane strainsawere not eliminated» from the 
experimental stress couple calculations at o.= SO 

6 = 45° in Fig. 5.24. This gave larger magnitude experi- 
mental stress couples which agree quite closely with the 
theoretical values, Theoretical and experimental results 
agree reasonably well at the other positions, however for 
this blast load initial experimental viscous damping 
appeared to be greater than 10 per cent of critical and 


the edge clamping value of Bra /0D was less ¥thane0, 523. 
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CIRCUMFERENTIAL STRESS COUPLE 


a 
Mo ony 





RADIAL STRESS COUPLE 


a 
M. 9 "D 





— 73-3.379 —— EXPERIMENT 
op 0S epee THEORY Bia/D=0.523,6; / w;=0 
a = 0.2024 Se B,o/D= 0.523, 0; / w; =01l 


FIG. 5.22 Plate Experimental and Theoretical Flexural Stress Couples vs. 
Time ate = 0.0471, @ = 0° andp = 0.0707, 6 = 0°, 
respectively, for a 1.5 1b. Shot 
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CIRCUMFERENTIAL STRESS COUPLE 





RADIAL STRESS COUPLE 


Qa 
M, :) “D 





D Jy = 3.379 —— EXPERIMENT 
m qa i 
AS ee THEORY Ba/D= 0523,6,/ w; =0 
do” = 0.2024 
D ; wa Ba /D=0.523,0;/ =O. 


FIG. 5.23 Plate Experimental and Theoretical Flexural Stress Couples vs. 
Time at p = 0.3771, @ = 0° for a 1.5 1b, Shot 
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CIRCUMFERENTIAL STRESS COUPLE 





RADIAL STRESS COUPLE 





a 
Meg ‘D> 


0 
Dias ——— 
1 44 = 3.379 EXPERIMENT 
m - 
3 


eeeeneeee= THEORY 'a/D = 0.523, 6; /#;=0 


—2— =02024 — — — B'a/D = 0.523, $;/%; =0.1 
D t 


FIG. 5.24 Plate Experimental and Theoretical Flexural Stress Couples vs. 
Time at.p = 0.5657, @ = 45° for a 1.5 lb; Shot 
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CIRCUMFERENTIAL STRESS COUPLE 





dD ~s = 3379 EXPERIMENT 

He -----=----- THEORY g'a./D = 0.523, f,/w,= 0 
3 | 

se = 0.2024 ——- Bia /D =0.523,f/ w,= 01 


FIG. 5.25 Plate Experimental and Theoretical Flexural Stress Couples vs. 
Time at P = 0.7778, @ = 45° for a 1.5 lb. Shot 
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These results indicate that satisfactory 
theoretical predictions for stress couples can be made 
for the entire plate surface, even when edge transverse 
displacement does not vanish completely and some membrane 
strains are present, providing viscous damping coefficients 
and edge spring constants relating the normal slope to 
the tangential flexural stress couple vector can be 


estimated with reasonable accuracy, 
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CHAPTER 6 
SHELL NUMERICAL RESULTS 


6,1 Transverse Vibrations of Thin Shallow Shells with 
Homogeneous Edge Conditions 

Symmetric eigenvalues for shallow spherical 
shells with regular polygonal plan were obtained for 
various values of the geometric parameters a/h and a/R 
and a range of edge conditions which were satisfied 
using the boundary collocation technique. Four edge 
conditions were satisfied at (L-l1) boundary collocation 
points while only the first three edge conditions were 
satisfied at the remaining collocation point which was 
always chosen at the shell corner, Consequently, an 
equal number of simultaneous boundary equations and 
unknown integration constants was obtained with ity 
collocation points on the boundary of the characteristic 
segment when the series solutions for nt and et eG ca 
(Sie aynande (3,270) Le truncated with n ranging from 
0 to (L-1). 

Modal participation functions were computed 
for some of the shells for uniformly distributed and 
central point loads, Numerical solutions were obtained 
using the digital computer program given in Appendix 


C.2 and machine computation using 16 figure accuracy, 
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Some eigenvalues were obtained using only 10 figure 
computational accuracy, 

The sets of edge conditions satisfied and 
the associated terminology, well established only for 
the simply supported and clamped shells, are listed 


below: 


a) Clamped shell: 


c aw(p, 6) . 
w(6,8) = on = Ff (6,9) = € (6,8) ea, 
V VS ss 


b) Simply supported shell: 


w(6,0) = M (6,6) = F (6,6) = € (,6) = 0. 
vs Vs SS 


ec) Clamped shell with sliding free edges: 


_ dw 6,6) 2 
w(6,8) = poe =oFo (6,0) =F (p,6) e200; 
Vv VV 


d) Clamped shell with sliding clamped edges: 


aw(6, 9) 7 
w(p,6) = Sao ari Be DSP iaa cy (6,0)"8= 0. 
Vv VV s 


(6.1) 


(6.2) 


(653) 


(6.4) 
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nly a4, 
e) Simply supported shell with sliding free edges: 
w(5,6) = M (6,6) = F (6,6) = F (6,6) = 0, (6.5) 
f) Simply supported shell with sliding clamped edges: 
w(6,0) = M (6,0) = F (6,0) = © (8,8) = 0. (6.6) 


The terminology applied to edge condition 
Eqs. (6.3) through (6.6) in each case describes 
successively the edge conditions in transverse dis- 
placement and stress function. Thus, for example, a 
shell satisfying edge condition Eq. (6.3) is referred 
to as a 'clamped shell with sliding free edges' where 
"clamped' refers to the conditions w(p, 8) = aw(B,6)/av = 0, 
'sliding' refers to the condition F (6,8) = 0 and 
'free' refers to the condition Fa(6,9) = 0, 

In the numerical analysis, age Bal tallioe 
Eqs. (6.1) and (6.2) were never satisfied at 6 = 0° 
since, if they were, the determinant value of the boundary 
equation matrix would vanish for all k as Be LS -ZEDo 
by derinition (seetbqs, (3..14b);,. (3.17c) and) (3.19c)) 
along the radial line ® = 0°, However, edge condition 
Eqs. (6.3) -and.(6.5) could be satisfied at ® = 0° since 


the fourth edge condition was never satisfied numerically 
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Satisfaction of the edge condition 
e,,(B,8) = 0 in conjunction with the condition 
w(6, 9) = 0 is equivalent, for symmetric shell response, 
to satisfaction of the edge conditions u(6,8) = w(6, 9) = 0, 
where u. represents in-plane displacement tangential to 
the boundary. Satisfaction of the edge condition 
F__(6,9) = 0 in conjunction with the conditions 
ug (6, 6) = au_(6,8)/av = 0 is equivalent, for symmetric 
response, to satisfaction of u(B, 8) = u (6,6) = 0, Here 
UL represents in-plane displacement perpendicular to the 


boundary. Consequently, specifying only 


w(B,6) = F (6,6) = €,.(5,6) = 0 (6.7) 


is not equivalent to specifying 
w(6,0) = u_(6,6) = u_(p,6) = 0 (6.8) 


anaevneseage condition Has. (6.1) and (6.2) are not 
equivalent to the usual conditions specified for clamped 
and simply supported shells. Physically, satisfaction 
of edge condition Eq. (6.7) eliminates edge elongation 
and in-plane shear but permits edge curving while 
satisfaction of edge condition Eq. (6.8) eliminates 

edge elongation and curving but permits edge in-plane 


shear, As far as is known, the edge condition of the 
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function du _/av has never been considered analytically. 
However, it is felt that the effect of satisfying 
either Eq. (6.7) or Eq. (6.8) is almost equivalent 

for symmetric shell response, except possibly for small 
values of p, since values of du (6,8) /av should be 
small. This near-equivalency is confirmed in Sections 
6.1.2 and 6.1.3 for simply supported and clamped shells 
with 12 and 15 sides. 

Oniashvili [88] and Nowacki [33], in 
analysing the vibrations of a shallow spherical segment 
with a circular boundary, used the edge conditions 
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as being the equivalent of 
ow , 
u =u =w=-—e= OQ, (6.10) 
0 6 ap 
However, this equivalency does not appear to be 


3 
correct, even assuming that U8 = 0 at the boundary, 
3p 


since integration of the stress resultant-displacement 


relations at the boundary yields 





oF 
u =c— + f(@) (6.11) 
p ‘ap 
and 
oF af( 6) 
Cie ae Cae ot , (6.12) 
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In these equations f(@) is an arbitrary function of 6 
and Cis Cos and ¢, are constants, A numerical check 
on the correctness of edge condition Eq. (6.9) applied 
to a shell with a polygonal boundary having 12 sides 
as an alternative to edge condition Eq. (6.10) applied 
to a shell with a circular boundary, appeared to give 
unsatisfactory modal frequencies when compared to 
values given by Kalnins [56], 
Weinitschke [89], in analysing the 
asymmetric buckling of a clamped shallow spherical shell 
with a circular boundary, used the edge conditions 
r) ’) OW 
€ =r—[e J]-e -—[{e J=wete—e=0O COR 13) 
60 are G9 Pr 205mg or 
as being the equivalent of edge condition Eq. (6.10). 
Although these two edge conditions in strain are the 
equivalent of satisfying the edge conditions u =u = 0 
for a clamped shell with a circular boundary, pney fee 
not applicable for a simply supported shell since they 
must be satisfied in conjunction with the two edge 


conditions w = = = 0 in order to be valid. 
r 


In the numerical results which follow, varying 
degrees of difficulty were encountered in obtaining 
fundamental eigenvalues for different edge conditions, 


shell geometries and boundary collocation sets. It is 
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postulated that essentially it was the insensitivity of 
AacOmecrnanges ein katorusmall k which made it. difficult to 
find accurately the lowest values of k for which the 
various boundary equation matrices had vanishing 
déverminants values. welhe decreasing sensitivity of A, to 
changes in k as k* becomes small in comparison to 


12(1 - v*)a*/(Rh)? is apparent from the relation 





where 


h = Eh/(maRw?). 

The function A is present in boundary equations 

involving the transverse displacement and its derivatives. 
Inaccuracies in calculating derivatives of Jf ke) and 
ey DOnek>~leueinee ower order derivatives, of these 
functions further amplify the problems caused by this 
PNSCHSLULVLUY «OL A, 

This postulate explains the variation in the 
degree of difficulty encountered in the numerical 
analysis in solving s for the different edge conditions. 
Fundamental eigenvalues associated with shells 
satisfying edge condtion Eqs. (6.1) or (6.2) were 


found relatively easily since the shell edges were 


reasonably well restricted against in-plane normal 
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boundary motion and consequently had values of Kk’ which 
were large in comparison to 12(1 - v*)a*/(Rh)?. Values 
Of; kK associated with shells satisfying edge condition 
Eqs. (6.4) or (6.6) were found for many shell geometries 
since u_ (8,9) was still partially restricted and con- 
sequently values of Kk" usually were not too small in 
comparison to 12(1 - v?)a*/(Rh)?, However, values of kK 
associated with edge condition Eqs. (6.3) or (6.5) were 
not found for most shell geometries since u_ (8,0) was 
unrestricted and consequently values of k" must have 
been small Sir .comparison tod - (1. —'y7) a" /CRh)*, 

This postulate also explains the increasing 
difficulty encountered in obtaining e as “p increased, 
since the magnitude of K generally decreases with 
increasing *p,. ‘at least» for the shells for which on was 
obtained, while 12(1 - v?)a"/(Rh)? is independent of p. 
inivadditiom, ait explains why the value of kK for’ any 
particular shell boundary shape and edge condition was 
solved more accurately for smaller magnitude values of 
the product of the shell geometry parameters = and =. 

Whesagiificulty in solving kK for edge condition 
Eqs. (6.4) and (6.6) for larger values of p might also 
be explained partially by the fact that, as p in- 
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gecreasesand finally vanish for a circular boundary. 
But the edge conditions 


w(6,0) = F (6,6) = e« ($,6) = —8($,6) = 0 (6.16) 
vs ss OV 


imply that in-plane displacement normal to the boundary 
vanishes so that F (5,6) cannot vanish as well. Thus, 
it becomes increasingly difficult to satisfy Lali) = 0 
as the polygonal boundary approaches the circular shape 
Hiees Coe 6 e000. 


ss 


Syeatieaa | Simply Supported Shells with Sliding 
Clamped Edges. The only previous analysis of the 


Oranteverse Vibrations, of a shallow shell of noncircular 
planform, as far as is know, has been given by Reissner 
[50], who solved the natural frequencies of a shallow 
thin paraboloidal cap with a square base. Neglecting 
the effects of longitudinal and rotatory inertia and 


transverse shear and employing the edge conditions 
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Teo. (0, lone trjemel« 31 5... Lor symmetnicomodes and 
H is the shell height or sagitta, 

Although Reissner states that Eq. (6.15) applies 
Momdmapner oalimidale surface, ithe surface to which 1¢ 
actually applies is slightly paraboloidal since the quantity 
(H/a)* has been neglected with respect to unity in the 
curvature calculation. For a truly spherical middle 
surface, Eq. (6,15) becomes 


W SSS Sr 


Rect] TO =e ets i 
2 2 SSS oe serene | (6.15*) 
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td 892. (1-2) —1y--a" 
Solutions were obtained for a shallow spherical 
shell having the same H/h and a,/h values as the paraboloidal 
shell. Edge conditions given in Eq. (6.6), the equivalent 
of those in Eq. (6.14), were satisfied at four and five 
collocation points on the boundary of the characteristic 
segment, Good numerical convergence was obtained for as 
Reweaspiourm collocation; points. jhe firstafdve| symmetric 
frequencies calculated using the present theory agree 
closely with the values calculated from Eq. (6.15) as 
shown in Table 6.1, and they agree with the values calculated 
from Eq. (6.15*) to at least four figure accuracy. 
It should be noted that the eigenvalues kK, given 
in Table 6.1 for the spherical shell are identical to 
the eigenvalues given in Table 4.2 for a simply supported 
square plate. According to Koplik and Yu [60], higher 


modal frequencies for clamped shells with circular plan 
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Table 6.1 


Symmetric Eigenvalues and Modal Participation Functions for Simply Supported 
Shallow Spherical Shells with Sliding Clamped Edges having 4, 5 and 6 Sides ° 


col as a 
| = 0.390625 0.1 0.390625 0.5 0.5 
3.14159 
7.0248 
9.4248 


11.327 
12.95 


0.3932(0.4093) * 
0.4515(0.4655) 
0.5645 (0.5758) 
0.7064(0.7154) 
0.8631 (0.8705) 


-0.19869 
-0.07923 
-0.0221 
-0.0290 
-0.0269 





v = 0.33 for. all calculations. 
Edge conditions satisfied at 5 collocation points at 0 = 36°, 30°, 20°, 40°; Gee 

Edge conditions satisfied at 4 collocation points at 0 = Bes 15 oho be S0u5 

Frequencies in brackets were obtained from Eq.(6.15) with i,j =1, 3, 5... This equation was given 
by Reissner [ 50] for shallow parabolodial shells. 
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converge to the frequencies for clamped circular 
plates as the curvature effect decreases with higher 
modes. Consequently, it is likely that the eigen- 
values for this shell would be identical to the plate 
eigenvalues for all modes. Mode shapes computed for 
this spherical shell are identical with the simply 
supported square plate mode shapes shown in Fig. 4.2a, 
as well, Variation of the thinness parameter a/h 
between 10 and 66.667 and the shallowness parameter a/R 
between 0.1 and 0.5 had no effect on the eigenvalues 
and mode shapes, 

In the remainder of Table 6,1, eigenvalues 
are given for spherical shells with 5 and 6 sides 
having different a/R and a/h values. The eigenvalues 
found for these shells are almost identical to the 
values for simply supported plates with the same number 
of sides which are given in Table 4,2, with the 
exception of k, for the 5-sided shell which likely is 
tieer ror ve  esWOUundary COLLOCaULON point SeUs Useduaor 
these shells with 5 and 6 sides are noted in the table 
since fundamental eigenvalues were not found for every 
set used, 

It appears from the results given in Table 6.1 
that eigenvalues for the edge conditions in Eq. (6.6) are 
almost independent of the values of a/h and a/R within 


the limits of thin shallow shells and that they agree 
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closely with eigenvalues for simply supported plates 
with the same boundary shape. 

Modal participation functions, given in Table 
6.1 for uniformly distributed and central point loads, 
were obtained using the eigenvalues given in the same 
table and their associated eigenvectors At where at ard, 
Integral values for the coefficients Drea ae were 
determined using the generalized Simpson's formula with 
46 increments of 0 for the shell with 4 sides and 36 
increments of 0 for the shells with 5 and 6 sides. 
Since eigenvalues and eigenvectors were nearly indepen- 
dent of a/h and a/R values, it appears that for these 
edge conditions the fundamental mode contributions will 
dominate shallow shell response for most uniform dynamic 
loads and that contributions of modes higher than the 
Poa oewl bepe winston tacany for-practical purposes, [he 
reasonably close agreement between the values of (xt /x9) 
and Xi/X, for most modes indicates that the modal 
participation functions computed from Eq. (3.34a) should 
Deprel table satwleast. for these boundary «shapes .aior 
modes higher than the first, modal participation functions 
appear to be almost independent of a/R and a/h values, 


at least for p = 6. 
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Sele Simply Supported Shells, Kraus and 


Kalnins [66] analysed the transverse vibrations of a 
simply supported shallow spherical shell with circular 
Plane icludingevne erfect of longitudinal inertia. The 
shell had a shallowness parameter R/h = 20, a thinness 
parameter a/h = 5.17638 and a Poisson's ratio v = 0.3. 
Employing the edge conditions 


wWla,0) = M %(a,0) = u (a,0) =u (a,9).= 0, (6,16) 
ro sip ) 


they obtained the first twelve symmetric frequencies 
and associated static center modal displacements for a 
load uniformly distributed over the shell surface area, 
Theoretical solutions were obtained for 
shells with regular polygonal bases having l2 and 15 
sides. These shells had the same shallowness parameter 
R/Ameand SPOtsson”s ‘ratio’ vas “the circular’ shell >*however 
the thinness parameter a/h was varied so that the shells 
with circular and polygonal bases had the same surface 
area as calculated from the relations 
A oT ne [2 - / = (a/R) Cow) 
Circular 


and 


p 
A - an? [1 - = [ 1 - 
polygonal 1 
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Tneeqs.) 40.17) and (6.18), A represents the 
circular 


surface area of a spherical shell with a circular 


plan and A represents the surface area of a 
polygonal 


spherical shell with a regular polygonal plan. 


A was evaluated numerically using Simpson's 
polygonal 


formula, Edge conditions given in Eq. (6.2) were 
satisfied at three and four boundary collocation points 
on the characteristic segment. Convergence of the 
eigenvalues and static center modal deflections to two 
or three figure accuracy was obtained for all modes 
except the fundamental. Poor convergence of the funda- 
mental eigenvalues may have resulted partly from the 
Poisson's ratio effect previously noted for simply 
supported plates and partly from errors in the higher 
orders of the Bessel functions required (orders up to 
45 are required for the 15-sided shell with L = 4), 
Frequencies’ and static modal displacements 
obtained for the three collocation point sets are 
compared in Table 6,2 with the values given by Kraus 
and Kalnins. In contrast to the results for three 
collocation points, fundamental modal eigenvalues for 
the shells with 12 and 15 sides obtained for four 
collocation points were slightly smaller and fundamental 


modal deflections were slightly larger than the values 
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given by Kraus and Kalnins, However, for modes higher 
than the first, convergence of frequencies and displace- 
ments with increasing number of sides p to the values 
for a circular boundary was consistent for either 
collocation set, The results given by Kraus and Kalnins 
for the shell with the circular boundary show that the 
largest static center modal displacement due to 
longitudinal inerta effects (i = 3) has approximately 
the same magnitude as the fifth modal displacement due 
to transverse inertia effects (i = 8), Consequently, 
longitudinal inertia effects can be ignored for these 
shells for most uniformly distributed transverse dynamic 
loads, 

The dependence of eigenvalues and modal 
participation functions upon the geometric parameters 
a/R and a/h for simply supported shallow shells is 
shown in Table 6.3. For the thinner shells (a/h = 66.667 
and 62.0459), the third and fourth modes dominate the 
response to uniform loads, while for the thicker shells 
(a/h = 5.2) the fundamental mode dominates the response. 
Of course, increasing the shallowness parameter a/R 
would increase the relative importance of the higher 
modes as well, as pointed out by Kraus and Kalnins [66] 
who compared modal displacements for a shallow and a 


hemispherical shell. In contrast, the relative 
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Table 6.3 


Symmetric Eigenvalues and Modal Participation Functions for Simply 
Supported Shallow Spherical Shells with 6, 12 and 15 Sides 








0.3 
0.26485 
5.297 






















0.390625 
66.667 


0.36355 
62.0459 


0.26265 
5.253 


-0.0496 
-0.016 
-0.038 
-0.033 
-0.001 







-0.0549 
-0.0153 
-0.0194 
-0.0371 
-0.0106 
-0.023 


number of boundary collocation points 


equal surface area equal surface area 
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importance of the higher modes to the response of 
shells subjected to central point loads increases with 
decreasing values of a/h. The decrease of eigenvalues 
i with increasing number of sides p for modes 1 and 2 
and the increase for modes 4 and 5 (except for p = 15) 
for these shells follows the pattern exhibited by plates 
in Tables 4.5 and 4.6. The increase in k, with p 
decreasing from 12 to 6 for mode 3 was not observed for 


polygonal plates, however, 


6.1.3 Clamped Shells. Kalnins [56] obtained 
the nonsymmetric frequencies of free transverse vibration 
of some clamped shallow spherical shells with circular 
plan including the effect of longitudinal inertia, 
Employing the edge conditions 

dw(a,0) 
w(a,0) = ———— =u (a,6) = u (a,6) = 0, (6.19) 
or r 6 
he obtained natural frequencies for shells with v = 0.3, 
a/nie= "10 and &/Ri=_ 0.1, 0.3.and 0.5, 

The symmetric frequencies given by Kalnins 
are compared in Table 6.4 with theoretical values for 
shells with regular polygonal bases having 12 and 15 
Sides. All shells had the same shallowness parameter 
R/h and Poisson's ratio v, however, the thinness 


parameter a/h was varied so that the shells with 
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circular and polygonal bases had the same surface area. 
Edge conditions given in Eq. (6.1) were satisfied at 
three and four boundary collocation points. Good 
numerical convergence was obtained for all frequencies. 

Frequencies obtained for these polygonal 
shells were slightly larger than the values given by 
Kalnins for the circular shells just as they were for 
polygonal and circular plates. The largest frequency 
differences occurred for the fundamental mode and the 
largest a/R value, 

Eigenvalues and modal participation functions 
are given in Table 6.5 for shells of regular polygonal 
plan with 6, 12 and 15 sides and a range of a/R and a/h 
values, The increasing relative importance of the 
higher modes for uniform and central point loads for 
increasing values of a/R is shown for shells with. 12 
and 15 sides and a/h = 10 and for shells with 6 sides 
and a/h = 66.667. The significant increase in the relative 
importance of the higher modes with increasing values of 
a/h for uniform loads and the decrease of their importance 
with increasing values of a/h for central point loads is 
shown for shells with 12 sides and a/R = 0.33. The pattern 
of variation of the eigenvalues for shells of equal surface 
area with 6 and 12 sides is the same as it was for the 


simply supported shells, at least for the first five modes. 
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Table 6.5 


Symmetric Eigenvalues and Modal Participation Functions for 


Clamped Shallow Spherical Shells with 6, 12 and 15 Sides 


0.33 0.33 0.3 0.3 0.3 0.3 : 0.3 
0.390625 | 0.36355 | 0.102331] 0.306993] 0.511655] 0.101481 0.507405 
66.667 62.0459 | 10.2331 10.2331 10.2331 10.1481 10.1481 


-0.0921 | -0.0743 
-0.0741 | -0.0645 
-0.0528 |-0.0526 


number of boundary collocation points 


equal surface area equal surface area 
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Nodal lines for the first four symmetric 
modes of response of a clamped shallow shell with 12 
Sideseandea/hae) 0.702331) a/h = 1052331 and v= 0.3 
are shown in Fig. 6.1. The nodal areas agree closely 
with the values for a clamped plate with 12 sides 
Piven in Fig. 4.e2d, 

Nodal lines for the first seven symmetric 
modes of response of a clamped shallow shell with 6 
Sides and a/R = 0.390625 and a/h = 66.667 are shown in 
Fig. 6.2, The first four mode shapes are quite different 
from the mode shapes for a clamped plate with 6 sides 
Shown in Fig.4,2b. The fundamental mode shape has an 
interior nodal line which was observed also by Chow and 
Popov [90] for a shallow shell with a circular boundary 
and a/R = 0.4422, a/h = 24.04 and wv = 0.3. 

Hoppmann and Baronet [91] studied free 
rotationally symmetric vibrations of simply supported 
and clamped shallow spherical shells silty eeoae en 
boundaries. They found that, for values of a/R of 0.1655 
and 0.399 and a/h of 24 and 48, the interior nodal circle 
of the fundamental mode approached closer to the boundary 


with decreasing a/R(a/h) for the same a/h(a/R), 


6.1.4 Other Edge Conditions, Eigenvalues 


and modal participation functions for some shallow 


spherical shells of regular polygonal plan are given in 
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MODE 4 
FIG. 6.2 Symmetric Mode Shapes for the Characteristic Segment 


of a Clamped Shallow Spherical Shell with 6 Sides 
and a/R = 0.390625, a/h = 66.667 and v = 0,33 
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Table 6.6 for simply supported shells with sliding 
free edges and in Tables6.7 and 6.8 for clamped shells 
with sliding free and clamped edges, respectively. These 
edge conditions, which are given in Eqs. (6.3) through 
(6.5), were difficult to satisfy numerically for the 
lower modes of vibration, Fundamental eigenvalues were 
not obtained for the edge conditions in Eq. (6.5) and 
were obtained for the edge conditions in Eqs. (6.3) and 
(6.4) only for some values of p, a/R and a/h and 
particular collocation point sets, Eigenvalues given 
without associated modal participation functions were 
obtained on the IBM 1130 computer using 10 figure 
computational accuracy while all other values were 
obtained on the IBM 360/67 computer using 16 figure 
computational accuracy, 

Comparison of eigenvalues given in Tales 
6.1 and 6.6 for simply supported shells with sliding 
edges having 4 and 6 sides shows that eigenvalues and 
frequencies for the second and third modes increase 
when the sliding edge is clamped against extension while 
higher modal values are unaffected. This increase 
becomes more significant with decreasing number of 
boundary sides, Fundamental eigenvalues were not 
obtained for the edge conditions in Eq. (6.5), however 
these results indicate that there would be a larger 


increase in the fundamental eigenvalues than observed 
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Table 6.6 ; 


Symmetric Eigenvalues and Modal Participation Functions for Simply Supported Shallow 
Spherical Shells with Sliding Free Edges having 4, 5, 6 and 12 Sides? 


0.390625 0.390625 0.5 0.390625 2 0.5 0.363557 0.5 
66.667 66.667 10 66.667 10 62.0459 





number of boundary collocation points 


1 v= 0.33 for all calculations. 


2 These shells have the same surface area. 
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Table 6.7 
Symmetric Eigenvalues and Modal Participation Functions for Clamped Shallow 


Spherical Shells with Sliding Free Edges having 4, 5, 6 and 12 Sides? 


Oe eae eee eee eee eee 
0.390625 0.5 0.390625 0.5 0.390625 0.5 0.36355? 0.5 
66.667 1 66.667 10 66.667 10 62.0459 





Sie he | Oe eee eee 
number of boundary collocation points 
1 y= 0.33 for all calculations. 


2 These shells have the same surface area. 
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Table 6.8 


Symmetric Eigenvalues and Modal Participation Functions for Clamped Shallow 


Spherical Shells with Sliding Clamped Edges having 4, 5 and 6 Sides? 


Peewee eee. eee | 
0.390625 0.5 0.390625 0.5 
66.667 10 66.667 10 





® 
4e) 
{fe} 
= 


v = 0.33 for all calculations. 

é Edge conditions satisfied at 5 collocation points at 9 = 154 Boe 23°) 10°, Oo: 
3 Edge conditions satisfied at 5 collocation points at 0 = LIAS Clee lees alee) eer. 
. Edge conditions satisfied at 5 collocation points at 0 Bo gore oe, 6257, 4900. 
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for the second and third modal eigenvalues when the 
sliding edge is clamped against extension, Since values 
of kK could not be obtained for any of the simply 
Supported shells with sliding free edges analysed, it 
would be logical to assume that kK is smaller for this 
edge condition set than for any of the other edge 
condition sets analysed, for any particular shell geometry, 
providing that the arguments given in Section 6.1 are 
valid. Eigenvalues for the lower modes of these shells 
were more sensitive to variations in the values of a/R 
and a/h for simply supported shells with sliding free 
edges than they were for simply supported shells with 
sliding clamped edges, 

Comparison of eigenvalues for clamped shells 
with sliding edges having 4, 5 and 6 sides in Tables 6.7 
and 6,8 shows, as well, that the eigenvalues for the 
lower modes increase when the sliding edge is clamped 
against extension. This increase appears: to be most 
significant for the fundamental mode and affects only 
the first two or three modes, depending upon the values 
of a/R and a/h., Eigenvalues for the clamped shallow 
shells with sliding clamped edges agree to two or three 
digit accuracy with eigenvalues for clamped plates with 
the same boundary shape for modes higher than the first. 


Eigenvalues for the lower modes are more sensitive to 
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changes in values of a/R and a/h for clamped shells with 
Sliding free edges than for clamped shells with sliding 
clamped edges, 

Nodal lines for the first seven symmetric modes of 
response of a clamped shell with sliding free edges having 
6 sides are shown in Fig. 6.3. These mode shapes are quite 
similar to the shapes for a clamped plate with 6 sides 
shown in Fig. 4.2b for the first five modes only, 

Mieeresults given in Tables 6,1, 6.6, 6.7 and 6.9 
indicate that the two stress function edge conditions ee = 


ave = 0 can be satisfied numerically more accurately for 
the lower modes of vibration than can the edge conditions 
ee = ioe Oke lyevecause the former edge conditions 
are associated with higher eigenvalues than are the latter. 
edge conditions. The accuracy of satisfaction of both sets 
of edge conditions in the stress function is dependent upon 
tnewvalues Ofea/R and a/h, especially for the edge conditigns 
aie = ee = 0, and decreases with increasing value of the 
PEeUuct .Ol/a/ Rh and <a/7h;, | 

For polygonal ane with sliding edges and as few 
as 4 sides, these stress function edge conditions usually have 
little effect on eigenvalues for modes higher than the third. 
Consequently, higher modal eigenvalues can be calculated from 
the uncoupled plate frequency equations in transverse displace- 


ment for these edge conditions. This uncoupling of the trans- 


verse displacement and stress function solutions does not occur 
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FIG. 6.3 Symmetric Mode Shapes for the Characteristic Segment 


of a Clamped Shallow Spherical Shell with Sliding 
Free Edges having 6 Sides and a/R = 0.390625, 


a/h = 66.667 and v = 0.33 
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for the simply supported or clamped shallow shells 
studied in Sections 6.1.2 and 6.1.3, which satisfy the 
stress function edge conditions Sees = €., = 0, until 
the fifth or sixth mode for shells with 6 sides and the 
fourth mode for 12 sides, 

Values of modal participation functions for 
shells with sliding edges having 6 and 12 sides indicate 
that the fundamental mode governs the response of these 


shells to uniform loads while contributions of modes 


higher than the third are relatively insignificant. 


Grae Edge Condition Satisfaction 


In the numerical solution of the transverse 
vibrations of plates and shallow shells with regular 
polygonal boundaries, edge conditions have been satis- 
fied approximately using the boundary collocation 
technique. Using this technique, edge conditions were 
satisfied exactly only at discrete boundary points rather 
than along the entire boundary. The effect of the 
resultant numerical errors in the prescribed edge con- 
dition functions will be small in the interior for a 
stable solution, according to St. Venant's principle, if 
these errors or edge disturbances are small relative to 
maximum interior functional values, 

It was found that the accuracy of satis- 


faction of a prescribed edge condition for a particular 
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function depends not only upon the number and location 
of the boundary collocation points and the degree of 
deviation ofarnempoundary from avcircle butjalso upon 
the interior spacial variation of the function and the 
functional derivative order when expressed in terms of 
the transverse displacement or stress function solutions. 

Numerical satisfaction of homogeneous edge 
Gondicvions 15 1limstrated in Figs. 6,4) through 6.7 for 
the first three modes of vibration of a shallow spheri- 
cal clamped shell with hexagonal plan having a/h = 66,667, 
a/R = 0.390625 and v = 0,33. Eigenvalues and modal 
participation functions were given for this shell in 
Table 6.5. 

Plots of ratios of prescribed 'vanishing' 


~ W ~ ~ 
edge functions w(6,94), “- C60)", Fg (P00) and ¢€ 
Vv 
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(p, 98) 


ss 
ow 
wee F 

: Se ( vs max 


Bn se. 5) » respectively, for 3, 4 and 5 collocation 
SS max 


to maximum interior values (w) 


ax max? 


points on the boundary of the characteristic shell 
segment show that edge condition satisfaction generally 
improves with increasing number of boundary collocation 
points. Of course, since each solution is found by 
inversion of a boundary equation matrix whose order 
depends upon the number of boundary collocation points 
used, there is a practical limit to the number of 


collocation points which can be satisfied. This limit 
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BOUNDARY 
COLLOCATION POINTS 


=——=55 at 0 = 2°.15°. 30° 


— — 4at 8 =2°,15° 25°, 30° 
—— 4at 8 = 4°,14°, 22°, 30° 
Sat 8 = 2°, 8°, 15°, 23° 





DEGREES 


Fig. 6.4 Numerical Satisfaction of Vanishing Transverse Edge 
Displacement for a Clamped Shallow Spherical Shell 
with 6 Sides and a/R = 0.390625, a/h = 66.667 and 
y= 0.33 
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ae BOUNDARY COLLOCATION 
POINTS 

—-——- 3 at 8 = 2°, 15°, 30° 

— — 4at 8 = 2°, 15°, 25°,30° 

—-— 4at@ = 4°, 14°, 22° 30° 


5 at 8 = 2°,6°15° 23° 30° 
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DEGREES 
Fig. 6.5 Numerical Satisfaction of Vanishing Normal Slope for 


a Clamped Shallow Spherical Shell with 6 Sides and 
a/R = 0.390625, a/h = 66.667 and v= 0. 33 
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Fig. 6.6 Numerical Satisfaction of Vanishing Tangential Edge 
Shear for a Clamped Shallow Spherical Shell with 6 
Sides and a/R = 0.390625, a/h = 66.667 andv= 0. 33 
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BOUNDARY COLLOCATION 
POINTS 


-10°3 2 at Os 2°, 15° 
— —3 atQ= 2° 15°25° 
—— — 3 at@= 4°, 149, 22° 
4 at Q@= 2°, 8215°23° 
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DEGREES 


FIG. 6.7 - Numerical Satisfaction of Vanishing Tangential Edge Strain for a 
Clamped Shallow Spherical Shell with 6 Sides and a/R = 0.390625, 
a/h = 66.667 and v= 0.33 
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is dependent upon the relative 'conditioning' of the 
boundary equation matrix and the computational 
accuracy used in the numerical solution. 

Plots of edge condition satisfaction for 
the first mode show that the degree of satisfaction 
generally decreases towards the boundary corners for 
equiangular spacing between collocation points. Thus, 
best edge condition satisfaction is achieved in general 
by decreasing the angular spacing between collocation 
points towards the boundary corners. In this example, 
the set of 4 boundary collocation points at 6 = 2°, 15°, 
2D jeoUcegaverssmaller maximum deviations for ali 
functions except displacement than did the set of 4 
Polrtemawevete lo A222 4930". Consequently, for 
decreasing number of boundary sides and associated 
increasing boundary deviation from the circular shape, 
the total number of boundary collocation points satis- 
fied over the complete shell boundary must be increased 
to maintain the same degree of edge condition satis- 
faction. 

For the 5 collocation point set, the maximum 
edge condition errors for the first 3 modes occur in 
the second mode for the transverse displacement and 
normal slope functions and in the third mode for the 


tangential in-plane shear and strain functions. These 
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errors reflected the relative magnitude and proximity 
to the boundary of interior functional variations. In 
general, the degree of edge condition satisfaction 
expressed in terms of maximum interior values will 
decrease with increasing mode number, 

In this example, maximum edge condition 
functional errors for the first three modes were 0,24 
per cent of the maximum interior value for displacement, 
Weaseper cent’ for™=normal "slope, ~5-°5 per cent for 
tangential in-plane shear and 18.5 per cent for tangential 
strain. Since the maximum error for tangential strain 
occurred at the shell corner where strain was not 
prescribed to vanish, maximum boundary errors for strain 
are not indicative of the accuracy of the boundary 
col. locatwonmmevnoa. However, ‘since slopewris™a first 
order derivative of transverse displacement and in-plane 
shear is a second order derivative of stress function, 
it is apparent that the degree of satisfaction of edge 
conditions tends to decrease with increasing order or 
the functional derivative of the transverse displacement 


BanGescress function solutions. 


6.3 Series Convergence 


The convergence of the series solutions for 
the functions describing each mode of response is 


illustrated for some edge condition functions in Table 6.9, 
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0.38699E-03 
0.38675E-03 
0.38675E-03 
0.38675E-03 
0.38675E-03 


-0.16798E-03. 


-0.16798E-03 


0.18926E-05 
0.19378E-05 
0.19378E-05 
0.19378E-05 
0.19378E-05 


q, = 1 psi, a/R = 0.390625, a/h = 66.667, D = 0.49315 E + 5 lb.-in. 


Table 6.9 


Functional Series Convergence for the Third Mode of Response 
of a Uniformly Loaded Clamped Shell with 6 Sides 


0.38699E-03 
0.38722E-03 
0.38722E-03 
0.38722E-03 
0.38722E-03 


-0.14518E-03 
-0.14489E-03 
-0.14489E-03 
-0.14489E-03 
-0.14489E-03 


0.49707E 01 
0.50840E 01 
0.50840E 01 
0.50840E 01 
0.50840E 01 


0.29105E-05 
0.28885E-05 
0.28885E-05 
0.28885E-05 
0.28885E-05 


position p,6 


0.30796E-03 
0.20447E-03 
0.20508E-03 
0.20508E-03 
0.20508E-03 


0.61413E-04 
0.42922E-04 
0.43360E-04 
0.43359E-04 
0.43359E-04 


0.54364E-05 
0.58404E-05 
0.58556E-05 
0.58558E-05 
0.58558E-05 


0.30796E-03 
0.41144E-03 
0.41206E-03 
0.41206E-03 
0.41206E-03 


0.53186E-04 
0.69199E-04 
0.69579E-04 
0.69580E-04 
0.69580E-04 


0.44755E-05 
0.45561E-05 
0.45645E-05 
0.45644E-05 
0.45644E-05 


0.52606E-04 
-0.72300E-05 
0.68555E-05 
0.55351E-05 
0.55532E-05 


-0.50153E-04 
-0.17342E-04 
-0.16095E-04 
-0.17138E-04 
-0.17106E-04 


0.73852E-06 


0.20617E-06 
0.60507E-06 
0.53542E-06 
0.53312E-06 


0.41145E-04 
0.32116E-05 
-0.11849E-04 
0.53700E-07 
0.13821E-05 


0.36952E-04 
0.12937E-04 
-0.13727E-04 
-0.67057E-05 
-0.52990E-05 


0.28051E 01 
0.11217E 01 
0.35879E 01 
0.92675E 00 
0.66602E 00 


15 437E-05 
«22866E-05 
~61585E-05 
.13056E-05 
.13912E-05 


ne Nace 0.33 
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These functions were calculated for the third mode of 
response of the uniformly loaded clamped shell with 6 
sides and a/R = 0.390625, a/h = 66,667 and v = 0,33, 
Eigenvalues, eigenvectors and modal participation 
functions were obtained using 5 collocation points on 
the characteristic segment boundary. Consequently, the 


A and iz, given in Eqs. (3.27a) 


series solutions for n 
and (3.27b), and their functional derivatives were all 
truncated with n ranging from 0 to 4, 

Finite ‘exact! solutions were obtained at the 


1 and Et from the zero-order terms only, 


origin for n 
since all higher order terms vanish. Functional 
derivatives of ni and Ei were difficultstoeobtain 
numerically in a small region about the origin. The 
number of terms in each series required to obtain a 
specified functional. accuracy increased with radial 
distance from the origin, as these results show, Com~ 
parison of functional series convergence at p = 0,84, 

6 = 0° (0.027 from the boundary) and p = 0,96, ®@ = 30° 
(0.020 from the boundary) in Table 6.9 shows that 
convergence in the boundary region was less rapid towards 


the corners. Convergence of derivatives of the 7A 


and 
Ei functional series appeared to be slightly less rapid 


than the convergence of the functional series, themselves, 
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For this example, series solutions 
converged within two or three terms to an accuracy 
sufficient for most practical applications over most 
of the shell interior, Results did not converge so 
rapidly near the boundary, especially in the corner 
region, however the degree of accuracy obtained with 
five functional series terms is probably at least 
comparable to the consistency of results which would 
DesolUcal neceinepractice, In géneral, the number of 
series terms necessary to obtain comparable functional 
accuracy increased with decreasing number of boundary 
sides p. It. was felt also that the number of series 
terms required would increase with shallowness and 
thinness parameters a/R and a/h, respectively, although 
results for the shells analysed in Section 6.1 did not 


indicate any Significant differences for these parameters. 


6.4 Auxiliary Edge Conditions 
The auxiliary edge conditions given in Eqs. 
Goel and tg, 23)eetorn the ifunct ions z* (0,8) and h” (0,6) 


which were assumed in the form 
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in Eqs. (3.25a) and (3.25b), were satisfied using the 
boundary collocation technique, The DRES IBM 1130 
digital computer was employed for these calculations 
using 10 figure computational accuracy. The results 
are summarized in Table 6.10 and the coefficients 
oe Beet. pis 


pn pn pn 
cnrough! O1lo. 


sand F's” are given in Tables 6.11 
pn 


Convergent solutions for n*(p,0) were not 
obtainedurfor all edge conditions, especially for smail 
values of p. Better results might be obtained for these 
cases by increasing the accuracy of the computer 
calculations, thereby permitting larger auxiliary edge 
condition matrices. to be inverted accurately, or by 
choosing some different functional form for A” (p,6) 
which will converge more rapidly. 

The solutions obtained for £*(p,6) can be 
used for both plates and shells with regular polygonal 
boundaries which are subjected to time-dependent edge 
conditions involving transverse displacement, slope 
normal to the edge or flexural stress couple tangential 
to the edge. 

The function z*(p,0) which is associated with 
time-dependent transverse edge displacement n(p,9,¢) is 
given by an exact one term solution in Table 6.10. The 
functions h}(p,0) and h2(p,8) which are associated with 


time-dependent edge stress resultant Fy ( 6s 0st) and 


aay | 
Off MHL BRAC eT .eupinnoss notgenollos yisbaved 

oS ate 

solisluolas seed sot Geyolqme eaw ceduqmoo Tastgkd 
tives ad? ,ysaivoe8 Leneltdstuqmos eusggtt OF sat bass 
cinstolttesos sft bas Of,0 sidst nt besbsemmse ‘ots 
fi.0 soltvaT af asvig sts es, bas , * 4%, * ay *s 
.0f.2 davord? 

son stew (6,9) a Tot enobjuloa Insgtevned 
sme tot yifstoeqas ,edottibaes e_be iis to? bentsédo 
see ror bentaddo ad Jdulm atinest weslgea .q i te) seulay 
satuqmoo sit lo yostvoos Seiy aniiasetont yd eeaeo 

a 

sabe vieliixus tepszel golesinieg ydoted? ,anoltefuolss 
~yistetyeos bedteval eo oF asolidam noltdibmos 
(O,q) 5 cot mot Leaottoau? dmoreTItb emoa antsooHo 
yibtas 1 atom sgtevno0o Liitw doldw 

sd mao (A.q)*% ot bentsatdo enoltuloa ent 


[anogyiog isiugst ridiw elleds be sejsig djod tol Dbeev 


enbs tmebasgsb-amts of betost{due eta dotdw esiisbavod 


eqote ,dnemessiqalb sarevenstd gnivloval enolsibnos © 


tay betleltse oxew ,Cd2S.£) bas (ec.€) vepa at 


fistidegnat ef{qvoo anette Lawxell 1 egbe ent od Lamon 






















Ail2 


AA 9 ss SA h € SA z I 
"4 Wate pues uate ot" wuSTM 9 $Ae/Ue UTM, gu" Wwouare nee ieee St 1 


340N 
















9T°9 aTaeL aes f(a) aI = g ozeym Sud Bee ot & d a] 
t 


Ase a ay (999) 20) To = [(6°o) 2a] 1G = iikoeghs alee 
h+UG Z7°2 z4+ud z : 












ud ud Q=U s 9 
ST°9 eTaeL eas S(n*d) y= g aaaym ‘gud en d a+ Cat | {= (099) U] 0 = C(OS) UIT Sum = [(0fg) WIT 
if iff ceso@l Or, ztud 11 rlseoel 7 ; 2 ; 












deynduoo [Te4t3tp oO€Il War euy SuTSN pauteyqo suoTt4ntos 4ueSreauoo ou 









ud ud ud Q=u 
£(d) = azoeum ‘gud eee d ae et 0 a | Ke = (0°9),y4 
atud z°z Gesbiel BOT I-71 iad 











GN johbiey el qrqnie? Bx: 29 





(e6d),u 














ies) Dirck, OSs 29) oi al areyum* od wa 
Zant Zen z* 1 





(a6d),u 














ud ud ud ud Q=U 
SiO ShGP TD OoSm (ted) Yo= 4A aedeum *ogud me d a + d a] a (0°90) yu 
a¢ud rz ego! 30%, I-11 








ud ud ud 























ud Osu 

eT°9 9TQeL ees f(d) q = q atoaum * oud a 0 a + 0 q ] nC = (969),3 
oe fl Za z+ud z%2z rel Bei T-1 

a pue d [T® dojJ [T = (6°d),3 
ud ud ud ud Q=u 

TI°9 eTqey ees f(d) q= q aroym* gud SE d eet d a | x = (6°90) 2 
Zell aut z+ud 2&2 ud zt T=1 








a pue d [Te doj T 


(e°9),3 


SUOF4FDUOD aspq AdeTTTxny s0J (900) a pue (Oae ye AJOJ SUOTANTOS 








SUOT4TPUOD a3pq ATeTTTxny 





OT°9 STQRL 


‘ 
enolsibnod ogba yislilavs to: (0,4) d bas de.) 


ke | grotdiened sgh 


nolsylos So } a 





°G°0 = & 404 peutezygo sem (Ody 8 jo eoUuesd4eAUOd 


°S°O = “ GuIsn peulezgo o408M “ 


Za 


S@149S }s8eq ‘A jo jUepusdepul eae ysnoyu} Ly 3 4Oy ueATs saenzten A 


213 


LO-4 2°0- 
LO=3 ¢°0O- 
LOzaaG Os 
eOa4et. Om 
00 3LE*O- 


LOzs c-O 
LO-3 2°O 
LO-Z ¢°0 
c0r3 7°0- 
00 30€°0O 


LO=3 1°07 
E0s4a 7.0 
cO-9 9°0- 
e074 S°0 
00 48¢°0- 


LO-3 4L°O 
c0-3 1°07 
c0-Z 8°O 
LO-J7L°O- 
00 44¢°0O 


LO-J 
LO-3 
€074 35.07 
cO=4) cag 
OO0FAc9¢ -0= 


6.02 
L105 


LOs3c<0 
LO-3 «L°0 
LO-3 «L°O 
LO=3 ¢c*0- 
00 49€2°0 


EO=a— 0 
c0-9 «COLO 
0-3 7°0 
c0-3 0¢°0- 
c0-J €€°0 
00 4892°0- 


Nereep Iker 
JOs4 GOs 
20-4 6S9°0- 
€0=5).99 °C 
LO-3 L€*0- 
00 3%¢z°0 


e074 S °0- 
20-3. Ss €°0 
COz4q"  e-0- 
cozy tO 
COst ele Os 
cO-3 20 
00 4z€Sc°0- 


S°O 
2307 
58, 


cO-J 
cO-3 
cO-J 
cO-9 SE°0= 
CO=4 ¢f,°0 
LO-4 S6€°0- 
00 3 S02°0 


* €079 


ZO-J «8 0- 
€0-3  6°0- 
€0-3. _-4°0- 
GO=4nenci0 

€0-3 _7°0- 
€0-3  9°0 

00 340S2°0- 


8°0O 
esOn 
c°O 


€0-J 
20-3 
cO-gs €°O- 
20-3 SL°0 
[0545999 0- 
00 38281°0 


d Sepis jo sequin 


TT°9 eTQeL 


‘6074 


5°0- 
c°*0- 


20-3 
20-3 
7073 eLeO 
COs4 anc 0m 
CO=4 0 ClO 
COs4aCLOs 
Ons umc.G 
€O=3 C20 
00 38672°0- 


ve0 
[0 
£:0 
2°0- 
€9°0 
8L°0- 
SL°0 
LO-3 OL6"0- 
00 32L7L°0 


20-3 
e073 
€0-3 
€0-J 
€0-4 
e073 


00 30000S2°0- 


00 3499991 °O0- 
LO-FEEEEES*O 


t 
‘oO = [(e*°9) 23] T SuOoT4;TpUuoD aspq AXetTyTxny aoj (d) 


% 


ud 
A ie 


Or NOM tM OR 


8 
fe 
9 
ee 
4 
€ 
A 
L 
0 





qa SyUSsTFOTSJe09 








‘11.3 efosT — 2 


es 
= ((0,4)*s}, a= <0 = [(5,8)"9] 4 anotstiqcd oghd yustitxuh 20% oe af 
) fs re ae 





! 











q aebie To vedind 


00 306.0 | Gb 3rs.0 | 00 Sees.0 | 0D 3ccs.0 [oo 2 s08.0 | 00 Sater.0 © 9etar.o eee 
$6.3 §.0-| rO<3gr.0-.| t0-3 G2.0-| MED Te.oe| to-3.aee.0- | 1O-3 ake. | 19-9 OFF.0-) 00 MYRodaT ee 
fo-3 $0 | $0.3 8.0 | fod £.0 | sod 28.0 hege .0 [| c0-3 *.0 | S03 a.0 








13 s.0 } sO-3 1.0) TOS 1.0 gO~3 C820- | SHI 2f.0- | S03 £.0- | 50-3 8F.G- . 
' : : 
135.0 |) fo5 f.0 [ro3 s.0 | 0-3 ed ip md £10 | HS 6,0 [tos 63.0 
G3 7.0-| 503 3.6 | 2 to | o-3 <0. 
| ° 
. . S03 2.0 | 969 6.0 63 EO | 
| . | 0.3 10] 
Cha4 aed 
—-——— a ae ee —— —_———— = - _ 
m ———— Se a eS ee ee 
WR. 0- | Mi Has.0- | OO 38S. 5- | SSERS.O~ | 00 aes oe] on aay vs : 
1S0-3 1.0-/ 26:3 2.0 [S02 $.0 | 50 22.0 | 50-5 0 cI da.0 |) t= .0-4 
r.3 ¢ 0. ca.39 2¢,0- 43° B06 3 Co. 4 rc } 4 A thea | co.3 CG | 
| t.3 <.0- Md AO PTR 1.0) 8-5 Si 3-3. tT. | 60 20 | PJ 5.0 
. . 2 ~ | ra ae | 
> & | 0 = ? r 3 4 * -¢ r< on amd on > ’ 
| M3 fF.0 | oo t o0-|265 £0 
{ : : : | of a to. ao | e0-3 5 | 
: | ey 
i : : : j — so 
. ; 1 SQe3 ri 
\" ———— — a ne _— - ina — ee — 
3 
- iolx an <4 ; ; 7 . 
geizee feed .¢ lo Inevrecesdni @ na aigucdi!: -( , goieu bents?do ater “2 wit 


2.0 = y 10) banietds sow (GegP Site 


Table 6.12 


Coefficients ES" (p) for Auxiliary Edge Conditions L (87(6,6)]= 0, L (a*(5,6)] a | 


0 
1 
2 
3 
4 
5 
6 
7 
8 
7 
ce) 
1 
2 
3 
4 
5 
6 
i, 
8 
9 


0.23577E 00 
-0.661 E 00 
0.287 E:00 
-0.163 £E 00 
Oni 12 al 
-0.8 E-01 
-0.1 E 00 
-0.6 E 00 
-0.2 E 01 
-0.4 E 00 


-0.68245E 00 
0.446 E 00 
-0.166 E 00 
Ont5 ae 00 
-0.9 E-01 
0.9 E-01 
0.4 E 00 
0.2 E 01 
0.5 E 01 
0.3 E 01 


0.3337E 00 
-0.376 E 00 
0.16) E00 
-0.1 £00 
0.7 €E-01 
-0.7 €&-01 
-0.1 E 00 
-0.1 £01 
=O.20 1 Of 
-0.2 €E 01 


-0.5527E 00 
0.25 E 00 
-0.1 E00 
0.9 €E-01 
-0.7  £E-01 
0.6 €E-01 
0.2) E 00 
0.2)).E.01 
0.4 .;E—01 
0.2 .E.01 


Number of sides p 


een re icin ae ee 
eed ee ee ee re | eee 


0.384E 00 
-0.258E 00 
Ont) 2eeo 
-0.5 E-01 
0.2 E-01 
0.7 E-02 
0.1 €E-01 
0.2 E00 
0.4 E 00 


-0.519E 00 

0.177E 00 
-0.8 E-01 
0.4 E-01 
-0.9 E-02 
-0.2 E-01 
-0.6 E-02 
-0.3 £E 00 
-0.5 E 00 


0.4147E 00 
-0.192 E 00 
0.1 E 00 
-0.6 €-01 
0.2 £E-02 
-0.2 £E 00 
O56) ge.00 
On7, GE 00 


-0.507 E 00 
On Sage 00) 
-0.9 £-01 
0.6 €E-01 
0.3 £E-01 
Oo 00 
O;6neme 00 
oe? [Boie) 


0.44E 00 
-0.15E 00 
0.7 E-01 
-0.6 E-01 
-0.2 E-01 
-0.2 E 00 


-0.51E 00 
O.1 E00 
-0.7 E-01 
0.5 E-01 
0.2 E-01 
0.2) E00 


0.46 E 00 
-0.15 E 00 
-0.2 E-01 
-0.3 E00 


-0.52 E 00 
0.10 E 00 
0.2 E-01 
033) E00 


0.47E 00 
-0.8 E-01 
(esi) te<Co}s| 


-0.49E 00 
0.6 E-01 
-0.1 E 00 
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0.48E 00 
-0.6 E-01 
0.1 E 00 


-0.49E 00 
0.5 E-01 
-0.1 E 00 
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Number of Sides p 








zal 








Poisson's ratio 
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,[87(B,8)] = 0 

















































































































































































OF] 5On106 E01 OHNE Oi | OnlS(So | Ol 2 Ol) O24 Ie oi 
Wl O12 Ts oe 0.91 E-01 | 0.97 E-01 | 0.115E 00 | 0.136E 00 
2 | One Teli 0.33 E-01 | 0.414E-01| | 0.561E-01 | 0.72 E-01 
A Osh Tex0H On ONE—Ollu On2O7E=Om ORS E-Oln Oe SnE=01 
AN Opi TESy 0.5 E-02 | 0.172E-01 | 0.237E-01 | 0.32 E-01 
5 He rs, ze, Os 5m E025 OnGME=028| Olds E=On 
<|6 g, Si 6, 0.2 E=02 
2 is : 
iS a 3 8 -0.292E 00 | -0.32 E 00 |-0.34 E 00 | -0.34 E 00 
= a oe =OnC9nE=Ofn | -O}76,E=-OnmMieOL94, E=Oln l= On1 12.100 
a a A -0.30 £-01 | -0.388E-01 |-0.542E-01 | -0.709E-01 
F -0.16 E-01 | -0.262E-01 |-0.38 E-01 | -0.50 E-01 
ie O55 [2-02 | —O- iSO |Sonaeeso)) || 0.5%) EO: 
=Os5 am e=02)| “ORG 2-02) Onl mEE=O1 
oe =0.3 E02 
Number of Sides p 
12 | 15 
Poisson's ratio Vv 
0.0 0.1 0.33 0.5 0.0 0.1 jos = 
o | 0.2260 01 | 0.246€ 01 | 0.313E 01 | 0.401E 01 0.292E 01 | 0.318E 01 | 0.407E 01 | 0.517E 01 
1 0.473 E-01 | 0.521E-01] 0.664E-01 | 0.811E-01 0.346E-01 | 0.387E-01 | 0.510E-01 | 0.605E-01 
2 0.141 E-01 | 0.164E-01) | ©.235E-01 | 0.317E-01 0.988E-02 0.117E-01 | 0.173E-01 | 0.20 E-01 
po Bl @a6us HE=O2 0). 7775 = OSMIMNOL T1101 On 5) E=01 0.501E-02 | 0.587E-02 | 0.863E-02 | 0.5 E-02 
be TORS en E=02) | Oss" EeO2n|Os5 8-02" 1"0%6, E-02 0.29 E-02 | 0.34 E-02 | 0.51 E-02 
B || oe EOE oss" 15-02) || One eRe || Onis Exer Q.11 E-02 | 0.13 E-02 | 0.16 E-02 
(2 |Eaaie 
$I -0.4175E 00 | -0.425E 00 | -0.444E 00 |-0.461E 00 -0.518E 00 | -0.528E 00 |-0.552E 00 |-0.564E 00 
5 1 | -0.377 E-01 | -0.423E-01 | -0.556E-01 |-0.692E-01 -0.285E 00 | -0.323E-01 }-0.437E-01 | -0.526E-01 
251 2 | -0.128 E-01 | -0.151E-01 | -0.222E-01 |-0.302E-01 -0.907E-02 | -0.108E-01 |-0.164E-01 | -0.19 E-01 
F Seen Ghee e202" |=Ok75E=02 Ont E-O1 9-515) 2=01 -0.479E-02 | -0.565E-02 |-0.842E-02 |-0.5 E-02 
i Tee Ons E025 0-36 F-02505 E-026|=0,0m E02 -0.28 E-02 | -0.33 E-02 | -0.50 E-02 
al) 7?) Tatoey oge (exes og] E4e2 Ilo.) 0? -0.11 E-02 | =0.12 E=02  ||-0.16 E02 








* Convergent series for h'(p,6) 





were not obtained for p < 7. 


Table 6,14 





12 Ps. oe 
Coefficients Fr (v)* for Auxiliary Edge Conditions L Ch* (6,6) ] = 0, L Ch*(6,6)] =) 














0.625000 


0.3 


OgSheat 





0.714286 


0.746269 


0.4 


0.833333 





2 rin} 
is given approximately by - 


= 0,18393 exp[exp(1.0531V) J 
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ek 1 A % s 
Coefficients F (p,v)* for Auxiliary Edge Conditions —L[h'(6,6) J =i. L Ch*(5,6)] = 0 
pn Eh 





Number of Sides p 
I 





5 | 6 


Poissonts ratio vy 








[oo | 0.1 0.33 

























































































































































































































































0} 0.1767E 01] 0.1705E 01 | 0.1562E 01] 0.1448E 01] 0.170E 01} 0.163E 01] 0.1476E 01) 0.1363E 01 
1| 0.94 E-01| 0.83 E-01| 0.61 E-o1] 0.47 €-01] 0.413€-01| 0.377E-01| 0.282 E-o1] 0.211 E-01 
2 | Ome) LOW) | -O.7 E=01 || —O.23 0 Eaot =O.1) =O) |/=059) E-02)\=0,8  E-021|-0.70 E=02 
3 0.7 €E-02] 0.8 E=02 0.9 €E-02 -0.2 E-02]-0.4 E-03] 0.1 0213)" E02 
20 4 a/h, (ao) || sO [0 -0.9 E-02]-0.6 £-02| -0.3 Ose) NE=02 
Gi Ove Tor) as 0.2 E-02 1-0.1 E-01]-0.8 £E-02/| -0.3 =0.1  E-02 
6 | -0.1 -0.1 =0.8 im) E-03 =03)) E-01|/-0,21) E-07| |||=00'7 =0.4 E-02 
7 0.2 £-02 }-0.2 £-01}-0.1 E-01 | -0.4 -0.2 

= 8 )|| 0-5 E-02 

2 —— 

alte) -0.357 E 00| -0.2864E 00 || -0.387E 00 | -0.351E 00 | -0.2684E 00 | -0.2052E 00 
1 -0.29 E-01]-0.18 €-01 | -0.296E-01 | -0.267E-01 | -0.180 E-01 | -0,127 E-01 
2 0.19 E-01] 0.24 E-01] 0.8 E-02] 0.7 E-02| 0.7 £E-02]'0.63 E-02 
3 [0.8 a 2021-0. 796 6-02 [10.3 E-02 110.9) E03) -0.6  (E-03)|/-0.11 E02 
4 O.3) | E=02)|) 10-20" E-02)]/ 0.11 E-o1) || 0.4 [E=02)|| 0-35) E02) 052: “E-02 
5 -0.2 £-02|-0.1 -02] 0.2 E-01| 0.9 E02] 0.4 £E-02] 0.2 £-02 
6 Os Gm E=031|(O-S ie E-O3) 0,3) E=O1| Oe 2mE=On)| OsGaue E0211 0.5) uE=O2) 
7 -0.2 E-02]-0.1 £-02] 0.2 E-01] 0.1 E-01] 0.5 £02] 0.3 E-02 
8 -0.2 €-02|-0.9 E-03] 








Number of Sides p 


Poisson's ratio v 
ite lee eae ae 
0.161 E 01 0.155 E 01 0.140 E 01 0.157E 01 0.152E 01 0.139 E 01 0.129 E 01 
0.27 E-01 0.24 E-01 0.18 E-01 0.14 0.18 E-01 0.17 E-01 0.13 E-01 0.99 £E-02 
0.1 E-02 O:T E-02 O53) E-02 0.4 0.1 E-02 0.1 £E-02 0.1 E-02 0. 
0.4 E-02 0.4 E-02 0.4 0.4 ! 0.3 Ons) E=02 0.3 €E-02 0 
0.4 E-02 0.4 E-02 0.5 0.6 0.4 0.5 E-02 0.4 £E-02 0. 
0.5 E-02 | 0.5 E-02 0.6 0.7 0.4 O.4 sE=02 0.4 §E-02 0 


-0.279E 00 | -0.210 E 00 


-0.309E 00 
=0.14 E-01 | -0.13 E-01 | -0.95 E-02 
EOed) E-O2: | =One-Orn|-Onl e-02 
-0.3 E-02]-0.3 E-02]-0.3 £-02 
=Oef E=O2)|=O5) E-O2N=0.5)  E=02 
=0.4 E-02|-0.4 E-02)-0.4 E-02 











Number of Sides p 
12 | 15 
Poisson's ratio V 


Bae Bits: Lee ee EE ee ee ee 


0.155 E 01] 0.149 E 01 0.137 £ 01 | 0.128 E 01} 0.154E 01 | 0.149E 01] 0.136 E 01] 0.127 E 01 


































































































































0.3 €-02| 0.3 E02] 0.2 E-02| 0.1 E-O2]} 0.6 E-03| 0.5 E-03] 0.3 £E-03| 0.2 £E-03 
-0.6 E-03/-0.6 E-03] -0.5 E-03|-0.4 E-03 | -0.7 E-03 | -0.6 

=-0.1 E-04]-0.1 E-04] -0.2 E-04|-0.1 E-03 |-0.3 E-03 | -0.3 

-0.3 €-03|-0.3 E-03|/ -0.2 E-03|-0.2 £&-03 E-03 |-0.2 E-03 


-0.186 E 00 






-0.268E 00 









-0.242E 00 




























-0.2 &-02 -0.6 E-03}-0.5 £E-03 -0.2 £E-03 
0.4 €-03 0.7 E-03|] 0.7 E-03 0.4 €-03 
0.2 E-O4 0.3 E-03| 0.3 E-03 O«2 £03 





0.2 €-03 0.2 E-03| 0.2 E-03 E-03 





* Convergent series for h'(p,®) were not obtained for p < 5. 
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strain creer, yt), respectively, are given by one term 


solutions in Table 6,10, also, when either edge condition 
is used in conjunction with shear ee Gince as the second 
edge condition in stress function. The evaluation of the 
modal participation functions aaa associated with Hato se) 
Trend (5 .5.¢)_) and fas associated with @'(p,8) in Eq. 


(3.31d), simplifies considerably for these time-dependent 


feavescond) Lions. 
For the time-dependent edge condition 
oR ee Eq. (3.34d) becomes 


i dL: 
rei Re o +B 9 +C 6 ; 
3 0 13 0 14 0 15 


x 


Thus os is given simply by the zero-order terms of Xi 


For the time-dependent edge condition 
Semin e; 2) used in conjunction with the edge condition 
Fig6Ps9), Eq. (3.34c) becomes 


Meh Sey 8 eae fie + 1 : 
we a | 0 Pas = ae oe I 


Thus i tn is simply a constant 4 bas camecuere Zero— 

order terms of a where Bein: is given in Table 6.14, 
For the time-dependent edge condition 

F(B,8,t) used in conjunction with the edge condition 


BUD 5 0)% Eq. (3.34c) becomes 


(6.20) 


(Gedy 


(2S 
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is simply four times the zero-order terms 
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CHAPTER 7 


EXPERIMENTAL COMPARISON WITH THEORY FOR A SHALLOW 


SPHERICAL SHELL ENCLOSING AN HEXAGONAL BASE 


7.1 Procedure 

(,1,1 Experimental Shell. "A shallow spherical 
Spun aluminium alloy 65S-T4 shell with hexagonal planform 
was mounted in the rolling section of the DRES six-foot 
diameter shock tube using the two support arrangements 
SnOWiein Figs. /.la and-qeib, Thessheldethickness eh; 
nominally 0.375 in., was checked with calipers at a few 
positions near’ the boundary and appeared to be 0.375 in. 
(270, 0050ing). For Trialil,sthe shelienadeasbase circle 
Cincumschibingediameter 2a of 50.0. ins (270.02 inw) sand a 
Spherical radius R of approximately 64 in., and for Trial 
2 the diameter increased to 50.20 in. (+ 0.05 in.) and 
the spherical radius ranged between 64 are (05 in, Or 
theoretical comparison purposes, the experimental shell 
geometric shallowness and thinness parameters a/R and a/h 
will be quoted as 0.390625 and 66.667, respectively, in 
reference to nominal values for a, R and h of 25 in., 
64 in. and 0.375 in., respectively. Young's modulus E and 
Poisson's ratio v for the shell were assumed to be 0” psi 


and 0.33, respectively. 
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FIG. 7.1 b SHELL TEST ASSEMBLY 
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The present theory neglects the effects of trans- 
verse shear and rotatory inertia (thickness-shear 
effect). Theoretical results given by Reismann and 
Culkowski [68] show that the thickness-shear effect can 
be ignored for symmetric modes of response lower than the 
fourth mode for clamped shallow spherical shells with 
CancuLarenianlormeror a/n o> 30 andi a/R as larve ag .0,5, 
Consequently, neglecting the thickness-shear effect when 
computing the experimental shell's transient response 
likely is valid as long as the lower modes governed the 
response, 

In the first trial, an attempt was made to obtain 
the shell's transient response subject to the edge 
conditions for a clamped shell given in Eq. (6.1) and in 
the second trial, subject to the edge conditions for a 
clamped shell with sliding free edges given in Eq. (6.3). 

For the first trial, 2 x 6 in. aluminium (65S-T4) 
Support blocks were attached with epoxy adhesive to the 
outer 2 inches of the shell surface, These edge support 
blocks were clamped between steel blocks which were 
bolted and welded to a steel circular box-beam, Shell 
in-plane edge motion was restricted further for some 
tests by shimming between the shell edge and a steel 
circular member which had an inner boundary shape 
matching the shell's hexagonal boundary, Also, for some 


tests, 350 lb. of lead weights were attached to the 
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circular box-beam to reduce support ringing. In this 
way, the shell was clamped reasonably well against edge 
transverse displacement and normal rotation, and at 
least partially clamped against edge extension parallel 
to the boundary as well as against normal in-plane 
motion. The shell and mount assembly was positioned on 
its side in the shock tube, resting on a wooden block. 
The circular box-beam was bolted to the 1 inch steel 
front plates as shown in Fig. 7.2. These steel front 
plates were used for the plate experiment and are shown 
in Fig. 5.1. The space between the shell steel outer 
ring and the tube wall was closed with 5/7 Opi nye Stee. 
plate and all icracks were sealed with adhesive tape as 
shown in Fig. 7.3. The six 60 degree shell segments are 
labelled 1 to 6 in this figure for future reference, In 
this way the air blast was reflected completely at the 
Snelietest section, For one series of tests, aluminium 
shims were placed between the 5/8 in, steel plate and the 
shell steel “outer ring in an attempt to obtain better 
edge clamping, 

For the second trial, the 2 x 6 in, aluminium 
edge support blocks were cut away in alternate 1 inch 
blocks, The remaining legs of the support blocks rested 
on the steel circular box-beam, Heavy duty Krylon 


Silicone lubricant was applied to the contacting sliding 
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FIG. 7.3 Front View of Shell Mounted in Shock Tube for Trial 1 
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surfaces. Edge piston-type clamps prevented the shell 
from falling off the circular box-beam when the shell and 
mount assembly was positioned on its side in the shock 
tube. By means of this arrangement, the shell was clamped 
reasonably well against transverse edge displacement and 
normal rotation. Edge extension parallel to the boundary 
WaSmrecUrlcLecepy. tne outer cinch portion of the whell, 
itself, and normal edge in-plane motion was resisted by 
thesinertial forces of the boundary support mass and 
frictional forces. The boundary support mass which 
followed the edge in-plane motion of the shell weighed 
approximavelyss4 lb, including’ the 11l™1b. outer 2 in, 
portion of the shell, while the shell weighed 62 lb. The 
shell and mount assembly was positioned on its side in the 
shock tube, resting on a wooden block as shown in Fig. 7.4. 
The shell assembly was not connected by bolts or welds to 
the shock tube walls in any way. Instead, the shell 
assembly was restricted from transverse motion by four 
clamping angles which are shown in Fig. 7.4. In this 
manner, shock tube expansion and contraction under blast 
loading had little effect on the shell boundary motion, 
The test section was completely sealed against the air 
blast, however the portion between the shell boundary and 
the shock tube wall was not a smooth continuous surface as 
it was for Trial 1. Resistance of the shell boundary to 


edge transverse motion and normal rotation was reduced from 
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FIG, 7.4) Front View of Shell. Mounted in Shock Tube for Trial 
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Proiaieiesince ther l inch steel plates supporting the 
circular box-beam were partially removed as shown in Fig. 
TQ: 

Previously, Walkinshaw, Riley, Siddal and 
Oravas [92] obtained the static response of this shell 
subjected to uniformly distributed loads using the edge 
conditions for a clamped shell with sliding clamped edges 
as given by Eq. (6.4), using heavy steel edge rollers and 
Support clamps. The mounting arrangement used in this 
Static test was considered impractical for tests in the 
shock tube since the shell and mount assembly had to be 
PGs CU onedeone1 Us side. JFurthermore, the inertia ofthe 
heavy boundary roller-support mass likely would have a 
considerable influence on the shell's dynamic response 
characteristics which could not be accounted for using the 
PEESENU THeEOry. 

The shell's spherical curvature was formed by 
spinning to obtain a nominal shallowness parameter a/R of 
0.390625 which was renin the accepted limits for shallow 
Sie lvsmOmOuccdy Nast 0,5... opinning Vas poinvedmouumpy, 
Hossack [69], introduces residual stresses and variations 
in thickness, 

The removal of portions of the aluminium support 
blocks for Trial 2 permitted some stress relaxation in the 
shell with an associated shell ‘flattening’ or an increase 


in R, near the boundary. A plot of the shell profile 
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measuredsalcng one of its radial lines in Fig. 7.6 shows 
Cilsm tlauvening’., The variation of the shell radius of 
curvature, measured along two radial lines with a Cenco 
spherometer, is shown in Fig. 7.7. This spherometer 
measured spherical radius of curvatures over 1 inch radius 
circular segments on the shell surface to an accuracy of 

+ 3 per cent. The radius of curvature measured on the 
convex surface of the shell appeared to oscillate three or 
four times between the apex and the boundary with maximum 
variations of approximately + 25 per cent about the average 
value of R which varied from approximately 65 in. near the 
apex to 70 in, near the boundary. 

Hor ¢rial 1, the average spherical middle surface 
radius of curvature was approximately 64 in., however 
ValilavuLonseiinscurvature over the shell surface were not 
measured. It was felt that the extent of the curvature 
Varltavions was similar to that-measured for Trial 2, 

The effect of surface curvature variations or 
surface ‘rippling’ on the transient response of shallow 
shells has not been studied either experimentally or 
theoretically in the shallow shell literature, as far as is 
known, and measured variations have not even been quoted, 
Okubo and Whittier [74], as a result of measurements for 
the static buckling of clamped spherical shells using the 


nondestructive technique, concluded that the lowering of 
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FIG.7.7 Variation of Experimental Shell Radius of Curvature R 
with Radius p for Trial 2 
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their experimental frequencies compared with those pre- 
dicted by theory might have resulted from curvature 
variations, However, they did not include measurements 
of the actual curvature variations present in their 


experimental shells, 


Intuitively, 'rippling' of the shell middle 
surface should decrease the shell's rigidity against 
in-plane stress and increase it against transverse 
bending stress relative to, the rigidity.of a perfect 
spherical surface. It was felt that the significance of 
this variation-in-curvature effect upon the transient 
response of shallow shells would increase with the 
number and magnitude of the surface 'ripples' and the 
values of shell geometric parameters a/h and a/R. The 
Significance of this effect likely depends upon the edge 
conditions to which the shell is subjected, as well. 

One hundred and fifty Tatnall strain gauges, type 
Cl2-121, having gauge factors of 2.09 (+ 0.5 hb), gauge 
Lesistance Of .tc0 ohms and: gauge Tengin of 0.125 sin ee were 
mounted on both sides of the shell along the radial lines 
Ge= 70°50", 20° and 30° av the positions snown win Pig, 
7.8a. Radial bending and membrane strains were measured in 
one shell characteristic segment at 36 different positions 
including the shell apex. These strains were used in 
conjunction with circumferential bending and membrane 


strains measured at related positions in another shell 

























we sk 
Jy Co 
HE 
‘Sear oh ; A; 
betaqmes seatoneupeT? lednemfteqta Ledd 
c bavivess oved sagion ross vd as 
gon bb yon? ,tsvowoH ,enotislisv 
sk: isv studsvatto [sugos edd Io 
elleta Ladnomiteqxs 7 
\ tuatlieate’ ,ylevivigdwal 
} ae | ge10eb bluods soatiwe 
; yrs 1s sestte ensiqent 
itait ett o wvitelet esetian ynibned 
jfe% ssw JI ,sostive Leolrenqe” 
t *s9%%9 erutsviwo-ni-nolsvalaisv ain’ 
fuow slie wollisde Io seenogqaest 
' gostave odd Yo sbutiagsm base tecmun 
" stitemoeg Ifede to seulsyv 
J gostte etdd to sorsoliingtls 
oe rg oe tds t Liste sat doldw os anoltenos 
STys ,&9: | ce Di sibavd en0 
OnueD . 1 Oo ¢) @0.S to atovosl segues gotived LSL-S£9 
_ 
otew ,.mt @Si.0 to dvgnel syua_ dns emo OSL To conadeleey: 
gentl iabba: oi? gnols Lfede ond to eobie dod a0, bedmom 


vai at nwordte enokiteog edt ts COE bas “08 ® OL 170»  € 


= : 7 an a a. 
ot bewasen stew existe ene tdmem bas _antbasd featbsn 88, 

wee 
et tz leoq. SneneItLo 8& 49 eee svekes sae 


> 
aeseid onew emis 5% ext 
one a aie were | Pie 
ae ee 





















‘+ 
paeer es 


235 


PLAN VIEW OF SHELL'S CONVEX (FRONT) SURFACE 


SECTOR 4 


—=— 
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DEFLECTION GAUGE MOUNT, 
TRIAL 2 
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DEFLECTION GAUGE MOUNT, 
TRIAL | 





Pressure-gauge, Trials 1 and 2 

Accelerometer, Trial 1 

Accelerometer, Trial 2 

Linear potentiometer deflection gauge, Trial 1 

Bently Nevada deflection gauge, Trials 1 and 2 

Bently Nevada deflection gauge, Trial 2 

Radial and circumferential strain gauges front and back, Trial 1 


Notes A 
1. All dimensions are horizontal. 
2. Radial distances refer to middle 


surface. Corrections were applied 
to locate strain gauges. 


Obooe@eosx 


Radial and circumferential strain gauges front and back, Trial 2 


Fig. 7.8a Shell Instrumentation 
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DEFLECTION GAUGE MOUNT, TRIAL 1 (not to scale) 


3/4" Dia. Brass Bolt 















Bentley Nevada Deflection Gauge 


(non-contacting) F 
Linear Potentiometer Deflection Gauge 


(glued to shell surface) 


Aluminum Shell 


Aluminum Support Block (Part C) 


Steel Support Block (Part D) 
Outer Steel Ring 
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Steel Box Beam 
21/2" x 2.1/2" x 1/8" thick 
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Steel Support Block (Part B) 


DEFLECTION GAUGE MOUNT, TRIAL 2 (not to scale) 


Bentley Nevada Deflection gauge 
(non-contacting) 













Aluminum Shell 


Fillet Welds Aluminum Dome Supports 
Spaced 1" 


3/4" Dia. Brass Bolt 








21/2 51 1/2" x 1 1/2" 
Tapped Steel Block 






Fillet Weld 
PLP LL 


Steel Circular Support Beam 
(Part A) 
Fillet Weld, one side 





Note - Refer to Fig. 7.1 for detail 12" x 5" 31.8 lb/ft I beam 


of shell test assembly 


Fig. 7.8b Shell Instrumentation 
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characteristic segment to obtain normal stress resultants 
and flexural stress couples at the 36 positions assuming 
perfectly symmetric shell response. The shell's symmetry 
of response was checked at 4 of these 36 positions by 
measuring radial and circumferential strains in two other 
segments, 

7.1.2 Natural Frequencies, Shell natural 
frequencies were obtained using acoustic excitation with 
sound pressures ranging between 100 and 130 decibels at 
the shell surface. These sound pressures were applied 
over a large portion of the shell surface area. In this 
manner, the suppression of some resonant frequencies, 
which McConnell [93] showed could occur for concentrated 
excitation forces such as those given by shakers, was 
avoided. The shell base plane was placed perpendicular to 
Gietax1 SsvoOtfetherexponential) thorn in order -tovincrease the 
intensity of pressure loading by using standing waves to 
advantage; a procedure pointed out by Hubbard and 
Houbolt [80]. Placing the axis of the acoustic horn 
parallel to the shell base plane would have provided a 
more uniform sound pressure level distribution over the 
shell's front (convex) surface, however the pressure level 
obtained was too low to excite some of the lower shell 
modes. Excitation frequencies ranged from 200 to 3,000 cps. 


Strain and acceleration responses at different positions on 
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the shell surface indicated the various resonance values, 
Nodal lines were established for several symmetric modes 
of response of the shell subject to the edge conditions 
obtained in Trial 1 from the phase changes observed for 
the Lissajous figures for acceleration signals recorded at 
many different positions on the shell's convex surface. 

In Trial 1, three Electro-Voice model 848, 30 
watt, 10 x 12 in. rectangular horn-type speakers with 
Frequency Trangessore)50rv0r 10 ,000%ecps»weresused tosexcite 
Chew sneun symmevricainvy.s IneTrial» ec y= symmetric: excitation 
Wasr ODvalned*with a Geloso’ model 257955200 watt, 29 in. 
diameter circular horn-type speaker with a frequency range 
OrelcUelorl es uoOF cps. 

The sweep-frequency was obtained using a 
Nelson-Ross model PSA-021 spectrum analyser in conjunction 
with a Nelson-Ross model 601 synchro-sweep generator. The 
sound pressure level at the shell surface was kept constant 
over the frequency range for some of the tests by regulating 
the amplitude of the sine wave generated using a General 
Radio model 1569 automatic regulator in conjunction with a 
General Radio model 1565-A sound level meter having an 
operating range of 44 to 140 decibels. The sine wave 
Signal from the regulator was amplified by a General Radio 
model 1308-A power amplifier and then fed to the horn-type 


speakers, 
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In order to maintain a constant sound pressure 
level over the entire frequency range using these 
speakers, the sound pressure level had to be kept at 
approximately 100 decibels. However, some lower shell 
resonances were not distinguishable from system noise at 
this pressure level. Consequently, for many tests 
constant power amplifier output voltage rather than 
constant acoustic pressure level was maintained over the 
frequency range giving pressures as high as 130 decibels 
adveubeom lower trequencies, ..dowever, using. this technique, 
response amplitudes for different modes were not 
indicative directly of the relative modal participation 
factors as they were when constant sound pressure was 
maintained over the frequency range. 

Acceleration signals were obtained using a 
Kistler model 818 piezoelectric accelerometer, Strain 
signals were amplified with DRES designed amplifiers 
having Cains Of 15007025, 000 and frequency nesponse Of 0 
to 20,000 cps. Acceleration and strain signals were 
converted from a.c. to d.c. using a Hewlett-Packard model 
4OOE voltmeter and then recorded on the y-axis of an 
Autograph model 65 two-axis recorder. At the same time, 
the sweep-frequency produced by the synchro-sweep 
generator was recorded on the x-axis. This frequency was 
monitored by a Hewlitt-—Packard model 3734A electronic 


counter, as well. Lissajous figures, obtained from accelero- 
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meter signals, were displayed on a Tektronix model 
502 Oscilloscope. 

The acoustic excitation equipment is shown in 
Fig. 7.9. The Geloso horn speaker and the General Radio 
SOumOmlevel meter iare shown in-Fig. 7.4, positioned on ithe 
convex side of the shell which is mounted in the shock 
tube for Trial 2. A complete instrumentation block 
dlasramelseschown iin Filey 7.10. 

Since the shell edge supports were too stiff to 
be excited acoustically, natural frequencies were obtained 
by striking the shell edge with different rubber and metal 
hammers while recording acceleration and strain signals on 
a Tektronix model 564 storage oscilloscope. 

7.1.3 Transient Response. The shell was 
positioned at the same section in the six-foot diameter 
shock tube as the square plate was tested previously, with 
the central axis of the shell aligned with the axis of the 
shock tube. The characteristics of the shock tube at this 
location are described in Section 5.1.3. For this set of 
experiments, the test section was subjected to air blast 
waves having peak normally reflected overpressures qn and 
positive durations Kes of approximately 10 psi and 32 msec., 
respectively. All blast overpressures resulted from the 
detonation of 1.5 lb. charges of 60/40 RDX/TNT, Higher 
overpressures were experienced at the test section for 


this charge weight than during the plate experiment as a 
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Two-axis recorder 
Oscilloscope 

Spectrum analyser 
Synchro-sweep generator 


Equipment 
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result of modifications to the shock tube compression 
chamber, 

the transient loading function for the plate 
experiment had the form q = q, Q(t), where 
Q(t) Ramu eat (t+) expe - t/t.). However, for the shell 
experiment, the transient loading function was described 
by a more complex function having the form 
GPa qa )ar*(p,t), as a result of the shell surface 
curvature, The assumption of axial symmetry of the load 
was exact for Trial 1 since geometric symmetry was 
maintained over the complete shock tube cross section, as 
shown in Kigvey.3 vin Trial 2, geometric asymmetry was 
introduced at the shell's polygonal boundary, as shown in 
Fig. 7.4, however its effect on the symmetry of the shell 
surface loading function was thought to be insignificant. 

The effect of surface curvature on the reflection 
of an air blast wave is demonstrated in the schlieren 
DNOUOST apis Jl tte, febs and in’ the pressure-time records 
Tiesto lest tor al cyLindrical surface whichenadeagragiiacs 
Of curvature Opel. ii, andathe same Jang lelot opening or 
23 degrees as the spherical shell. This test was conducted 
in the DRES 2.x,J2 in, rectangular shock tube for an 
incident air blast overpressure of approximately 10 psi at 
the shell surface. This shock wave had a peak normally 
reflected overpressure of 26 psi and a shock front velocity 


Ob l/ *300 Tps"at the cylindrical surface. 
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(a) t = 0.017 msec (b) t = 0.140 msec 
Fig. 7.11 Schlieren Photographs Showing Air Shock Front Reflection by a Cylinderical Surface 
ina 2x 12 in. Shock Tube 
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Fig. 7.12 
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The rectilinear incident and curvilinear 
reflected shock fronts, which appear as dark lines in 
scmitenen “photographs, “are shown in Fig. 7.lla ‘at 
t = 0.017 msec. after the incident shock first reached 
the shell apex. The curvilinear shock fronts, 
reflected from the shell surface and the rectilinear 
shock tube walls, are shown in Fig. 7.11lb at t = 0,140 
msec. 

The overpressure magnitudes of the shock 
fronts travelling across the shell surface are a 
function of the spatial variation of the obliquely 
replectedstiock wave magnitude, ~ Theoretically, the 
reflected overpressure magnitude generally decreases 
with increasing angle between the shock front and the 
réflectingvwsurfiacesy rand the relative magnitude’ ofthis 
decrease in overpressure generally increases with 
incident blast overpressure magnitude. Calculations 
from equations given by Bleakney and Taub [94] show that 
the magnitude of a normally reflected overpressure of 
26 psi at the shell apex would decrease by 1.5 per cent 
at the shell edges. Extrapolating from curves given in 
the "Effects of Nuclear Weapons" [95], this decrease 
appears to be approximately 2 per cent. 

The four pressure-time histories in Fig. 7.12, 
which were recorded perpendicular to the cylindrical 


surface, indicate that there was a larger decrease in the 
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magnitude of the initial reflected overpressure with 
increasing angle of incidence of the shock front than 

is predicted by [94] or [95]. However, the pressure 
records are distorted by gauge 'ringing' which makes 
initial overpressure values difficult to estimate, 
Consequently; these records serve merely as an indication 
of the general comportment of the transient blast 
overpressure for this curved surface. These records 

show that pressure disturbances travelling across the 
cylindrical surface were not important after t = 0,4 msec. 
Thiseimpdilesethauysateleast for thisgsheliecurvaturejand 
blast overpressure, the transient overpressure was 
distributed almost uniformly over the shell surface after 
the incident shock wave had travelled past the shell apex 
to the shell boundary, reflected to the shock tube wall 
anduthen’travelied back across thershell surfacento,its 
apex which was located at the central axis of the shock 
tube. 

In the shallow spherical shell experiment, the 
shock front arrived at the shell apex 0.34 msec. before 
reaching the shell boundary corners which were 5.1 in. 
further downstream, Reflections from the shock tube wall 
crossed the shell apex approximately 3 msec, after the 
incident wave first reached the shell. However, 
according to equations given by Bleakney and Taub L9O4], 


the normally reflected overpressure at the shell apex of 
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10 psi would decrease in magnitude by only 0.4 per cent 
at the shell boundaries. Curves given in [95], Fig. 
3.71b, p. 147, show that the maximum decrease in the 
reflected overpressure for this shell would range 
Beuween Upardeo. sper cent, respectively, for incident 
overpressures between 2 and 5 psi which are associated 
with normally reflected overpressures between 4,25 and 
11.5 psi, respectively. It was felt that pressure 
disturbances travelling across the shell's surface 
would not be important for such small variations in 
reflected overpressure magnitude. Consequently, radial 
variations in the pressure magnitude G(p) were ignored, 
The importance of the radial variation in the 
time of shock arrival of up to 0.34 msec. which can be 
represented by a function Q*(p,t), depends upon the 
values-of the shell's governing periods of response, 
Generally, for most transient loads it is assumed that 
structural response can be determined with sufficient 
accuracy using values of the load impulse which are 
averaged over time intervals as large as one-tenth of 
the natural period of a one-degree system (see, for 
example, Biggs [96]). Therefore, it is reasonable to 
assume that if the shell governing periods of response 
were greater than 34 msec. (i.e., shell governing fre- 
quencies were less than about 300 cps), then the function 


Q*(p,t) could be replaced by the function Q(t). Other- 
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wise, radial variations of the loading time function 
should be included in the theoretical analysis, For 
normalemode theory. the. shell's response.to.a load 
function Q*(p,t) can be obtained only approximately, 
This is done by superimposing the shell's response to 
loads applied on different portions of the surface 
associated with different loading time functions Q(t), 

which combine to give an approximate representation of 
Q*(p,t). 

Trievnesshe ligt: lois leandeceapressures .0is— 
placements, strains and.,accelerations were recorded _at 
the locations shown in Fig. 7.8 with three Ampex 14- 
Urackesctape recorders, gusing the techniques, describedvin 
Section 5.1.3 forjthe plate experiment. .Face-on 
reflected overpressures were measured at four positions 
around the shell boundary using DRES shear tube piezo- 
electric pressure transducers. Shell central aeniaees 
ment was measured with a Bently Nevada non-contacting 
inductance transducer. This transducer had a model 
3500N detector-driver ee model 308 proximeter coil 
with a range of 0,250 inches. Central displacement also 
was measured with a linear potentiometer in Trial 1. The 
transverse accelerations of the mounts for these central 
displacement transducers shown in Fig. 7,8b were measured 
with an Endevco piezoelectric accelerometer, For Trial 2, 


in-plane boundary motion of the shell was measured with a 
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Bentley Nevada transducer having a model 302 proximeter 
coil with a range of 0,100 inches, For Trial 1, in-plane 
boundary acceleration was measured with an Endevco 
accelerometer, Strains were measured with 150 strain 
gauges at 36 unique r,O middle surface positions in 
Trial 1 and with 34 strain gauges at 9 unique positions 
rote are, , 

The instrumentation block diagram shown in 
Fig. 5.5 for the plate experiment is representative of 
the instrumentation set-up for the shell experiment 
except that the Bentley Nevada displacement transducer 
Signals were conditioned with DRES strain gauge bridge 
balances and amplifiers, Also, some data recorded on 
magnetic tape for Trial 2 were read directly into the 
DRES IBM 1130 digital computer, This was accomplished 
with an SYS model 2113 data acquisition and control 
system having a Miniverter 8 channel, 13 bit analogue- 
to-digital (A-D) converter multiplexer, The signals, 
recorded at 60 ips on magnetic tape for frequencies of 
up to 20,000 cps, were played back on the converter at 
3 3/4 ips giving a 16 to 1 speed ratio. The 8 tape 
channels were read sequentially at a rate of 10,222 
samples per second. Thus, 20,444 samples per second of 
actual shell response were digitized and stored in the 
computer for each tape channel. Assuming two samples are 


required to define one response cycle, this gave an A-D 
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digitizing capability for response frequencies up to 
Ole care DSz 

Other magnetic tape data were digitized with 
the chart reader and plotted in the manner explained for 
Ene piace .expeminents in Sect tonn5.ts 3eiwixpentmental: 
strain data and theoretical and experimental sectional 
resultants were represented graphically for 8 msec. with 
responsesinequency reproduction «wf jup> tosh ,222ncps fon 
Trial 1 and for 10 msec, with frequency reproduction of 


uplitos J0 J000meps, fomiinrial 2% 


7.2 Theoretical Shell Enclosing an Hexagonal Base 


Theoretical modal frequencies and static 
central deflections, computed from eigenvalues and modal 
participation functions given in Chapter 6, are given in 
Laven eer Ore une) i. 0St Six or seven symmetric modes of 
vibration of uniformly loaded shells with hexagonal plan 
satisfying the six different edge condition Eqs. (6.1) 
Through (6.6). Fundamental eigenvalues @ould not be found 
when satisfying Eqs. (6.3) through (6.5) for values of 
a/R and a/h of 0.390625 and 66.667, respectively, as 
explained previously. 

It should be recalled that the evaluation of 
the modal participation functions was simplified by 
integrating over the shell plan area rather than the 


actual surface area. However, the errors introduced by 
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using this simplification were not large since the shell 
plan surface area of 1624 in,” is only 3.29 per cent 
smaller than the actual surface area, assuming a = 25 in. 
and R = 64 in, 

Comparison of the results given previously in 
Table 5.2 for the natural frequencies of a square plate 
with those in Table 7.1, reveals that the lower modal 
frequencies of a shallow shell do not vary as much with 
the type of edge conditions satisfied as they do for a 
plate, However, changing the geometric parameters a/R 
and a/h causes significant variations in the shell lower 
modal frequencies, as shown in Table 7.1 when a/R changes 
from 0.390625 to 0.25. Comparison of values of central 
modal displacements indicatesthat for a shell with 
a/R = 0.390625 and a/h = 66,667 satisfying edge condition 
Hgs.|(6.1) or (6.2), modes 2 through 5 govern shell 
response with mode 3 dominating, while it appears that if 
any of edge condition Eqs. (6.3) through (6.6) are 
Satisfied, modes 1 through 3 govern response with mode 1 


dominating. 


7.3 Theoretical and Experimental Results for Trial l 
7.3.1 Natural Frequencies and Mode Shapes. Some 


of the shell's natural frequencies for Trial 1 are 
indicated inekigs., (13 through 7,15 by the various 


resonances recorded with strain and accelerometer gauges 
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Notes 


1. Five 1 in. wedges equispaced 
on each boundary side. 





2. No lead weight at boundary. 
3. Two 30 watt speakers 
4. Constant voltage excitation. 
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Notes 

1. Wedged all around. = 
2. 350 Ib. of lead weight. 

3. Three 30 watt speakers. 

4. Constant voltage excitation. 
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1. Five 1 in. wedges equispaced 

2. 350 Ib. of lead weight at boundary 
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1. No wedges at boundary 
2. No lead weight at boundary 
3. Two 30 watt speakers. 
4. Constant voltage excitation 
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1. Wedged all around. 
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FIG. 7.13. VIBRATION OF ACOUSTICALLY EXCITED SHELL AT APEX FOR TRIAL | 
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Notes | 

1. No wedges at boundary. 

2. No lead weight at boundary. 
3. Two 30 watt speakers. 

4. Constant voltage excitation. 
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(B) TRANSVERSE ACCELERATION 


FIG. 7.15 VIBRATION OF ACOUSTICALLY EXCITED SHELL AT BOUNDARY OF SECTOR 3 
FOR TRIAL |. 
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for acoustic excitation frequencies ranging between 200 
and 3000 cps. Although the abscissa scales have not 
been calibrated in these figures, the excitation 
frequencies have been noted for most of the significant 
resonances in cps and are accurate to approximately 

* 5 cps. Ordinate scales were not calibrated either 
since sound level pressures were not kept constant for 
most tests. Modal response amplitudes cannot be related 
between graphs even approximately since different signal 
amplifications may have been used. 

The large number of higher resonances obtained 
may have resulted in part from increases in the sound 
level overpressure occurring when standing waves were 
BOvmeimavpeGrvelerenbeecxCavation Trequencies, a0! course, 
shell longitudinal inertia modes were excited in addition 
to the transverse inertia modes predicted by the theory. 
However, the larger number of resonances obtained at 
positions away from the shell edge, as shown in Fig. 7.14d 
for example, indicates that antisymmetric as well as 
Symmetric modes may have been excited. The occurrence 
of nonsymmetric resonances could have resulted from 
asymmetries in the sound pressure distribution as well as 
in the shell geometry and edge fixity. 

The attempt to attain better edge clamping by 
shimming and adding mass at the shell edges in order to 


simulate the edge conditions in Eq. (6.1) produced some 
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changes in the modal frequencies, The resonant frequency 
ranges obtained for the various clamping attempts are 
summarized in Table 7.2 along with comparable modal fre- 
quency ranges for shells with a/R = 0.390625 and 

a/h = 66,667 satisfying edge condition Eqs. (6.1), (6.2), 
CO. 4) and? (6,6). 

The experimental resonant frequency ranging 
between 400 and 460 cps was not predictable for this shell 
geometry and these theoretical edge conditions. Oniashvili 
[88] pointed out that eigenvalues k, must be real as the 
frequencies of a shell which can resist bending are always 
higher (thaneits membrane frequency, Theoretically then, 
the lowest possible shell frequency occurs for k = 0 and 
has the value w = 1/R VEE? Assuming R was 64 in., the 
minimum shell frequency for any edge conditions would be 
495.1 cps. Therefore, it is felt that the occurrence of 
this low resonant frequency resulted from either coupling 
between the shell and its mount or an effective increase in 
the shell radius of curvature R, 

Results trom Jriali2, which will be given in 
Section 7.4, indicate that the occurrence of this low 
resonant frequency resulted from properties of the shell 
itself rather than the mounting system. However, it should 
be noted that the dominant resonant frequency of the shell 


aluminium support blocks was about 230 cps when the steel 


yes 























yonsupert saanoset ed? .astomevpett Isbom odd nt segnans, 
o1s udqmest¢a gatgmelos evottsv sd tol beatsido segnst 
-97° Labo eldsisqmoo nviw amofa S.¥ efdeT at bestzsmmue 
bis @Sa0er.0 = fie dtlw elfedte tol segnst Yoneup 

($,d) ,(£.3) .apt mobdinmoo sgbe gaiylalise yad,00 = d\s 


a 


aiigto: yoneypert tnenoset ([sdnemiteqxs eT 

(fede etd¢ 10% sidazolpbetg Jon saw ego 03 bers 004 neewted 
iLiviesinO .anekttbnos egbe Igeoltetoedt seedd Ons Yaseniosg 
ed? 25 {set sd Jeu ? aculevaegls gadd suo bedatog (88) 
vewie 91s untbaed satest apo dotdw Ileds & To sstonsupett 
edd ylisoltetosdT ,yoneipetl sasidmem att csi? tsfiaid 
bas 0 ® x tot etiopo yomeupet?t leds sldlesog geswol ons 
ait amt id sew A goimwesA » yNay A\l = w» eulav ods eed 
ad bivow aenotsibaos exbs yas tol yoneupert Ilede muminia 
40 sonotruosve ond sade 2fe% ef 2 ,etotetesn? .ago Lae 
enliques 19itte mort betiwest yoneuypet? .Jnsnoses wor atti? 
qt sesevont evicoette as to tevom ett bas Ilefe odd asowsed 
fl esudeviwo to eutbss ciode edt | 


at hovig ed titw dotaw ,S feieT mor? pee if 
| ay tar 
wol efi? to exnsrsmpso od ted? svaotbat Per 
7 te CA a 
12, bed, onsuper? Insnoe 
ASAE 96? 3p Rape INNS Oe | sage ge epee a 


ages! 












Position on Experimental Shell Surface where Resonances were 
Measured 


Table 7.2 


Shell Experimental and Theoretical Natural Frequencies for Trial 1 
















Mode Numbers and Frequencies of Theoretical 
Shells with a/R = 0.390625 and a/h = 66.667 
which have Frequencies Comparable to 
Experimental Shell 


Edge Conditions Satisfied 
Eqs. (6.1) and (6.2) J|Eqs. (6.4) and (6.6) 
Mode Frequency | Mode Frequency 
Range (cps) | Range (cps) 
410 to 460 Lei cl lc hall leans A 
536 to 554 11,2 LST te fa25 
511 to 568 


535 to 560 (2 resonances) 
720 to 760 | 722 to 760 ‘| 3 720 to 797 
898 to 948 944 to 1000 | 4 B59) £0 955 








Approxi mate Exper imental 
Modal Frequency Ranges 
ae 


x 
2 1064 to 1242(4 resonances)| 5 1101 to 1239 | 5 1099 to 1226 
1340 to 1395 1398 to 1537 | 1372 to 1518 


eine 
eee eee eee — ———_—_——_——— ————— ee 
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higher resonances 
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clamping blocks and wedges were removed, These blocks 
were excited by a hammer striking the shell surface 
rather than by acoustic sound pressure, in order to 
obtain a measureable resonance, 

PueeerheCorve sneer. radius Of Curvature Rk 
would have to be approximately 80 in. in order to pre- 
dict a natural frequency as low as 400 cps using Eq. (6.6) 
Whrie TUewouldsnave torbe 100 in, to predict this low 
Brequency uUsing™Eq.*(6.1). “It was felt=that Rk” approxi- 
mated 64 in. much more closely over the shell surface in 
trvaueeechierm 167010 in irial-c although actual variations 
Were NOU measureda ia wiolglemmererore, it is postu 
lated that the resonant frequencies between 400 and 460 
cps were introduced as a result of surface 'rippling' 
which was discussed in Section 7.1.1. 

The dominance near the apex of a resonant 
frequency of approximately 540 cps and the occurrence of 
another resonant frequency of approximately 500 cps 
indicate that there might have been a homogeneous elastic 
restraint condition relating F (Ps 9) and F_ (6,0), and 
also that edge condition Eqs. (6.4) and (6.6) were more 
applicable than Eqs. (6.1) and (6.2). 

This closer approximation of the experimental 
shell edge conditions to the roller rather than the 


clamped edge conditions is confirmed by the mode shapes 
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shown in Fig. 7.16. These mode shapes agree more closely 
with the theoretical mode shapes of a shell satisfying 
edge condition Eq. (6.6) shown in the same figure than 
they do with the mode shapes of a shell satisfying edge 
condition Eq. (6.1) which were shown in Fig, 6.2. 
Complete experimental nodal lines could not be obtained 
tie Gasca wmeapeciaily for those near’ the boundary 
where the vibration amplitude was smallest. The 
asymmetry of some of these mode shapes probably is 
indicative of the nonuniformity of the edge condition 
satisfaction around the boundary. 

7.3.2 Transient Pressures. The average peak 
normally reflected overpressures qd, and associated 
durations toe as measured Tor eachy oigthesseven I beibs 
RDX/TNT shots in Trial 1 by three shear tube pressure 
gauges located around the circumference of the shell 
mount assembly, were 10.3 psi (+ 0.3 psi) and 32 msec, 
eceieecum respectively, 

It was not considered worthwhile to include 
variations in the time of arrival of the shock front 
over the shell surface since theoretical solutions could 
not be obtained for some of the roller-edge conditions 
which were thought to best simulate the experimental edge 
conditions. Therefore, the transient load history over 


the complete shell surface was assumed to have the form 
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a) Experimental Resonant Frequency — 540 cps b) Experimental Resonant Frequency — 750 cps 
Theoretical Frequency, Mode 2 —— 541 cps Theoretical Frequency , Mode 3 —— 720 cps 
c) Experimental Resonant Frequency —- 900 cps d) Experimental Resonant Frequency — !098 cps 
Theoretical Frequency , Mode 4 — 859 cps Theoretical Frequency , Mode 5 —— I!099 cps 
Notes 
t+ Three 30 wott speakers used for excitation. ——— Experimental nodal lines 


2 Ib- of id weight at . 
ee ad nal amet sed Sel eal Mac ----—-— Theoretical nodal lines for a simply supported 


3: Wedges all around shell edge. shell with sliding clamped edges and a/R = 
0.390625, a/h = 66.667 and v= 0.33 


FIG. 7.16 Some Experimental and Theoretical Shell Mode Shapes for Trial | 
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q(t) = qo (1 = Vee exp c= t/t.) ee bees toe GES) 
where q, and ve were -10,0 psi and 32 msec., 
respectively. The negative overpressure sign signifies 
that the transverse or en coordinate was chosen as 
positive inea direction outward from the front or 
convex surface of the shell. The value for dg was 
ehosen as = 10.0 psi rather than -10.3 psi to account 
for the slightly lower obliquely reflected overpressures 
Which would be initially experienced over much of the 
shell surface in comparison to the normally reflected 
overpressure experienced by the pressure gauges. 

Typical experimental blast overpressure and 
SHelimesir alle necordssare Shown) in Pigs. /.1/. 

7.3.3 Transient Response. Solutions for 
transverse displacements and sectional resultants, 
formulated in Appendix A.2, have been computed for a 
clamped shell and a simply supported shell with sliding 
clamped edges. These solutions, obtained by summing 
the first seven symmetric modes of response, will be 
rererrea to as ineory, 1 and Theory, 2, respectively. 
Geometric parameters used in the calculations were 
Bee edi, eet rOd its and, hb € 0, 3/58 1g tne sex peri 
mental sectional resultants were computed by relating 
radial strains in sector 1 to circumferential strains in 
sector 3 using the assumption that each of the six 


sectors of the shell surface responded identically. 
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This assumption of symmetric response appeared to be very 
good on the basis of comparison of the radial strains at 
p = 0.64, © = 0° and p = 0,80, © = 0° in sectors 1 and 4 
and the circumferential strains at p = 0,64, © = 30° and 
p = 0,80, 9 = 30° in sectors 3 and 6. 

Damping was not included in the theoretical 
response calculations as it was felt that the edge 
conditions of Theory 1 and Theory 2 did not simulate the 
experimental edge conditions closely enough to make its 
inclusion worthwhile, Damping effects, although they were 
not negligible, did not appear to be as significant as 
for the plate experiment. 

Dimensionless theoretical static and maximum 
dynamic transverse displacements are plotted along the 
Padia lines ©  =f0- Fin fic.) 7.108. 4 loeistinteresting to 
note that, because jot beating, the ratlogor maximum 
dynamic-to-static displacement exceeds 2 in the region of 
the shell apex. In contrast, this ratio would not 
usually exceed 2 for plates since for most transient 
loads the fundamental mode dominates plate response. The 
inclusion of viscous damping for the shell would cause a 
greater reduction in the maximum values near the apex 
than in the values closer to the boundary since maximum 
values near the apex occurred several oscillations after 


maximum values were attained near the boundary. 
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TRANSVERSE DISPLACEMENT 


Clamped 


Shell 





a Static Load 


-—O-—- 0-— Blast Load, cy = Ol; qe a -10 psi, a = 32 msec, 


FIG. 7.18 Shell Theoretical Static and Maximum Dynamic 


Transverse Displacements on @ = 0° 
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Dimensionless maximum dynamic experimental and 
theoretical sectional resultants are plotted along the 
wadial ‘Wines 9 = 0%, 10°, 20° and 30°°in Figs.) 7.19- through 
7.22. The edge conditions which were satisfied by the 
experimental shell at 0 = 0° can be readily inferred from a 
comparison of the theoretical and experimental values of 
the sectional resultants shown in Fig. 7.19 near the shell 
boundary (p = 0.866, © = 0°) since the r,® coordinates 
coincide with the v,s coordinates along 9 = 0°, Deviations 
from the prescribed edge condition values in this figure are 
larger for Theory 2 than Theory 1 primarily because the 
numerical solution for Theory 1 was obtained using a 
collocation point at 0 = 0° while the solution for Theory 2 
did not use a collocation point at 0 = 0°, It is apparent 
from a comparison of values of the radial stress resultants 
that the boundary of the experimental shell satisfied more 
closely the roller edge condition Fy (By) = 0 than the 
clamped condition u_ (8,8) = 0. Eliminating edge extension 
by setting E(B 59) = 0 is equivalent to setting 


~ 


Fe 5639) =v ELA Comparison of the experimental radial 
and circumferential stress resultants reveals that the shell 
boundary was not well restricted against edge extension, 
Comparison of the values of the radial flexural stress 


couples indicates that the experimental shell boundary 


satisfied the condition $¥(5,6) = 0 quite well. Setting 
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STRESS COUPLES 







—P—P—) Ss Experiment 


Theory, t4= 0 


an Os On Clamped Shell 


---@---@--- Simply Supported Shell with. 
Sliding Clamped Edges 


FIG. 7.19 Shell Maximum Dynamic Experimental and Theoretical 


Sectional Resultants on 0 = 0° for Trial i 


a“ 


ot. al 


7 tw rw} oe 






Tes 


weiseo TeaN7 






oS 






noteewoert? ene Jargnesagsd olny, ’ 
~~ hn 
i fate? 424% “v= se ww Adan ueOn Lae hoe 


n 7) “ar 4 


268 


STRESS COUPLES 











0-8 1-0 





STRESS RESULTANTS 









—pb—P —_ Experiment 





Theory, oy = 0 
—-o— «+-o— Clamped Shell 


--@---@-- Simply Supported Shell with 
Sliding Clamped ldges 




















FIG. 7.20 Shell Maximum Dynamic Experimental and Theoretical 


Sectional Resultants on 0 = 10° for Trial 1 
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STRESS COUPLES 











STRESS RESULYANTS 









—p ——_- PD — Experiment 





Theory, b= 0 
—-O—-O—.- Clamped Shell 
--@-- @e-- Simply Supported Shell with 


Sliding Clamped Edges 

















TIG. 7.21 Shell Maximum Dynamic Lxperimental and Theoretical 


Sectional Resultants on 6 = 20° for Trial 1 
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COUPLES 
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FIG. 7.22 Shell Maximum Dynamic Experimental and Theoretical 


to} 
Sectional Resultants on @ = 30 for Trial 1 
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w(6,0) = 0 results in the relation between the edge flexural 


stress couples that May (5,0) = = Vv mM. (6,9). Comparison of 


Ss 
the values of the radial and circumferential flexural stress 
couples at the boundary shows that the condition of edge 
transverse displacement was satisfied quite well for the 
experimental snell and for Theory 1 but not, quite iso, well 
for Theory 2. 

It appears from this discussion of edge condition 
satisfaction that a theoretical solution for a clamped shell 
with sliding free edges would provide better agreement with 
the experimental results than the predictions of either 
Pie ory eeOre  NeOmy oa HOWwever,. i tmprobably would.pe 
necessary to use elastic restraints for both edge conditions 
in stress function in order to obtain the theoretical 
solution most compatible with the experimental results, at 
least with reference to the satisfaction of homogeneous edge 
conditions. 

Theoretical values of the twisting stress couples 
Me = - Moo showiein Figs,../s20 and {.21 areslarcersin 
magnitude for the edge conditions of Theory 1 than for those 
of Theory 2, however .they are smaller than the values of+the 
flexural stress couples Myre and Mor for either set of,edge 
conditions. Theoretical values of shear stress resultants 
F.9 = Fon shown in the same figures are generally smaller in 


magnitude than the values of the normal stress resultants 


oe and F although the magnitude of f approaches that 
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oh Foo near the boundary for Theory 2, 

Experimental and theoretical central transverse 
displacements vs. time are compared in Fig. 7.23. The 
experimental displacement record is not too reliable as it 
is composed of some gauge mount vibration as well. Experi- 
mental and theoretical sectional resultants vs. time are 
compared at six different positions on the shell surface in 
Figs. 7.24 through 7.29. It is apparent that there were 
experimental frequencies present which are lower than any 
predictable theoretically for the shell geometry assumed, 
Although the experimental stress couples responded at 
Prequerncaes similar to those for Theory 2 over most of the 
shell's surface except near the boundary where they 
responded at a frequency of about 400 cps, the stress 
resultants were governed by frequencies ranging between only 
350 and 450 cps which generally decreased towards the 
boundary. It is from these results that the hypothesis 
regarding the significance of variations in the middle 
surface radius of curvature in the form of ‘surface rippling’, 
CUS CURSCUMMiInoe CULO si...) ALOSC. 

A spectral analysis was performed on some of the 
combined membrane plus bending strain data recorded in the 
radial direction on the front (convex) surface of the shell 
using a digital computer program developed at DRES which 
employs the fast Fourier transform, Record lengths of 50 


msec. were analysed at a rate of 20,776 samples per second 
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CENTRAL DISPLACEMENT 
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FIG, 7.23 Shell Experimental and Theoretical Central 


Displacement vs, Time for Trial 1 
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FIG: 7-24 Shell Experimental and Theoretical Sectional Resultants vs: Time 


at p= 0,@=0 ond p= 0-06, O= O°, respectively, for Trial | 
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STRESS COUPLES 
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FIG: 7-25 Shell Experimental and Theoretical Sectional Resultants vs’ Time 


ot p=0:24, 9= 30° for Trial | 
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STRESS COUPLES 
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FIG: 7-26 Shell Experimental and Theoretical Sectional Resultants vs- Time 
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STRESS COUPLES 
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FIG: 7:27 Shell Experimental and Theoretical Sectional Resultants vs: Time 
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FIG: 7-28 Shell Experimental and Theoretical Sectional Resultants vs: Time 


at p+0-72, 9*20 for Trial | 
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FIG: 7-29 Shell Experimental and Theoretical Sectional Resultants vs: Time 
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giving a maximum frequency reproduction of 10,388 cps, an 
accuracy of + 20 cps and a resolvability of 40 cps, Con- 
sequently, the presence of frequencies associated with the 
first two modes of vibration of the experimental shell for 
sliding edge conditions, theoretically separated only by 
approximately 40 cps (see Table 7.1), might not be 
detectea" from’ the” analysis of a single strain record. 

Some of the more significant response frequencies 
in cps can be summarized in groups of decreasing importance 
with respect to strain at different positions on the shell's 


surface as follows: 


p = 0, 0 = the 1) S449 

PoC O aT ocee Looe 

3) 488, 894, 264 

4) 610, 1687, 2378, 346, 122 
Biba evel GytO Eel OCs wl oma! 0.7) 00529 

2) 366, 122, 752, 1138, 305, 468 
ee ace pan a 07 

2) 142, 529, 285, 752, 894 

2) e920 35m 610, gllos 


Frequencies underlined were considered to arise 
from '‘'beating' of certain modes, For example, the fre- 
quencies of 122 and 142 cps likely represent beating between 
the modes with frequencies of 407 and 529 cps, and 407 and 


549 cps, respectively. These results show that there were 
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dominant response frequencies at HOT, 529 and 549 eps 
(+ 20 cps). This analysis of the spectral representation 
of the shell's forced response characteristics agrees 
reasonably well with the acoustic excitation measurements 


SHOW Dairies. f2l3, through 7.15 


7.4 Theoretical and Experimental Results for Trial 2 

7.4.1 Natural Frequencies. Some of the shell's 
Maturalerrequenciesstorsirial 2 are indicated in Fig." 7.330 
by the various resonances recorded over the shell surface 
with strain gauges for acoustic excitation frequencies 
Taverne veyweenecU0 land 1200 cps, 

The governing frequency of response over most of 
Chershiell surface cecurred at 530 cpS with other significant 
resonances occurring at approximately 580, 710 and 950 cps. 
HNOWeVver wos riot = lawien™a-relatively low trequency wor 
around 400 cps dominated response near the shell boundary, a 
frequency of around 300 cpsedominated response in the 
boundary region for this trial. This low resonant fre- 
quency must have Seon from proper Vvics OlmuNeuscacam 
1tseli rather than of the mount, since in-plane edge 
rigidity of the mount had been virtually eliminated for this 
brial. 

Accurate experimental mode shapes could not be 
obtained for this trial since the boundary was not well 


restrained against free transverse vibration, However, the 
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(A) FRONT CENTER STRAIN (B) FRONT RADIAL STRAIN AT r=6", 6=45° 
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(C) FRONT RADIAL STRAIN AT r=i4" @20° (D} FRONT RADIAL STRAIN AT r=21" @#0° 
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FIG. 730 VIBRATION OF ACOUSTICALLY EXCITED SHELL FOR TRIAL 2 
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mode shape for the 530 cps resonant frequency appeared to be 
nearly the same as that for the 540 cps resonant frequency 
Oe Fis elewnichewas: snown in fic, 7,16a. 

7.4.2 Transient: Pressures. Transient over- 
pressures on the shell as measured by two shear tube pressure 
gauges were unchanged from Trial 1 and consequently were 
represented theoretically by Eq. (7.1) where thas and See were 
=10,0 psi and@se msectryerespectively. 

7.4.3 Transient Response. Experimental trans- 
verse displacement and sectional resultants are compared in 
this section with the theoretical solution for the simply 
Supported shell with sliding clamped edges used previously 
fOtemeu Tabet) ae, 

Experimental and theoretical central transverse 
displacements vs, time are compared in Fig. 7.31. This 
experimental displacement record is considered to be 
reliable since the gauge mount vibration was almost com- 
pletely eliminated for these tests. The maximum experimental 
deflection-to-shell thickness ratio of 0.05 was well within 
the linear response range. The experimental edge in-plane 
displacement normal to the boundary shown in the same figure 
was almost as large as the transverse central displacement. 
This edge displacement appears to have undergone considerable 
damping. 

Since strains were measured at only nine unique 


r,9 positions, it was not possible to plot variations of the 
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EDGE NORMAL IN-PLANE DISPLACEMENT 
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FIG: 7-31 Shell Experimental and Theoretical Central Transverse 
Displacement and Experimental Edge In- Plane 


Displacement vs. Time for Trial 2 
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values of the sectional resultants along radial lines as 
was done for Trial 1 when strains were recorded at 36 
unique positions, However, experimental and theoretical 
sectional resultants vs. time are compared in Figs. 7.32 
through 7.37 at the same six positions on the shell surface 
at which they were compared in Trial 1 in Figs. 7.24 
Cheoucheinco, respectively 

Comparison of experimental sectional resultants 
at the shell apex for two different shots in Fig. 7.32 
reveals, that there was some sensitivity, especially in 
values of stress resultants, to the changes in edge support 
conditions which could occur between successive shots. 

Comparison of the values of stress resultants 
near the boundary in Figs. 7.29 and 7.37 shows that, 
although more edge extension was permitted in Trial 2 
than in Trial 1, the normal boundary force did not change 
significantly for the two different. support arrangements, 
at least at 6 = 0°, Comparison of the values of radial 
stress couples in these two figures shows that the support 
arrangement of Trial 2 introduced an edge stress couple of 
opposite sign to the couple which occurs for a clamped 
edge with dw(p,0)/av = 0, Comparison of the radial and 
circumferential stress couples near the boundary in Fig. 
7.37 indicates that the condition of vanishing transverse 


edge displacement was not well satisfied since 
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at p= 0, 6+ 0 and p= 0-06, 8+ 0, respectively, for Trial 2 
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My .(B 9) # - vM_4(6,0). Transverse edge displacement was 
not eliminated completely for Trial 2 since the shell 
boundary was supported against transverse motion only in 
one direction over alternate one inch intervals as shown in 
Fig. 7.la. It appears from these results that the theory 
for a simply supported shell with sliding free edges would 
provide a better approximation of the experimental edge 
conditions than Theory 2, 

Thesresulte in Fies., 7.31 through 7.36 show that 
the shell transverse displacement and sectional resultants 
responded at frequencies comparable to the theory after a 
ReveivivetyacLowelirSu.duarver cycle of response, It is 
CoNsigeredsthavy this initial frequency response, noviced 
particularly for stress resultants which responded Perey 
aueerequencaes: even slower than those for Trial 1, was too 
much lower than any possible theoretical frequencies to be 
accounted for completely by the variation over the shell 
BurLace Olethe Gime Of arrival of the shock ‘front. This 
time of arrival had a maximum variation of only 0.01 
divisions on the time scale of these figures. Near the shell 
boundary the governing response frequency was about 300 cps 
asyseen fromaresultsein Wigs. This frequency: was 
predicted from the acoustic excitation tests, 

Some of the more significant response frequencies 
in cps obtained by spectral analysis can be summarized in 


groups of decreasing importance with respect to front radial 
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strain at different positions on the shell's surface as 


follows: 
pec 0 sso) = 0° 1) 529 
ye soe 
Bee OlO ee OLG 4 | 
4) 976, 264, 833 
Dom Ose seo nol, 1090. 1159 leu. 
1362 
OeeelC LO LO0 elOU) se aa5 
6 = 0,56, 6 = 0° 1) 285 
2) 325, 529, 61, 468 
Sih Saas 
We MOMS» MS SrA HME He NO 
Sie ols, tes mo oH 
p = 0.84, 6 = 0° 1) lye 
2 SSI aE) 


3) oe Ot 5 L034 5 GLU 
4) 386, 854, 1016, 1403 


Frequencies underlined were considered to arise 
from beating. This spectral representation of the shell's 
forced response agrees reasonably closely with the acoustic 


excitation measurements shown in Figs. 7.30a, 7.30c and 7.30d, 
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CHAPTER 8 
OBSERVATIONS, CONCLUSIONS AND RECOMMENDATIONS 


8,1 Observations and Conclusions 

8.1.1 Theoretical Analysis of Plates, Some of the 
more important results obtained in this thesis from the 
theoretical analysis of the linear dynamic symmetric response 
of thin elastic plates with regular polygonal boundaries can 
be summarized as follows: 

PeepoLeenvalues, jeligenvectors, mode shapes and 
modal participation functions, many of which previously were 
Mie veilos tem UmunewLateravure are Obtained for la number of 
different simply supported, clamped and elastically clamped 
polygonal plates subjected to uniformly distributed and 
central point loads, 

e, the trends apparent in the changes of eigenvalues 
and mode shapes with boundary shape often agree quite closely 
with the trends observed in these functions for shallow 
spherical shells, In fact, in many cases eigenvalues and 
mode shapes as defined by the nodal lines are identical for 
these two structural types for higher modes when they have 
the same boundary plan. 

3.. the pattern of change in the modal frequencies 
Ws with number of plate sides p was observed to be the same 
for both the simply supported and clamped edge conditions. 


It was found that for polygons of equivalent surface area, 
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increasing p has the following effect on the associated 


ut 
then decreases; W increases up to p = 7 and then decreases; 


Velie smote WwW .- Ww deereases ; Ww. increases up to p = 4 and 


W increases up to p = 9 to l2 and then decreases; and a 
increases, 

It was observed that the maximum modal frequency 
occurs for the polygon with the least number of sides for 
which the nodal lines do not deviate significantly from 
circles concentric with the origin, and is generally associated 
with the polygon having the largest difference between nodal 
areas responding in the same direction. This observation 
can be explained physically in the following way. As p 
decreases, tle relavive proportion of "Inactive’ =corner area 
Pncreaseo, 950 vit, as long as the nodal lines maintain their 
Crrcutatmonape ticle as alr associated ancrease in che plate 
stiffness and modal frequency. However, the deviation of a 
nNeoadalr sine Irom the circular shape in order’ to adjust to the 
polygonal boundary shape, is accompanied by a relaxation in 
plate stiffness. The increase in the nodal area function 
relacvive to that’ for a Aveodiar boundary is specified by the 
function ms, and=iasea maximum value Or aboul Super cente1or 
any of the first five modes, 

This correlation between the maximum modal 
frequency and the mode shape, vs. the number of boundary 
Sides, likely occurs for shallow shells as well, although 


enough shells were not investigated numerically to be certain. 
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4, The products of terms of different orders n 


i 
0 


integration over the plate surface area of the square of 


were omitted in the evaluation of xy which involves 
the displacement eigenvector, The omission of these 
terms considerably simplifies the numerical evaluation of 


oe and, it appears, still provides reasonably accurate 


a 


values of x5 


for boundaries with as few as four sides. 

5. An approximate method is introduced whereby 
mocgadsparticapacion functions are calculated for different 
polygonal plates using known values of these functions for 
ayparticular polygonal plate with large :-p.,. In»this 
method, changes in the areas between the nodal lines 
pertaining to the different boundary shapes are related to. 
changes in the modal participation functions. Reasonably 
accurate modal participation functions are obtained for 
plates with more than three sides using this simple 
technique in spite of the approximations made with regard 
to the modal geometry. 

6, The results indicate that plate response will 
be represented adequately for practical requirements for most 


uniformly distributed load transients when contributions of 


modes higher than the second are ignored, and for most central 
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point load transients when modes higher than the third 
are ignored, 

{. It was observed that the static central 
deflections of plates having equal surface areas increase 
with the number of boundary sides p while the associated 
values’ of the fundamental frequencies decrease, Likely, 
these variations can be attributed to the decrease in the 
relative proportion of ‘inactive' corner areas with 


increasing p. 


Oaliee Comparison of Theoretical Predictions and 
Experimental Results for a Square Plate, A comparison of 


theoretical predictions with the experimental results 
obtained when an elastically clamped square plate was 
subjected to blast loading revealed the following points: 
1, satisfactory theoretical predictions for 
flexural stress couples can be made for the Stiles oles 
surface even when some membrane strains are present, 
providing edge spring constants and viscous damping 
coefficients can be estimated with reasonable accuracy. 
2. High levels of damping which increased 
rapidly with stress level up to a limit, occurred in the 
experimental model, This nonlinear damping likely was 
attributable to structural interface damping in the plate 
mounting system, Air damping was considered to be 


relatively insignificant for most blast overpressures, 
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3. The fundamental frequency of the experimental 
plate initially decreased with increasing load even more 
than could be explained by the large values of viscous 
damping observed, It is felt that this experimental frequency 
decrease might have been an effect of edge conditions in 
transverse displacement and normal slope being time- 
dependent rather than homogeneous, since it could not 
Desasnonlinearseifiect ..heoretically, nonlinear .eomet ry 
effects will cause an increase in the frequency with load, 
However, this nonlinear hardening effect should have been 
practically unnoticeable in this experiment since the 
maximum displacement-to-thickness ratio was always less 
GHAanBUe > « 

4, The presence of membrane strains in the experi- 
MGhiwa apace omnOueaccOounLed for in the presentelinear theory 
nor is it predictable by other nonlinear theories ,at these 
deflections. Measurements of plate boundary motion indicated 
that these strains could have arisen from in-plane time- 
dependent boundary loading. It is possible that the in-plane 
boundary loading frequency was not constant and could have 
accounted for the initial plate frequency decrease with 
increasing load which was discussed above. However, the 
experimental results seem to confirm, as predicted by Dawe 
L47], that the effects of bending and membrane action were 
essentially uncoupled at the relatively low levels of 


membrane strains encountered, 
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8.1.3 Theoretical Analysis of Shells. Results 
obtained from the theoretical analysis of the linear dynamic 
Symmetric response of thin shallow spherical shells having 
regular polygonal planforms can be summarized by the 
following points: 

1. Eigenvalues, mode shapes and modal participation 
functions are obtained for a number of shells with different 
polygonal plans subjected to uniformly distributed and central 
point loads, for a wide range of geometric parameters a/R 
and a/h and edge condition sets simulating variations of 
clamped and roller edges, The only results previously 
aveltaplew asetaneasels known;econcern®the» free vibrations 
of shells with square and circular plan and the forced 
hespenserol shells twith circularsplanjyfor just avefewrof 
these edge conditions. 

2. Numerical problems which arose in the determin- 
ation of fundamental eigenvalues Ki, appear to be associated 
with an insensitivity of the circular frequency w to changes 
in k which occurs when k* is small in comparison to 
12(1 - v?)a*/(Rh)*?, Consequently, values of k, were not 
always found for small values of ky or large values of the 
product of a/R and a/h, The results show that smaller values 
of ky are associated with a relaxation in edge restraints and 
an increase in the number of boundary sides p, as would be 


expected, 
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3. EHigenvalues, eigenvectors and modal participation 
functions for shells with roller edges satisfying F (6,8) = 0 
are unaffected by edge conditions in stress function and its 
derivatives, usually after only the first two or three modes, 
This uncoupling of the stress and transverse displacement 
functions for the higher modes is associated with an 
insensitivity of the eigenvalues, eigenvectors and modal 
participation functions to changes in a/R and a/h, cIt ds 
associated also with an interchangeability of eigenvalues 
and mode shapes for shells and plates having the same 
boundary shape and subjected to the same edge conditions in 
transverse displacement and its derivatives. 

4, The importance of the contributions of modes 
higher than the first to shell response and the response 
sensitivity to values of a/R and a/h, both increase 
significantly for shells satisfying €,,(6s9) = 0 when the 
second edge condition in stress function is Z, 
rather than F (6,8) = 0, Physically, this implies that 
even though edge extension is completely eliminated by 
satisfying ES G70) = 0 “in conjunction with w(p,6) = 0 
so that in-plane motion normal to the edge will be very 
small in comparison to that associated with the edge 
conditions F (848) = F (6,8) =40), ‘formexampiel wathe 
additional 'wrinkling' permitted when Fy (658) = 0 is 


satisfied as the second edge condition rather than 


Fg (by) = 0, drastically changes the shell response 
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characteristics, 

5. The edge conditions LD = ug, (6,6) = 0 
were satisfied in an attempt to simulate the conditions 
uy (5,6) = me ome) = 0, The equivalence is not exact since 
the function du (6,6) /av is not necessarily zero for either 
sev ‘of edge’ conditions. ~ However’, the numerical results 
indicate that effects of satisfying either edge condition 
set may be equivalent for p = l2 and that differences might 
NOUmOC eat etcer Olmoiatler > 

It was discovered at a late date that expressions 
for in-plane displacements in terms of the stress function 
had been obtained for the static case by Reissner [97] and 
Fettahlioglu [98]. Therefore, it should be possible to 
check the effect of satisfying the two different edge 
conditions for any polygonal boundary. 

6, The numerical satisfaction of a prescribed 
homogeneous edge condition using the boundary collocation 
technique generally improves with increasing number of boundary 
sides p, with decreasing functional derivative order of w 
or F, and with increasing number of boundary collocation 
points L up to a maximum number, the value of which 
increases with the digital computational accuracy used, 
Best satisfaction is achieved when boundary collocation 
points are located with closer angular spacing towards 
the boundary corners where larger magnitude errors always 


tend to occur. However, it would be very difficult to 
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specify the one best spacing for each boundary shape and 
number of boundary collocation points used, since any 
particular collocation point spacing will give best overall 
satisfaction for only one edge condition function. This is 
true since the degree of edge condition satisfaction 
depends upon the interior spatial variations of the edge 
condition functions, and these functions are linearly 
independent, 

These observations and conclusions, made 
specifically from observing the edge condition 
savislacvLton of a clamped shell *should*apply-ain’general, 

Although the application of the boundary 
collocation technique is not new, even today few 
investigations have been made specifically into the most 
suitable choice of collocation points, as Collatz [99] 
Deluveasoul in 1959, [tC has been stated for the static 
solution of these shells that the number and location of 
thew collocation points is not of paramount significance 
[92]. However, in the present dynamic solution these 
factors are not always negligible. For example, both the 
number and the location of the points influenced 
considerably the ability to determine ky for some shells 
with roller edges. Consequently, in many cases the number 
and sometimes the location of the collocation points have 
been given along with the associated numerical results. 

Lo, Niedenfuhr and Leissa [100] have suggested 


that overall edge condition satisfaction can be improved 
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by obtaining more edge condition equations than unknown 
integration constants and then using a least squares 
procedure to minimize errors in the integration constants 
solved. This approach does not appear to be valid since 
it is impossible to simultaneously minimize the errors 

in two linearly independent functions with one condition, 
In addition, since the practical application of the 
boundary collocation technique depends upon the truncation 
of the functional series, it might be more advantageous to 
use the extra edge condition equations to compute 

eda vronalmvermns @inethesnunctional series (rdependingsupon 


the rapidity of the series convergence, 


Crete H Comparison of Theoretical Predictions and 
Experimental Results for a Shallow Shell with an Hexagonal 


Base, Comparisons of theoretical predictions with the 
experimental results obtained when a shallow spherical 
shell was subjected to acoustic and blast loading for two 
different sets of boundary restraints, revealed the 
following points: 

1, shell governing modal frequencies for uniliorm 
loads as measured from acoustic excitation tests provided a 
good estimation of the governing frequency components of the 


shell's response to blast loads, 
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2. Comparison of experimental and theoretical 
mode shapes along with natural frequencies provided a 
good indication as to the experimental edge conditions 
which actually governed shell response, The comparison 
of mode shapes is considered to be essential, since ‘for 
certain edge conditions some modal frequencies are not 
"well separated' and could be overlooked in an analysis 
of the frequency spectrum alone, Closely spaced modal 
Préequenciessoccur, ifonnexample, for ithe first; twoemodes 
COimomcue) Pywith robbentedges ealAlso; frequencies sof 
different modes of the same shell satisfying different 
bheoretcalvedgce conditions rcouldrbepalmost identical, 
making it difficuit=to ascertain. from a frequency 
comparison alone the actual experimental edge condi- 
Ghonsesatusiied;,* Mhis edifficulty arose for the -expeni= 
mental shell as the theoretical frequencies of the first 
three modes for the simply supported shell are almost 
identical tothe sfrequencies of modes <2; 3 and 4, res- 
pectively, for the simply supported shell with sliding 
clamped edges, 

3. Comparison of theoretical predictions and 
experimental results for the response of the shell to 
blast loading revealed only fair agreement, Some experi- 
mental frequencies of response occurred which were much 
lower than any predictable using the present theory. 


These lower response frequencies became more dominant 
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toward the shell boundary and affected the response 
of the stress resultants more than the stress couples, 

The discrepancies between theory and experi- 
menueoccurred | homiseveralyreasons, oiFirst ,» theoretical 
solutions could not be obtained for the edge conditions 
best simulating the two different sets of experimental 
edge conditions. Second, spatial and time variations 
in the blast overpressure over the shell surface were 
not accounted for in the theory, Third, damping was not 
included in the thoeretical calculations, although it 
Goum@amnavesbeen- quiveumcadPiv:, mefourth “the experimental 
shell geometry deviated from the properties of perfect 
symmetry, uniform shell thickness and constant middle 
Surface radius of curvature,all of which were assumed 
in the theory. Fifth, experimental edge conditions were 
not uniform around the boundary and were nonhomogeneous. 

The occurrence of the low experimental res- 
ponse frequencies was considered to be attributable 
principally to variations in the shell's radius of 
curvature and not to mount vibrations since these low 
frequencies occurred even when in-plane boundary 
rigidity had been virtually eliminated, 

Hae biatappeans «phat at» isi notepracti calsaito 
Simulate the clamped edge conditions in stress function 
u (6,8) = ape) = 0, at least for a shell having a 


surface area as large as that of the experimental shell 
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(i.e. 1670 in.*) which is subjected to uniform loads. 

It is felt that in most practical applications the edge 
conditions in stress function will more closely 
approximate the sliding free edge conditions 

F.(6,5) = F.(6,6) = 0 than the clamped edge conditions. 
8,2 Recommendations 

Certain limitations in the theory presented for 
the dynamic response of these plate and shell structures 
became apparent from the numerical and experimental 
ativesvulgations. However, the major “effort extended vo 
investigate the ability of the theory to predict the response 
of the wide range of structural shapes for which it is 
applicable to the variety of transient load distributions 
which are permissable, naturally curtailed detailed 
investigations in particular areas, Consequently, there 
are several problems which are recommended for future 
Puvesvieavion,. Nese are. 

1. Modal participation functions xi should be 
evaluated without neglecting the products of terms of 
different orders n, for modes having nodal lines which 
deviate significantly from circles concentric with the 
origin. This should be done for all modes of response 


of plates and shells with 3 and 4 sides, 














@ _ ' 
9 to 

,elbaol amotiny 0% betoetdve et dotriw Coat oar’ 9st) 

om hie 
eno tk, hiqas [gotdoetq teom al dsa7 tiet at 
: at 

[ftw nolsonwt seetsea at enotilbneo 

100 enbe sort gilitlea esd? sivamixorqas 

. 


. = P 7 ce - 4 ~ ~~ « 2 ae 
- : ‘ oo | 3 $0 a Sif: 0 . \ ¥,4 4 ay a (8 +4) y7 


anol *sboom rooeh Ss. 8 
a 


oe 


) eric tSV H ,anotisatiesvnt 
“ it @ i"titds edz egitseewnt 


Lawes a — th he ‘a <_< 
LBtisouvRTIa 0 Sans? sbiw sAg to 


2 tative vileauigen ,oldseetaisq sta cotdw 


xa teluotirsa at enol’salseovalt 


; tS fic mefdorq istavea ata, 
re%s saodT nottestvsevak 
6d biueds ~~ enottonut notssqihotinsaq LeboM ,f 
Oo emie7 ‘' suborq eddy gnivosisen Juensiw bavauisve 





‘ r ° : ’ 7 - va 
fiotdiw eentl Ichon goalvsed gebom sco!  ~n eTsbto tnetek? 5 


eas cds lliemeee aaa 
d2 diiw olstneonss selosto Mott viinsoltingte steh 
vi . 


genéqge0% to asbom {La vol smob ed 6 
_* - - | vy 
a : rae 
eon F (< 24 af. ‘ t b is t ag iw eit rf 





a bake 

a 

of | gia , ft. 
ne ‘oaks > a 








nates 


306 


2, The accuracy of the approximate solutions 
for these problems is dependent upon the degree of 
Satisfaction of edge conditions achieved using the boundary 
collocation technique. The accuracy of the solutions 
obtained in this thesis generally decreases with the 
number of boundary sides p, since the general solutions 
have been formulated in polar coordinates... However, errors 
in fundamental eigenvalues do increase with p for large p. 
Dnereflore wau snould prove beneficial to formulate solutions 
TOmunesce problems ein \Carvesian coordinates, especially for 
plates and shells having boundaries with 3 and 4 sides. In 
addition, the use of general solutions formulated in 
Cartesian rather_than polar coordinates should yield more 
readily approximate solutions for plates and shells having ~ 
rectangular boundaries with large aspect ratios and rhombic 
boundaries with small acute angles, This has been shown cM 
Che stavical case by Y..R. Kan [101] in the analysis of, 
uniformly loaded rectangular and rhombic plates, Kan 
obtained approximate solutions for these prob lens by 
applying the boundary collocation technique to the general 
solution of the biharmonic equation formulated in Cartesian 
coordinates. 

3. The effect of introducing time-dependent edge 
conditions should be investigated theoretically. It will 
be relatively easy to introduce the time-dependent edge 


conditions n(p,0,t) = 64(t) and F..(6,0,t) = f(t) 
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or €,,(6,6,t) = f(t) when either of the latter two edge 


~ 


conditions is associated with B68) = 0, since the 


Aah 
associated modal participation functions Xue A ist and 


wae are simply multiples of the zero-order terms of ier 
the modal participation functions related to the load 
distribution. The experimental results given for the 
elastically clamped square plate subjected to blast loading 
should provide a model for testing the practical 
applicability of introducing the time-dependent edge 
Conditions nip,e,t) = e'(t) and M (628s) ~ 

- B an(p,6,t)/av = e*(t), 

4, The theory for the dynamic response of plates 
should be extended to include the possibility of in-plane 
loading and membrane stresses, This will prove useful in 
practice even for cases where only the transient transverse 
bOsGoe a ee viOugiL GO be Significant, Since: in=plane 
loading will always occur to some degree at the edge of 
he Dlauc ac Une Dplave=will be elasticadly Clampea=tora 
mount which will be excited at its natural frequencies 
So wen, 

5. There are three areas related specifically 
to the dynamic response of shallow spherical shells which 
require further investigation, 

PIirstiy, Unere should be additional numerical 
studies in which the satisfaction of some of the various 
theoretical boundary conditions presented in this thesis 


are attempted in more detail. In particular, the effect 
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of varying the elastic clamping coefficients gt and gt 
Should be evaluated and additional effort should be made 
to obtain fundamental eigenvalues for shells with roller 
edges. Edge conditions in stress function which are the 
equivalent of the conditions (6,9) = u,,(5, 6) = 0 

and F (658) - piu (8,6) = 0, where gt is.an elastic 
clamping constant, should be formulated and satisfied 

for some numerical examples. 

Secondly, modal participation functions should 
be computed for loads distributed over annular rings 
Onevbnesshellesurface, so.that.different.load, time 
functions Q(t) can be applied simultaneously to different 
portions of the.shell.surface using, superposition. Then 
the, effect of variations in, time. of arrival. at..the 
ehelusslieeceroL a rectilinear, shock. front, for.example, 
can be studied and the maximum variations in load time 
functions over. the shell surface which can be ignored 
without introducing significant errors in the theoretical 
solution can be tabulated versus shell governing fre- 
quencies, 

Thipaly, pune,effect.of variations Jdnesheld 
geometry from the ideal, particularly variations in 
radius of curvature, should be studied both experi- 
mentally and theoretically in relation to values of 
the shell parameters a/R and a/h, 

6, Experimental values of viscous damping 


given in this thesis indicate the desirability of 





SOE 
2 


ka bas tg adnetotttees gniqmalo ottaals edt gatysay 20 
sbam ed divoda grotte Isnotstbbea bus bedsuleve ed biugda 
neffos ditw elfede tol eexisvaegios [ssnomabnut fitatdo of 
aid ots doldw aoltonw. esette al anotiibnos eyhd .éaghe 

0 = (6,4) uv = (0.3) » enoldtbnes saz to thelevinpe | 


v 
otinsie oa at by oTaiw 0 = (6,4) ape - (6,4) 9 baa 





helieliss brs besslumro? ed bivoda gasdenoco aniqmeip 
,eciqusaxe Lsoliemin omos 167 
biveds eaoltonut moltaqtoitisg Lebom ,yibnoose 


saatt taluane tevo betudtitgethb ebsol 102 basJugmos ag5 


i i 


emit beol dnetstith jad? coe ,sostive [feta sad fe 














inersttib o¢ yLevoensilumia bekiqgs sd mao (39 anoljonu? 
ted? .noktiseqisque gatey esatawe Llede edd 20 eaokisog, 
odt t8 Leviers to eats alt enotisiasy to tostie sig 
,eigusxe 10% ,¢nort Aoods tseatiigoss s to soataus Liens 
omits beol al enolisitev simixas edt bas beibude o¢ nao” 
stomyt ed aso dotdw sostuus [lede sdt steve enolsonut 
{soltetosdt eds at storie insoliingle gatouborJal duodtiw 
~stl antasevog Lleda auetev bossiuds? sd aso noteutos 
Lise at endttatrtey Lo Josette odd pe ea or 
ok aitottetiay: _ eee 

gee 

eh 


309 


including the possibility of nonlinear amplitude-dependent 
damping in the theoretical models, 

7. An elastic edge condition relating transverse 
displacement to Kirchoff's effective transverse shearing 
force should be formulated for both plates and shells. The 
applicability of this condition. can be tested by a theoretical 
comparison with the experimental results for the plate. 

8. A study should be conducted into the 
antisymmetric vibrations of these plates and shells, This 
can be done theoretically by replacing cos pn@ by sin pnod 
in the expressions for transverse displacement w and stress 
PUNCtI On wimancereplacing cos pm by sin pm@ an) the 
expression for the transverse loading function q. In this 
analysis, p will represent the number of antisymmetries in 
the load over the surface, while the characteristic segment 
for which the edge conditions must be satisfied numerically 
Woeowmeviweberdetined by 0 < @ < n/p. For example, pi = 2 
for a load which is antisymmetric about a shell plan 
diameter while p equals the number of sides of the regular 
polygonal boundary for a load which is antisymmetric about 
each of the radial lines passing through the midpoints of 


the boundary sides, 
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9. Predictions for the response of these plates 
and shells to central point loads using regular solutions 
do not satisfy the mathematical requirement of infinite 
transverse shear at the point load. Consequently, 
results using solutions which retain the singular terms 
should be compared to the present results. Inclusion of 
an inner boundary solution will be relatively easy as 


long as the outer boundary solution remains unchanged. 
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APPENDIX A 
SECTIONAL RESULTANTS 


A,1 Plate Sectional Resultants 
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Gime-denendent edge) condition Eqs. (36a) and (3,15b), 
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A.2 Shell Sectional Resultants 

The sectional resultants for shells subject to 
time-dependent edge condition Eqs. (3.15) and (3.16), 
ObvainecdsOyesupseicuGine Eqs. (38), (3.205), €3.27) and 
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APPENDIX B 


ORTHOGONALITY OF EIGENVECTORS 


yal Orthogonality of Plate Eigenvectors 
The orthogonality of the eigenvectors for plates, 
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If the motions represented by 1 (p,@) and 7 (p, 0) have 


different frequencies and y # y(p,6), then 
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proving orthogonality of the plate eigenvectors, 
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Substituting 
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[fen co, eas = AC 
Ss 


proving orthogonality of the shell transverse 


displacement eigenvectors, 
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APPENDIX C 


COMPUTER PROGRAMS 


C.1 Plate Program 


The plate digital computer program calculates 
the symmetric eigenvalues, eigenvectors, modal partici- 
pation functions for uniformly distributed and central 
point loads, transverse displacements and sectional 
resultants vs. time including viscous damping effects 
for regular polygonal plates subjected to various idealized 
transient load shapes and the elastically clamped 
homogeneous edge condition Eqs. (2.10a) and (2.,10b). 

The program has been written for use on the 
University of Alberta IBM 360/67 digital computer. Double 
precision (16 figure) accuracy is used in all computations. 
It has been modified for use on the DRES IBM 1130 digital 
computer on which it uses only ten figure computational 
accuracy. A high order computational accuracy and a large 
underflow-overflow exponent range ( 1077® on the IBM 360/67) 
are features required in order to invert accurately some of 
the larger boundary equation matrices which are ill- 
conditioned principally due to the presence of higher 
orders of the Bessel functions Jon and Ione 

Other computer programs which input 


experimental analogue data and output experimental and 
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theoretical strains, stresses, displacements and sectional 
resultants vs. time in the form of X-Y plots and determine 
theoretical mode shapes have been written for the DRES IBM 
1130 computer but are not included in this section. 

The plate program compilation from the Fortran 
Source deck,which uses approximately 52,000 core bytes, 
program listing and formation of an object deck require a 
total of 2.0 minutes execution time on the IBM 360/67. 
Program execution on this computer involving the search 
for one eigenvalue and the calculation of the response of 
the associated mode of vibration, excluding compilation 
Prmeeeregurres sopproximacveiy 0,5 minutes for aysolution 
using three collocation points on the characteristic 
segment boundary, 0.8 minutes for four collocation points, 
1.4 minutes for six collocation points and 4,8 minutes for 
benecollocation points. The largest portion of the 
computer time required for each of these solutions is used 
to solve the modal participation functions., 

A typical data set which causes the program to 
search for two symmetric modes of vibration of the partially 
clamped square plate discussed in Chapter 5 and compute the 
undamped modal response at two positions on the plate for 
an exponentially decaying blast load is listed by card 


number below: 
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A description of the data correlation to the 
general polygonal plate vibration problem is given in 
Table C.1. Problems of determinant value underflow 
(i.e. values Ais 0 on the IBM 360/67 or <1073° on the 
IBM 1130), are automatically corrected by the program 
which obtains a finite determinant value by incremen- 
ting the magnitude of the matrix element multiplier, 
data number {9a} by multiples of 10°, However, an 
initial choice of the matrix element multiplier which 
results in determinant value overflow (i.e. values >1076 
on the IBM 360/67 or >103° on the IBM 1130) will cause 
a computer ‘'hang-up', 

The program consists of a mainline routine 
which supervises the program execution according to 
Options specitied in. thei datvaysandeninessubroutines. 
Three of these subroutines DMINV, BESJ and BESI are 
modifications of standard IBM subroutines*« The principal 
functions of each of these routines are listed in Table 
C.2. A flow diagram of the possible program execution 
paths through the routines is shown in Fig. C.l. 

In order to facilitate data input and output 
on any computer, card input and line output variables 
termed IREAD and IPRIN, respectively, are defined in 
routine VIBRA and are used in all input/output Fortran 


statements, The computer line output is labelled 


appropriately with the exception that the values Ol to 
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Tabieu Cr. 


Explanation of Computer Data for Plates 


































Data 
Card 


Number Format 






Number of boundary sides p 
Density yg 


Circumscribing boundary radius 
Gali) 


Thickness h (in.) 
POLsSSson.s .racLo Vv 


Young's modulus E (psi) 





1,.E-50 for IBM 360/67; 
i beso LOC bia SU neCOnStant 
used for the backward recurrence 


series computation of J and I 
pn Pe 





Number of terms calculated in back- 
ward recurrence series, This number 
should be at least an order of 

10 greater than the highest order 
Bessel function required (i.e. 

Pye Gis) peas alas) 





i 


Number of boundary collocation 
points L 












F5.1,...} 9 locations of boundary céibocation 
points (degrees). Program is 
dimensioned for a maximum of 10 
collocation points on the charac- 


teristic segment boundary 


Dieta 
(t3a} 


numbers) 





Initial k/a value starting search 
for eigenvalues 


Final k/a value ending search for 
eigenvalues 
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Table C.1 (continued) 


Data Number Format 
Card 


Increment in k/a value 


Decimal place errorin eigenvalues, 
Kk, /a 


=l= clamped or partially clamped 
edge condition 
=2= simply supported edge condition 


Edge clamping factor 6 CDs) 


Viscous damping ratio of criticad 
damping 65/05 


Even number of boundary points used 
in, Simpson's formula Tor numerical 
invegration over 8 when calculating 
modal participation functions 


Maximum number of terms a used in 
calculating ®. and ae 


=l1= search for eigenvalues only and 
oli inee Cexsbarelcy USP yerey ASB. 5 a 

=2= use the sum of values {4a}? and 
{4b} as an eigenvalue and compute 
the associated eigenvector and a 
forced solution 

=3= search for eigenvalues, compute 
eigenvectors and forced solutions 


=O0= compute modal participation 
functions 7 
= any Value = use vonls Viawue aa xX, he 
and omit calculation of modal 
participation functions 


Controls the magnitude of the 
boundary equation matrix determinant 
value. Each element of the matrix 
is»multiplied by this value thereby 
eliminating some computer underflow 
and overflow. 





1 {4a} represents the first number on the fourth data card. 
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Table C.1 (continued) 


Data Number Format 
Card 


exponentially time- 
decaying load 

linearly time-decaying 
load 

flat impulse 

viscously damped structure 
Sub {ect svOba sti near Ly 
time-decaying load 
viscously damped structure 
subject. to an.exponentially 
time-decaying load 


Number of time instants for 
which transverse displacement 
and sectional resultants are 
Calculateduat each or 3G 
location 


Time interval (seconds) between 
calculation of transverse 
displacement and sectional 
resuilvants: aie e€achem.o 

location 


Positive duration of overpressure 
(seconds). Response cannot be 
calculated after this time 


Reference™overpressure qa. (psi) 


Number controlling output of 
transverse displacement and 
sectional resultants versus 
time to be on disk or paper, 
This feature was used on the 
IBM 1130 but is not shown in 
the program listing. Changes 
must be made in subroutine 
RESUL, Lorsits. use. 
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B45 


Table C.1 (continued) 


Format 


Number of r,6@ locations at 
which transverse displacement 
and sectional resultants are 
calculated 


Radius r off position of 
calculation: (in, ) 


Angle 6© of position of 
calculation (deg,.) 


Cards 10 to 14... are repeated 
once’ for eachsesymmetric mode 
of wibrationuptwhich will be 
found between the initial k/a 
value given by {4a} and the 
final k/a value given by {4b} 


=O0= program restarted with a 
new data set following 
=99= end of program execution, 


This card will be read when 
the k/a value exceeds the 
final k/a value given by {4b}. 
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Numbers in brackets { } refer to data as outlined in Table C.1. 


FIG. C,1 Flow Diagram of Plate Computer Program 
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O4 and 05 to O10 refer to the modal participation 
coefficients Je through ee and OF through oo» 
respectively... Calculation of the "sectional resultants “is 
inaccurate in the region near the plate center since the 
derivatives of the Bessel functions are calculated from 
their lower order derivatives, If values of sectional 
resulVanvomane desired in this region, it may be necessany 
to reduce the accuracy specified for the Bessel functions 
computed for subroutine RESUL by increasing the order of 
magnitude of the variable termed ERR in this subroutine, 


The Fortran source program is listed below: 
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C.2 Shell Program 

The shell digital computer program calculates the 
Symmetric eigenvalues, eigenvectors, modal participation 
functions foruniformly distributed and central point loads, 
transverse displacements and sectional resultants vs. time 
including viscous damping effects for shallow spherical 
shells with regular polygonal boundaries subjected to 
various idealized transient load shapes and the elastically 
clamped homogeneous edge condition Eqs. (3.22) and (3.24), 

The program has been written for the University of 
Alberta IBM 360/67 digital computer using double precision 
computation. The free vibration part of the program has 
been modified for use on the DRES IBM 1130 computer, Since 
the boundary equation elements are in dimensional form in | 
this» program, ill-conditioning results» from the presence of 
terms of the form r® as well as from the higher orders of 


the Bessel functions J and gin i 
pn pn 


Other programs used on the DRES IBM 1130 in the 
plate analysis and me eried in Appendix C.1 were modified 
for the shell analysis as well. In addition, a program 
calculating the functions e* and oe associated with the 
nonhomogeneous time-dependent edge conditions, Eqs. (seeia 
and (3.23),has been written for the IBM 1130, These 


programs are not included. 
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The shell program compilation from the Fortran 
Source deck, which uses approximately 84,000 core bytes, 
program listing and formation of an object deck require 
a total of 3.6 minutes execution time on the IBM 360/67. 
Program execution on the IBM 360/67 for the complete 
solution of one mode of vibration requires approximately 
one minute for a solution using three collocation points 
and two minutes for a solution using five collocation 
points. Accurate execution times for different solutions 
are not known since a solution using five collocation 
points took as long as 4.25 minutes when the computer 
monitor was at level 13 and as little as 1.15 minutes when 
the monitor was at level 16, 

A typical data set which causes the program to 
search for one symmetric mode of vibration of a simply 
Supported shell with 12 sides and solve the undamped modal 
response at three positions on the shell surface for an 
exponentially decaying blast load is listed by card 


number below: 
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Data 
Card 
1 12 97.5E-323267,2E=3 37.5E-2 3.3E-1 1 On+7 64,0E+0 
S 1,E-50 50 
5) eran ems. 0. 15.0 
4 0O,1E-1 3.,0E-1 0,2E-1 0,.5E-6 
5 6 Mel0B-10 0 ,0E+00 0,0E+00 
6 3652020 
fi 3 
8 0.,0E+0 
9 1.0E+0 
10 re 50 0,2E-4 30,0E-3+10.0 
nal af 
Ee 3 
H3 0.000 0.000 
14 5.000 0.000 
5 22.476 0.000 
16 99 


The data correlation to the general shell vibration 
problem follows the explanation given for plate data in 
Table C.1 except for the modifications which are given in 
Table Ces reas x 
oF (6, 9) oF (p, 8) 

The edge conditilonse———=— os/0yand/s7-—~., 6=70 
OV ds 


which are employed when data {5a} = 3, 4, 5 or 7, were 


postulated by Oniashvili [88] to represent vanishing 
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in-plane displacements perpendicular and tangential to the 
Shell edge, respectively. Satisfaction of these edge 
conditions did not yield satisfactory fundamental frequencies 


for clamped shells with 12 and 15-sided polygonal boundaries. 


Tabtesects 


Modifications of Plate Computer Data for Shells 


| Data Number Format 
Card 


Radius of curvature of shell middle 
surface R (in.) 


Number of terms calculated in 
backward recurrence series; 
should be 2 p(L-1) + 11 


bF.ee ORSSL 8 locations of boundary collocation 
({3a} points (deg.). Program is 
numbers) dimensioned for L < 7 


a i2 Number designating shell homogeneous 
edge conditions to be satisfied; 
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Table C.3 (continued) 


Number Format 


Edge clamping factor RB bey ime. 


=l= enter subroutine NONDIS and 
search for a determinant value 
which becomes progressively 
smaller in magnitude but does not 
change sign 


=0= do not enter subroutine 
NONDIS 





The shell program mainline routine and subroutines 
perform the same functions as the plate routines with the 
Same names as outlined in Table C.2 and the shell program 
execution essentially follows the flow diagram for plates 
shown in Fig. C.1l. In order to accommodate the shell 
program on the IBM 360/67 it was necessary to divide the 
Subroutine BLOND into five shorter subroutines called 
BLOND, BLOND2, BLOND3, BLOND4, and BLOND5,and the subroutine 
FOURI into two shorter subroutines called FOURI and FOURI2. 
One additional subroutine NONDIS was written for the 
shell program, This routine is entered only when data 
{7b} = 1. It searches for eigenvalues for which the 
determinant value magnitude becomes progressively smaller 
but does not change in sign. This routine was used in an 


unsuccessful attempt to find eigenvalues for clamped shells 















e 
ti _ 


=e 
Pa nl Lad 
(beun} noo) F.o aildsT 


Po et ee oe ee 6 RE INR oy Demet ttn epee 
: $amt07% to dmul 


“y - - te ote etn Sennen enemies ee 


£.0L8 b 


— 


) 
ele 3 S d 

; : moO resi 4 
r a? ; 
4 : 
ifane } 
- ad 4 

rs bila’ 
: 
m : 
Cc ve | 

: * vf 

he 4 =i 


a ee me rag A tl aba 


¥vilstinsees nolivosexs - 


cebio.aI «41.9 ,glt af mwods 


. 
No metRoX . 


~~ 
Lee) 
i! 
{ 
oa 
< 


i304 iettode evil otat TWOJE entsvouGusS 4 a 





entsvoudyve ot? bos .cI40U8 bas Mole ,€QvOTa ,SdvOJE ,atoOle , 
.SIHUOG bone TRUCT belles sentiuordue tstiona ow? ognt TAUOT | 
oft (92 deszinow asw 2TQMOM saidvordus fanoljibbe sn0 
; ee : : .) - F 
: atdT , ms vie ‘ff ii i 7 
; — 5 : Aw sf o 
_-  wa- : 


a¢ ih) 8 3 al i 7 
i wo 









‘pied nentw ylne betesme at enLuo 
Came 





: 60% dotov 102 easeisvasgte 20% 
Pon) _ 






PRE ee 





a 


nofiame ylovissstacrg 


7 
2t 
“J 


a. -_* a3 
ve ; - 


369 


using the edge conditions of vanishing in-plane dis- 
placements as postulated by Oniashvili. 

The shell computer output is labelled appro- 
priately with the exception that the values of O1 to 
012 and 013 to 018 refer to the modal participation 
coefficients %, through a and >, through oo 
respectively. 

The Fortran source program, except for the 
subroutines BESJ, BESI and SIMP which are unchanged 


from the plate program, is listed below: 
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